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ABSTRACT 


A  theoretical  study  is  given  for  store-and-forvard  communication 
networks  in  which  the  nodes  have  finite  storage  capacity  for  messages. 

A  node  is  "blocked"  when  its  storage  is  filled,  otherwise  it  is  "free." 

A  two-state  Markov  model  is  proposed  for  each  node,  and  the  fraction  of 
blocked  nodes  in  the  network  is  shown  also  to  have  a  two-state  Markov 
process  representation.  The  tiia^-dependent  probability  that  any  given 
node  in  the  '-.etwork  is  blocked  is  obtained  for  some  uniform  networks  of 
arbitrary  dimension,  and  various  results  describe  the  clumping  phenomena 
In  these  networks. 

Through  a  modification  of  the  basic  Markovian  network  model,  the 
fraction  of  blocked  nodes  in  a  computer-simulated  store-and-forward 
communication  network  is  predicted  with  reasonable  accuracy. 
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CHAPTER  1 


INTRODUCTION 


A.  Computer  Networks 

In  the  early  1960's  the  first  tune-sharing  facility  began  opera¬ 
tion.  Since  that  tire,  facilities  and  systems  have  grown  and  developed 
across  the  country  into  quite  sophisticated  and  unique  sites  each  having 
special  features  and  capabilities  in  the  form  of  exceptional  computer 
prograne,  data  files,  hardware  devices,  resources,  and  human  talent 
which,  in  general,  are  not  easily  transferable.  A  desire  to  share  these 
resources  has  led  to  the  development  of  computer  networks  which  permit 
the  separate  ocrputer  facilities  to  communicate  with  each  other. 

A  computer  network  is  a  collection  of  nodes  (computers)  connected 
together  by  a  set  of  links  or  lines  (cxxrmunicaticn  channels) .  Messages 
in  the  form  of  commands,  inquiries,  replies,  and  file  transmissions 
travel  through  this  network  over  data  transmission  lines.  At  the  nodes, 
the  task  of  relaying  messages  (with  all  proper  routing,  acknowledging, 
error  control,  queueing,  etc.)  and  inserting  and  removing  messages  which 
originate  and  terminate  at  that  node  must  be  carried  out. 

The  Advanced  Research  Projects  Agency  (ARPA)  Network  [1-5]  is  a 
store-and-forward  computer  communication  network  linking  approximately 
fifteen  research  centers  across  the  country  at  the  present  time  with 
approximately  five  more  scheduled  for  completion  by  the  end  of  1971.  In 
a  store-and-forward  network,  messages  are  broken  up  into  convenient 
sized  packets  that  individually  make  their  way  through  the  net,  "hopping" 
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from  node  to  node.  If  a  packet  cannot  be  transmitted  iirmediately  out  of 
a  node  on  its  way  through  the  net  because  its  designated  output  line  is 
in  use,  it  forms  a  queue  and  awaits  its  turn  to  be  transmitted. 

Western  Union  has  used  the  store-and-forward  concept  for  years  as 
has  the  United  States  Air  FOrce  in  its  Sage  Defense  System.  In  November 
of  1969  DATRAN  Corporation  proposed  to  the  Federal  Ccrruni cations  Comas- 
sicn  a  network  for  digital  ocmtuni cation  linking  35  metropolitan  areas 
from'Boston  to  San  Francisco  and  conp rising  240  tnicrowaye  relay  stations. 
Eventually  they  propose  to  make  it  a  store-and-forward  network  [6] .  We 
see  that  numerous  store-and-forward  networks  are  already  in  use  and 
others  ^re  being  planned.  ] 

B.  Structure  of  the  AHPA  Network 

'  Let  us  examine  the  structure  of  the  ARPA  Network  more  carefully. 

\ 

At  each  site  in  the  network  there  is  at  least  one  large  digital  computer 

called  a  HOST,  which  acts  as  a  source  and  terminal  for  messages  in\  the. 

\  1  ! 

network.  These  computers  are  basically  inoompatible  in  hardware,  soft¬ 
ware,  file  structure,  etc.,  and  hence  there  is  a  need  for  an  intermedi¬ 
ate  device  to  interface  these  HOSTs  to  the  communication  net  which  con¬ 
nects  them.  This  function,  and  others,  is  performed  in  the  ARPA  Network 
at  each  site  by  a  digital  computer  called  an  Interface  Message  Processor 
(IMP) .  The  IMPs  carry  out  the  message  handling  process  in  the  network, 

so  when  we  speak  of  the  nodes  in  the  net  we  are  actually  referring  to 

1  _  \ 

the  various  IM’s.  >  ' 

i 

An  IM?  in  the  ARPA  net  receives  messages  from  two!  sources: 

1.  Other  BPs  like  itself  over  fully  duplex  50  Kbit/sec.  leased 

\  1 

telephone  lines.  ) 
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2.  One  or  more  HOSTs  over  100  Kbit/sec.  fully  duplex  lines. 
Message  bits  are  sent  in  series  and  are  protected  by  error  detection 
schemes.  If  an  error  is  detected,  the  message  must  be  retransmitted. 
Compared  to  any  HOST  computer,  tne  IMP  is  a  small  machine  with 

i 

finite  storage  space  for  messages.  _  'Part  of  the  IMP  storage  is  strictly 
allocated  for  messages  which  are  relayed  from  neighboring  IMPs  and  which 
must  be  transmitted  to  still  another  IMP  before  reaching  their  destina¬ 
tion;  this  is  called  store- and- forward  traffic.  Part  of  the  remining 

•  i 

i  i 

storage  in  an  is  strictly  allocated  for  the  reassembly  of  nulti- 

'  .  \ 

packet  messages  destined  for  one  of  the  DP's  HOSTs.  ’(A  multi-packet 

message  is  one  which  is  too  large  to  be  transmitted  as  a  single  packet 

\ 

whose  maximum  size  is  1008  blits.  Multi -packet  messages  may  be  ip  to  8 

! 

packets  in  length,  and  each  of  these  packets  must  be  held  until  all  are 

received  in  the  final  node,  at  which 

\ 

original  message  and  delivered  to  the 
partitioned  in  the  HOST  into  many  multi-packet  messages.)  The  remaining 

storage  is  allocated  between  these  two  types 

\ 

all,  the  IM5  contains  storage  space  for  about  50  single-packet  messages. 


From  time  to  time,  during  periods  of  high  utilization,  the  IMP's 
storage  can  become  filled,  so  that  arriving  messages  must  be  refused. 

i 

then  this  occurs  we  say  that  the  node  is  "blocked."  Blocking  in  the  IMP 
can  occur  in  any  of  three  ways:  , 

1.  There  are  no  more  reassembly  spaces  available  for  HOST  traffic, 
and  packets  for  a  HOST  that  were  sent  by  other  IM’s  must  be  refused. 

2.  There  are  no  more  spaces  available  for  s  tore- and- forward 

\ 


of  traffic  as  needed.  In 


time  they  are  reassembled  into  the 
HOST.  Longer  messages  must  be 


1 


3 


traffic,  and  this  ncn-HOST  packets  must  be  refused. 

3.  Biere  are  no  more  spaces  for  arriving  messages  and  all  traffic 

must  be  refused. 

Certain  high  priority  messages  are  never  blocked,  e.g.,  space  is  always 
saved  for  positive  acknowledgments  sent  by  neighboring  IMPs  to  indicate 
that  a  message  previously  sent  by  the  IMP  has  been  received  without  er¬ 
ror  and  can  thus  be  discarded  by  the  IMP.  the  other  hand,  a  blocked 
message  is  ignored  by  the  IMP,  and  the  absence  of  a  positive  acknowledg¬ 
ment  tells  the  Bf>  which  sent  the  message  that  the  message  will  require 
retransmission . 

Selective  blocking,  as  in  points  (1)  and  (2)  above,  or  total  block¬ 
ing,  as  in  (3) ,  can  occur  in  this  network  if  the  input  rate  of  messages 
equals  or  exceeds  the  output  capacity  over  a  period  of  time.  We  would 
normally  expect  this  to  occur  only  during  peak  hours  of  the  day.  How¬ 
ever,  it  is  a  potentially  dangerous  situation  because  a  blocked  neighbor 
reduces  a  node's  message  output  rate  with  no  corresponding  change  in  its 
input  rate.  This  causes  its  storage  to  fill  at  a  faster  rate  and  in¬ 
creases  its  chance  of  becoming  blocked.  Thus  blocking  could  propagate 
in  both  space  and  time. 

The  purpose  of  this  research  is  to  gain  an  understanding  of  the 
blocking  behavior  in  a  message-switching  network. 
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CHAPTER  2 


THE  MODEL 


A.  General  Description 

Selective  blocking  is  a  very  difficult  problem  to  analyze.  The 
allocation  of  storage  between  store-and-forward  traffic  and  HOST  traffic 
is  equivalent  to  the  formation  of  two  distinct  queues  with  finite  wait¬ 
ing  room,  or  storage  space,  in  which  the  maximum  size  of  the  waiting 
room  for  each  queue  is  dependent  on  the  nunber  of  customers  (i.e.  mes¬ 
sages)  in  the  other  queue.  To  make  the  problem  mathematically  tractable, 
the  network  we  analyze  will  consist  of  nodes  having  a  single  queue  for 
messages.  If  there  is  an  empty  space  in  the  queue,  the  first  arriving 
message,  regardless  of  its  final  destination,  will  take  that  space.  If 
there  are  no  spaces  for  arriving  messages,  then  the  node  is  "blocked." 

As  soon  as  one  message  is  transmitted  by  a  blocked  node,  it  becomes 
a  "free"  node.  It  remains  in  this  state  as  long  as  there  is  at  least 
one  empty  spaae  in  storage  that  could  be  used  by  an  arriving  message. 
When  the  storage  fills  again,  the  node  re-enters  the  blocked  state. 

Figure  1  shows  a  simplified  model  of  such  a  node  in  the  terminology 
of  the  ARPA  Network.  The  IM5,  when  free,  accepts  messages  into  its  main 
storage  from  two  sources: 

1.  Other  IM>s . 

2.  A  single  HOST  which  generates  and  receives  messages  (as  a 
source  and  terminal) . 

A  message  in  a  message  buffer  is  queued  vp  for  transmissicn  over  an 
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FROM  HOST 

Figure  1.  Schematic  of  a  Node 


TO  HOST 


appropriate  output  line  to  some  neighbor  as  determined  by  the  final  des¬ 
tination  of  the  message,  and  is  then  transmitted  serially  to  that  neigh¬ 
bor.  Any  of  these  neighbors  can  become  blocked,  thus  preventing  the  use 
of  the  output  line  feeding  such  neighbors. 

Nodal  blocking  is  caused  by  the  finite  storage  room  for  messages  in 
the  IMP  and  the  overutilization  of  the  system.  By  overutilization,  we 
mean  that  when  the  node  is  accepting  messages,  its  average  arrival  rate 
equals  or  exceeds  its  average  service  rate  (which  is  the  total  output 
channel  capacity  divided  by  the  average  message  length) .  Elementary 
queueing  theory  [7]  shows  that  if  (1)  the  system  is  underutilized,  and 
(2)  there  is  storage  space  for  approximately  twenty  messages  or  more, 
then  under  fairly  general  conditions  there  will  be  essentially  no 
blocking. 

The  analysis  of  the  propagation  of  blocking  is  difficult  for  at 
least  three  reasons.  First,  it  involves  networks  of  queues  for  which 
only  stationary  results  at  best  can  generally  be  obtained.  Second,  the 
pertinent  stochastic  processes  are  dependent,  for  if  a  node  becomes 
blocked,  it  cannot  accept  messages  from  its  nei^bors  and  cheir  storage 
will  tend  to  fill  at  a  faster  rate.  Finally,  it  is  a  transient  queueing 
problem  and  even  the  sinplest  of  these  is  very  difficult  to  solve.  (For 
exanple,  the  queueing  system  with  Markovian  arrivals,  a  single  experien¬ 
tial  server,  and  unlimited  waiting  room  has  modified  Bessel  functions  in 
its  time  dependent  solution  [7].) 

B.  Related  Work 

A  nunber  of  topics  in  graph  theory  are  related  to  this  problem. 
Ignition  phenomena  as  developed  by  Rapoport  [8]  and  Allans  on  [9]  treats 
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vertices  (nodes)  which  are  "excited"  if  they  receive  a  certain  minimum 
nutber  of  stimuli  within  a  certain  amount  of  time.  This  excitation  is 
assumed  to  stimulate  d  other  vertices  to  which  it  is  randomly  connected. 
Stable  states#  i.e. ,  constant  fractions  of  vertices  being  excited,  are 
shown  to  exist  under  seme  conditions.  Ihis  model  has  immediate  applica¬ 
tion  to  neural  networks  because  of  their  essentially  random  connectivity 
and  the  nearly  deterministic  behavior  of  neurons.  However,  the  model 
cannot  be  reasonably  applied  to  conputer  networks  because  they  are  not 
randomly  connected  and  the  probabilistic  nature  of  information  transfer 
in  the  form  of  variable  length  and  time  of  arrival  of  messages  makes  the 
excitation  process  (i.e.  the  blocking)  very  much  non -deterministic. 

Percolation  Theory  [10]  considers  lattices  in  which  a  branch  be¬ 
tween  any  two  nodes  is  present  with  probability  p  or  deleted  with 
probability  1  -  p.  The  main  concern  here  is  the  minimum  value  of  p 
(i.e.  the  critical  value)  for  which  a  connected  component  of  infinite 
length  exists  in  the  lattice  with  probability  one.  The  relation  of  this 
theory  to  the  work  of  Gilbert  [11]  an  random  plane  networks  is  clear. 

In  the  study  of  probabilistic  graphs  [12] ,  branches  anchor  vertices 
are  deleted  in  some  random  fashion.  The  questions  raised  (and  answered) 
are  the  following:  what  are  the  probabilities  corresponding  to  various 
kinds  of  connectivity;  what  is  the  distribution  of  the  size  of  the  larg¬ 
est  connected  component,  etc. 

An  interesting  variation  on  the  network  vulnerability  problem  is 
that  which  considers  a  probabilistic  repair  time  for  vertices  or 
brandies  that  have  been  damaged  by  an  attack  from  scare  weapon  system. 

In  [13]  the  tine  varying  probability  of  connectivity  is  determined  for 
random  graphs. 
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In  none  of  these  areas  of  graph  theory  is  the  state  of  a  node  (or 
vertex)  ever  taken  to  be  a  function  of  the  states  of  its  neighbors, 
lhus  such  results  are  not  applicable  to  the  stucV  of  blocking  propaga¬ 
tion. 

Eden  (14]  and  Morgan  and  Welsh  [15]  studied  two-dimensional  Poisson 
growth  processes.  They  assured  that  "infection"  in  a  cell  network 
spread  from  cell  to  neighboring  cell  in  an  amount  of  time  taken  from 
some  probability  distribution.  These  authors  obtained  results  on  the 
shape  of  the  infected  area  and  the  rate  of  spread  of  the  infection.  In 
their  models,  once  a  cell  becomes  infected  it  remains  in  that  state  for¬ 
ever,  thus  their  work  cannot  be  taken  as  a  solution  to  the  blocking 
problem. 

Roach  [16]  studied  the  overlap  of  objects  placed  at  random  in  sane 
space  and  called  these  overlappings  "clunps."  He  treats  the  nunber  of 
clurps,  their  size,  their  shape,  and  the  spacing  between  them  for  a  num¬ 
ber  of  interesting  cases,  including  the  square  lattice.  He  assumes  the 
probability  that  a  lattice  point  is  marked  (i.e.  blocked)  is  the  same 
for  all  lattice  points.  Because  of  this  independence  assumption  and 
also  because  his  system  is  static,  we  cannot  utilize  his  results;  how¬ 
ever,  we  will  adopt  his  terminology. 

C.  The  Mathematical  Model 

The  blocking  problem  is  a  difficult  one.  Since  we  cannot  solve  the 
problem  exactly,  our  goal  is  to  make  good  approximations  that  allow  us  to 
analyze  the  system  and  characterize  its  blocking  behavior  in  seme  way. 

To  this  end  we  make  the  following  assumptions: 
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1. 


The  HOST  cannot  become  blocked  (it  is  an  infinite  sink) . 

2.  a.  Input  traffic  from  the  HOST  is  Poisson. 

b.  Traffic  on  all  lines  (including  the  HOST-IMP  line)  has 
the  same  average  rate  so  that  total  traffic  into  each 
node  is  o  messages/sec. 

3.  a.  Message  lengths  are  exponentially  distributed. 

b.  Service  (transmission)  time  on  any  line  is  therefore 
exponentially  distributed  such  that  for  a  node  with  k 
blocked  nei^ibors,  the  rate  at  which  messages  exit  from 
that  node  is  messages/sec.,  dependent  on  the  num¬ 

ber  of  blocked  neighbors. 

4.  The  probability  of  an  eitpty  queue  in  the  IMP  is  approximately 
zero  (since  the  system  is  assured  to  be  overutilized. 
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CHAPTER  3 


ANALYSIS 

A.  The  Nodal  Model 

Under  the  assumptions  in  "The  Mathematical  Model"  (Section  2.C) ,  we 
arrive  at  a  simplified  blocking  model  for  a  node  in  the  network  as  a 
two-state  Markov  process  (Fig.  2).  If  the  node  is  blocked,  i.e. ,  in 
state  b,  it  becomes  free  in  the  next  instant  of  time  At  with  proba¬ 
bility  v  ^  At  where  k  is  the  number  of  blocked  neighbors  it  is 
experiencing  at  that  time.  Similarly,  if  the  node  is  free,  i.e.,  in 
state  f ,  it  becomes  blocked  in  the  next  instant  of  time  At  with 
probability  xHt  where  k  is  again  the  nunber  of  blocked  nei^rbors. 
Thus  X is  the  rate  at  which  a  free  node  becomes  blocked  in  the 

(k) 

presence  of  k  blocked  neighbors ,  and  snculd  increase  with  k.  Mv  , 
on  the  other  hand,  being  the  rate  at  which  a  blocked  node  becomes  free, 
should  decrease  with  k. 

1.  Derivation  of  y  ^ 

Below  we  shew  the  appropriateness  of  this  model.  First,  we 
require  the  Laplace  transform  of  the  message  interdeparture  time  proba¬ 
bility  density  =D(s).  For  any  node  let  p  =  P  [non-enpty  node]  and  let 
the  Laplace  transform  of  the  probability  density  of  the  message  inter¬ 
arrival  time  process  be  A(s) .  Because  we  have  assumed  that  the  service 
time  is  exponential  with  parameter  y  ^ ,  we  knew  that  the  Laplace 
transform  of  the  departure  process,  conditioned  on  a  non-empty  system  is 
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Figure  2.  Blocking  Modal  for  an  Imp 
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b)  DUAL  QUEUE  STATE  TRANSITIONS 

Figure  3 
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pW/(s  +  vi  ) .  Therefore, 


01J(k) 

D(s)  =  Py . ,Lr  +  (1  -  P)  A(s) 


u 


(k) 


(1) 


By  assumption  (4)  we  have  p  *•  1 


D(S) 


,  (k) 


w 


s  +  u 


(2) 


which  says  that  the  departure  process  is  a  Poisson  stream.  See  Burke 
[17]  and  Reich  [18]  for  further  details  on  departure  processes. 

We  have  assumed  that  the  traffic  on  all  lines  has  the  same 
average  rate.  If,  for  exanple,  every  node  has  exactly  four  neighbors 
and  one  HOST,  then  there  are  five  output  lines  from  each  node.  All  of 
these  lines  are  equivalent  (except  that  the  HOST  cannot  become  blocked) 
and,  by  the  assumption  of  exponential  message  lengths,  the  departure 
process  from  each  output  line  constitutes  a  Poisson  stream  at  rate 
Vi  /5  when  that  neighbor  is  not  blocked  {and  at  rate  0  when  that 
neighbor  is  blocked) . 

y(k)  =  IzJL  y(0)  k  =  0,1,. ..,4  (3) 

5 

where  y  is  a  given  system  parameter  and  represents  the  maximum  mes¬ 
sage  departure  rate  from  a  node.  This  set  of  numbers  is  merely  an  illus 
traticn;  any  conbination  can  be  treated  by  this  model.  These  results 
show  that  we  can  approximate  the  time  spent  in  the  blocked  state  as  be¬ 
ing  exponentially  distributed  with  parameter  y  ^  •  Because  of  the 
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nenoryless  property  of  the  exponential  distribution,  the  expected  value 
of  the  remaining  time  to  be  spent  in  the  blocked  state,  given  that  k 
changes  to  seme  new  value  k^  while  in  the  blocked  state,  is  siirply 
1/M(kn).  By  Eq.  (3)  this  means  that  an  increase  in  k  should  tend  to 
increase  the  time  spent  in  the  blocked  state,  and  a  decrease  in  k 
should  tend  to  decrease  this  time.  We  would  expect  to  see  such  behavior 
in  a  reed  computer  network. 

2.  Derivation  of 

Hie  derivation  of  the  parameter  is  not  nearly  as  sinple. 

The  time  that  an  IMP  spends  in  the  free  state  is  distributed  as  the  busy 
period  in  a  queueing  system  with  finite  queueing  rocm  for  customers,  as 
we  now  show.  We  begin  by  first  considering  the  state  transition  diagram 
or  Markov  chain  model  for  such  a  single  node  finite  storage  queueing 
system  as  shown  in  Figure  3a.  The  nunbers  inside  the  circles  represent 
the  nutber  of  customers  (messages)  in  the  node.  We  assume  that  custo¬ 
mers  arrive  in  a  Poisson  fashion  with  parameter  o,  and  depart  after 
receiving  service  (exponentially  distributed  with  an  average  of  1/P 
seconds) .  A  busy  period  begins  when  a  customer  arrives  to  find  an  enpty 
system  (at  which  time  he  immediately  enters  the  service  facility) . 
Customers  arriving  during  his  service  time  form  a  queue  behind  him. 

With  each  arrival  the  system  moves  to  the  right  along  the  state  transi¬ 
tion  diagram  because  the  nuirber  in  the  system  is  increased  by  one ,  and 
with  vjh  service  oonpleticn  (i.e.,  departure)  it  moves  to  the  left. 
Customers  arriving  when  the  system  contains  N  customers  are  lost 
(i.e.,  depart  without  service).  The  busy  period  ends  the  first  time  the 
system  goes  enpty  after  initiation  of  the  busy  period. 
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For  the  IMP  model  we  now  consider  a  dual  queue  in  which  the 
roles  of  service  and  arrival  are  reversed,  and  the  mitbers  inside  the 
ti-i  oq  now  represent  the  nvirber  of  enpty  places  in  storage  that  could 
be  used  by  arriving  messages  (Fig.  3b) .  The  free  period  of  the  IMP 
begins  with  the  departure  of  a  message  from  a  previously  filled  system, 
i.e. ,  no  enpty  places  for  arriving  messages.  With  a  transmission  (de¬ 
parture)  the  system  moves  from  state  0  to  state  1.  It  continues  to  move 
to  the  right  with  each  transmission  and  to  the  left  with  each  arrival. 
Hie  free  period  ends  the  first  time  the  system  returns  to  the  0  state. 
Die  correspondence  between  the  primal  and  dual  queues  is  perfect;  thus 
any  results  obtained  for  the  busy  period  in  the  primal  system  are  appli¬ 
cable  to  the  dual  queue  free  period  in  the  IMP  sinply  by  substituting 
yM  for  o  and  a  for  y,  as  in  Figs.  3a,b. 

The  busy  period  for  a  finite  queueing  room  system  is  difficult 
to  obtain,  but  the  result  for  unlimited  queueing  room  is  well  kncwn. 

The  probability  density  of  the  length  t  of  the  busy  period  in  such  a 
system  is 

p(t)  ^  _L.e-(a+y)tii(2t Jov)  (4) 

tSfi 

where  p,  the  utilization  factor  =  (o/y)  <  1  and  I^(x)  is  the  modi¬ 
fied  Bessel  function  of  the  first  kind,  of  order  one  [19] .  If  the  size 
of  the  queueing  room  is  greater  than  20,  the  solution  for  unlimited 
queueing  room  is  a  good  approximate  solution  to  the  limited  queueing 
room  problem.  (This  follows  since  we  have  assumed  P  [enpty  IMP]  =  0; 
but  the  P  [enpty  IMP]  corresponds  to  the  probability  of  being  in  state  N 
(i.e. ,  all  N  spaces  are  enpty)  in  Fig.  3b,  and  thus  an  increase  in  N 
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will  not  seriously  affect  our  results.)  We  make  the  further  approxima¬ 
tion  that  Eq.  (4)  holds  when  a  varies  as  p^ ,  i.e.,  when  o  is  time 
varying.  Since  we  have  assured  overuti li zation ,  we  have  (vj^/o)  <  1, 
and  we  are  justified  in  substituting  this  (or  p  /o)  for  p.  Thus  we 
get  the  following  for  the  approximate  probability  density  of  the  length 
t  of  the  time  spent  in  the  free  state:  \ 

p  <t}  "  e~ <0+P  °C)  )txi  (2tVoy(k) )  (5) 


As  the  ratio  p /a  approaches  0,  i.e.,  as  the  system  becomes  more 

\ 

overutilized,  this  density  approaches  that  of  the  exponential  distribu¬ 
tion  except  out  on  the  tail  of  the  distribution  where  the  probability 

' 

density  will  be  assured  negligible.  To  arrive  at  ;  a  more  tractable 
model,  we  therefore,  approximate  the  free  period  distribution  by  an  expo¬ 
nential  distribution  having  the  same  meari  value.  The  mean  value  of  the 
busy  period  in  the  original  system  is  easy  to  obtain/  and  is  given  by 
l/p(l  -  p) .  Therefore ,  as  an  approximation  to  the  free  period  in  the 

»  i  . 

IMP,  we  take  an  exponential  distribution  with  mean  value  1 /(a  -  p ^ ) , 

*{k)  =  a  -  pW  (6) 


For  the  marginal  case,  0  ■  \p'u' ,  elementary  queueing  theory  shows  that 

i 

we  must  take 


o/N  for  a  =  p 


where  N  is  the  size  of  the  storage  capacity  of  the  IMP  (in  messages) , 


\ 


I 


|  Our  model  for  the  blocking  IMP  is  thus  a  two-state  Markov 
process  or,  in  the  language  of  renewed,  theory,  an  alternating  Poisson 
renewed  process  [20] ... 


B.  Derivation  of  the  Network  Model  , 

One  way  to  describe  the  dynamics  of  a  network  of  such  nodes  is  to 
examine  the  probability  that  any  given  node  is  blocked  at  sane  time  t. 

Fbr  a  network  let  us  enplqy  a  two-dimensional  integer  lattice.  In  this  \ 
way  we  can  have  a  large  system  and  yet  minimize  the  aorplexity  of  its 
description.  Consider  a  node  with  *L  four  neighbors  nunbered  1  to  4: 


and  let 


pk(t)  =  p[k  neighbors  blocked  at  time  t] 


p(t)  *  P[nyde  blocked  at  time  t]  (9) 

Then,  from  elerrentary  considerations,  we  have  (correct  to  within  o(At)) 

/ 

^  *4  \ 

p(t  +  it)  -  (1  -  p(«)g)Pk(t)Xtk,it  +  p(t)  (1  -  |>k(t)u(k)At) 


where  from  Bq.  (3) 


^(k)  =  v(0)  _  (k/SJy* 
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and  from  Eq.  (6) 


X(k)=  o  -  y(k)  =  o  -  y(0)  +  (k/5)y(0) 

for  o  >  y^.  We  will  assume  that  this  holds  for  a  =  y*°*  as  well. 
Uie  usefulness  of  the  results  that  we  will  obtain  will  justify  this 
approximation. 

We  also  note  that 

X(k)  +  y<k)  =o  (10) 


Thus, 

P(t  ♦  At)  -  P(t), ,  (i  -p(t))X;pk(t)x(k)  -  P(t)fpk(t)lJ(k) 

At  iS)  k=0 

Letting  At  approach  0,  we  have 

-P<t>y;pk<t><x<1'>  ♦  M(k))  +  Ei*<m(k) 
at  k=0  k=0 

=  -<jp(t)  ^(t)  +  EPk(t)  (o  -  y(0)  +  (k/5)y (0)) 
k=0  k=0 

fn\  . , (0)  4  i_ 

=  -op(t)  +  a  -  y  +  ^  kP  (t)  (11) 

3  k=0 

This  can  be  sinplified  by  noting  that 

ETnurber  of  blocked  neighbors  at  time  t]  =  £  kP^t)  (12) 

k=0 

where  E  denotes  expectation.  Define  the  indicator  function 
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II  if  node  n  is  blocked  at  tine  t 
0  otherwise 

Now  let 

p  (t)  *  P[node  n  is  blocked  at  tine  t] 
n 

then 

E[f  (t)]  =  p  (t)  (13) 

n  n 

Further,  from  Eq.  (12)  we  have  that 

£  kP^t)  =  E(£  f  (t))  =  £  E(fn(t))  (14) 

k=0  neM  neM 

where  M  is  the  set  of  neighbors  for  this  node  (wliich  we  nuttoer  1,2,3, 
4).  From  Eqs.  (13)  and  (14)  we  get 

4 

SkP^t)  =  p,(t)  +  p2(t)  +  p3(t)  +  p4(t)  (15) 

k»0 

Finally,  from  Eqs.  (11)  and  (15)  we  have  the  result 

-op(t)  +o  -  y(0)  +  4— (P,  (t)  +  p2(t)  +  p3(t)  +  p4(t))  (16) 

at  3  A 

It  is  interesting  that  this  relation  can  also  be  derived  from 
epidemiology.  We  will  adopt  the  notation  from  Bartlett  [21] . 

Consider  a  deterministic  epidemic  without  migration  of  individuals 
and  with  but  two  types  of  individuals,  infected  and  susceptible,  in 
which  "cured"  individuals  are  returned  to  the  susceptible  ranks.  Assure 
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that  the  number  of  individuals  in  an  infected  group  who  are  cured  at 
fjmo  t  +  At  is  equal  to  the  mxrber  of  infectives  in  the  group  at  tine 
t  multiplied  by  a  constant,  uQ,  diminished  by  the  nuntoer  of  infec¬ 
tives  found  simultaneously  in  surrounding  areas  weighted  in  sate  spatial 
manner.  Similarly,  we  will  assume  that  the  nvirber  of  individuals  in  a 
group  of  sus  cep  titles  who  become  infected  at  time  t  +  At  is  equal  to 
the  number  of  suscep titles  in  the  group  at  time  t  multiplied  by  a  con¬ 
stant,  \Q,  increased  by  a  spatial  weighting  of  the  infected  neighbors . 
Neighboring  infectives  will ,  therefore,  always  have  a  detrimental  effect. 
They  tend  to  increase  the  rate  of  infection  and  decrease  the  rate  of 
cure. 

Cortoining  these  assumption  yields  the  following  equation  for  the 
density  of  infectives  at  point  r  at  time  t  +  At 

f (r,t  +  At)  =  f  (r,t)  [1  -  At(y0  -  /  u( X  -  s)f  (s,t)ds] 

+  (n(r)  -  f(r,t))At[XQ  +  /  X( r  -  s)f(s,t)ds] 

where  n(r)  is  the  density  of  individuals  of  both  types  at  r,  u(r  -  s) 

is  a  scalar  function  with  a  vector  argument  that  gives  the  effect  of 

infectives  at  s  on  the  cure  rate  of  infectives  at  r ,  and  Mr  -  s) 

gives  the  effect  of  infectives  at  s  on  the  infection  rate  of  susoep- 

tibles  at  r.  Define  p(r ,t)  =  P[an  individual  at  r  is  infected  at 

f  (r,t) 

time  t],  then  p(r,t)  =  n |r^--  '  and 

p(r,t  +  At)  =  p(r,t)  [1  -  At(u0  -  fu(r-  s)n(s)p(s,t)ds) 

+  (1  -  p(r,t))At[X0  +  /X(r  -  s)n(s)p(s,t)ds) 
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9p(r,t)  , 

— —  =  +  +  "  §)  "  M(r  -  s))n{s)p(s,t)ds] 

+  Aq  +  /A  ( r  -  s)n(s)p(s,t)ds 

In  our  exarple  system  each  Individual  (node)  occupies  a  point  on  the 
integer  lattice ,  and  since  each  node  is  connected  only  to  its  four  near¬ 
est  neighbors,  we  have 

A  (r  -  s)  * 

and  y  (r  -  s)  * 

The  result  is  a  system  of  differential  equations  which  relate  the 
probability  that  any  node  is  bad  at  time  t  to  the  probability  that 
other  nodes  are  bad  at  time  t.  The  equation  for  a  non-border  node  at 
r  is 

3P(r,t) 

- —  *  -p(r,t)  IyQ  +  Ag  +  (A  -  y)  (pfs^t)  +  p(s2»t) 

+  p(s3,t)  +  P(§4»t))  J  +  Aq  +  A  (p(s^,t)  +  ptsgrt) 

+  P(s3,t)  +  pfs^t)) 

where  s^,  s^,  and  s^  are  the  four  nearest  neighbors  to  the  node  at  r. 
Values  for  the  parameters  Ag,  A,  yQ,  and  y  are  obtained  in  the 
following  way: 


A  for  |r  -  s|  <  1 
0  otherwise 

y  for  |r  -  sj  <  1 
0  otherwise 
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*  a  -  y(0)  +|  y(0)  «  XQ  +  kA 

(k)  (0)  k  (0)  - 

M  =  M  -  j  y  -  yQ  -  ky 


,  A- 


p(0) 

“5T 


Substituting  these  values  into  the  differential  equation  yields 

=  -op(r,t)  +  a  -  y*°*  +  Hg — (pfs^t)  +  ptsj/t)  +  pfs^t)  +  pfs^t)) 

thick  was  obtained  in  Bg.  (16)  from  a  strict  probabilistic  model. 

Adjacent  nodes  have  nearly  equal  probabilities  of  being  blocked. 
Consider  the  case  when  all  of  these  probabilities  are  exactly  equal  (as 
an  approximation) .  Then  from  Eq.  (16) 

*  -crp(t)  +  o  -  y<°>  +  i  y*0)p(t) 

=  -(o  -  iy(0))p(t)  +  o  -  y(0) 

D 

which  has  the  solution 


p(t) 


p(0)  - 


a  -  y 


(0) 


O-iu<0) 

5 


y°> 


it 


o  -  y 

77T-v 


(0) 


W 


(17) 


which  will  be  assured  to  hold  for  o  >  y^ . 

Now  consider  the  alternating  Poisson  renewed  process  shown  in  Fig. 
4.  There  are  two  states,  called  blocked  (B)  and  free  (F).  If  the 
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system  is  in  state  B  at  time  t,  it  goes  to  state  F  in  the  next  instant 
At  with  probability  (y(°)/5)At.  In  similar  fashion,  the  probability 
that  it  leaves  state  v  and  re-enters  state  B  is  (o  -  y^)At.  There¬ 
fore,  the  probability  that  it  is  in  the  blocked  state  at  time  t  +  At  is 

PB(t  +  At)  =  Pg<t)  (1  -  At)  +  (1  -  p^t))  (o  -  y(0))At 


-  Pb(0) 


-  o-M(0)l 


(0) 

T  . — 


o-iylO) 

5 


23 


This  is  the  same  as  Eq.  (17)  which  was  obtained  for  the  probability  that 
a  node  is  blocked  at  time  t!  In  a  large  homogeneous  network,  the  frac¬ 
tion  of  blocked  nodes  may  be  closely  approximated  by  the  probability  that 
any  one  rf  them  is  blocked.  Therefore,  the  fraction  of  blocked  nodes  at 
t-jire  t  in  a  large  uniformly  connected  (i.e.,  two-dimensional  lattice) 
network  is  approximately  equal  to  the  wrobability  that  the  two-state 
Markov  process  shown  in  Fig.  4  is  in  the  blocked  state  at  time  t.  Thus 
we  may  take  this  two-state  Markov  process  as  a  model  for  the  network. 

So  far  we  have  presented  only  aggregate  results.  To  obtain  the 
probability  that  any  given  node  in  the  network  is  blocked  at  time  t  we 
must  consider  a  system  of  equations  of  the  form  (see  Eq.  (16)) 

*l(t)  ,  -op.  (t)  +  o  -  Ii(0>  +  (t)  +  put)  +  pt(t)  +  Pm(t» 

dt  1  d  J 

for  each  node  i  in  the  network  with  neighbors  j,  k,  and  m.  These 
equations  are  obviously  of  the  form 

P(t)  -  AP(t)  +  C  (19) 

If  there  are  M  nodes  in  the  net,  then  P(t)  is  the  M  x  1  matrix 
whose  i**1  oonpenent  is  the  probability  that  node  i  is  blocked  at  time 
t.  A  is  an  M  x  M  constant  matrix  and  C  is  an  M  x  1  constant 
matrix.  The  solution  is  well  known: 

P(t)  =  eAtP(0)  +  A”1  (eAt  -  I)C  (20) 

For  a  small  net  this  solution  poses  no  difficulty,  but  for  a  large 
one  the  required  matrix  oonputations  rapidly  get  out  of  hand.  There  are 
seme  special  cases  which  are  solvable,  however,  and  we  obtain  the  solu¬ 
tion  for  one  of  these  belaf. 
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This  observation  holds  for  a  square  grid  with  any  nunber  of  nodes  n  on 
a  side.  (See  Appendix  A  which  gives  the  conplete  solution  for  P(t) 
with  arbitrary  n  for  this  network  configuration  and  two  others.) 

The  network  model  predicts  that  the  equilibrium  fraction  of  blocked 
nodes  is  zero  for  the  case  a  -  y  ^ .  For  an  infinite  value  of  N  (the 
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storage  size  in  the  D**)  this  result  would  be  obtained.  However,  for 

finite  N  the  equilibrium  fraction  of  blocked  nodes  is  non-zero.  To 

obtain  an  expression  for  this  equilibrium  value  we  must  look  at  the 

different  topologies  of  connected  blocked  nodes,  which  we  call  cluttps. 

As  a  by-product  of  this  analysis  we  will  also  get  the  clunp  size  distri- 

(0) 

buticn  for  the  case  a  -  p 

C.  Clutping  Analysis 

1.  Definition  of  a  Clurp 

ftor  a  lattice  network  in  which  each  node  has  exactly  four 
neighbors  (adjacent  nodes)  we  wish  to  define  a  clurp  of  blocked  nodes. 
TWo  blocked  nodes  are  in  the  same  clurp  if  they  are  adjacent  or  are 
linked  to  each  other  throu^i  a  series  of  adjacent  blocked  nodes.  A 
blocked  node  that  is  surrounded  by  four  free  nodes  is  a  clurp  of  size 
one. 

2.  Markov  Chain  Model  for  Clurp  Growth 

Suppose  o  -  and  that  the  expected  fraction  of  blocked 

nodes  is  very  low,  say  less  than  .1.  Then  the  probability  of  the  inter¬ 
action  of  two  cluips  is  very  small,  being  on  the  order  of  .01,  and  we 
are  justified  in  locking  at  the  growth  of  clumps  from  single  nodes  (as 
an  approximation) .  Thus  we  will  neglect  the  possibility  that  two  clumps 
aonbine.  The  simulation  results  (described  later)  indicate  that  this 
approximation  is  good  for  a  storage  size  N  >  50.  Also,  we  neglect  the 
effect  of  the  HOSTS  since,  by  assumption,  they  cannot  became  blocked. 

Consider  one  free  node  in  the  midst  of  many  free  nodes.  It  be¬ 
comes  blocked  in  a  Poisson  fashion  at  a  rate  X  ^ .  Then  we  have  the 
following  situation: 
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X  =  blocked  node 


0  =  free  node 


This  clutp  of  one  can  become  a  clurtp  of  two  at  a  rate  4A  ^ ,  again  in 
a  Poisson  fashion,  or  die  out  at  a  rate  y  ^ .  Sippose  it  becomes  a 
clutp  of  two,  then  we  have  the  following: 


0  0 

0  X  X  0 
O  0 

This  clunp  of  two  can  become  a  cluip  of  three  at  a  rate  6A  ^ ,  or  be¬ 
come  a  clutp  of  one  at  a  rate  2y  ^ .  Suppose  it  goes  to  a  clutp  of 
three,  of  which  there  are  two  forms: 


I)  0  0  II)  O  ,0  O 

OXXO  OXXXO 

0  X  0  0  0  0 

0 

Form  I  has  a  growth  rate  of  6A^  +  A^  and  a  death  rate  of 
2y^  +  y^ ,  while  form  II  has  a  growth  rate  of  8A^  and  a  death 
rate  2y  ^  +  y  ^ .  The  death  rates  are  obviously  equal  for  the  two 
different  forms,  but,  surprisingly,  the  growth  rates  are  also.  Recalling 
that 

x<w  +£„(0) 


and 


u(k)  =y(0)  _k  y(0) 
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we  have  for  form  I  (using  \j  to  indicate  growth  rate  for  form  I) : 
Xj  -  6X^  +  X*2)  *  6(o  -  y*0)  +  ^  y*0*)  +  o  -  y*0)  +  |  y*0) 

.  7(0  .„(«>,  +  3  <°> 


XII  =  8X 


8(o  -  y(0))  +  §y(0) 


Then,  for  the  case  o  ■  y  ^  ,  we  have  Xj  ■  X  =  8X  ^ . 


nodes: 


There  are  five  different  topologies  for  a  clunp  of  four  blocked 


0  0 


0  X  X  0 


0  X  X  0 


0  0 


0  0  0  0 


0  X  X  X  X  0 


0  0  0  0 


\u  -  10X ' 


yIX  =  2y(2)  +  2y(1) 


0  0 


0  X  X  0 


0  X  X  0 


0  0 


Xttt  ■  6X (1)  +  2X  (2)  *  10X (1) 
yin  -  2y(2)  +  2y(1) 
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IV) 


0 


0  0X0 

0  X  X  X  0 

0  0  0 


IV 

PIV 


8A(1)  +  X(2)  -  10A(1) 


=  2y(2)  +  2p 


(1) 


V)  0  0  0 

0  X  X  X  0 

0X0 
0 

The  growth  and  death  rates  are,  except  for  the  square,  form  I,  the  same 
for  the  different  forms.  So,  to  determine  the  growth  and  death  rates 
for  a  clump  of  four  blocked  nodes  it  is  approximately  sufficient  to  look 
at  the  straight  line  form,  form  II.  For  larger  size  clunps,  we  consider 
only  this  straight  line  form  for  determining  the  growth  and  death  rates. 
Such  a  simplification  is,  of  course,  necessary  since  the  number  of  dis¬ 
tinct  topologies  prohibits  exhaustive  treatment.  Simulation  results 
support  this  approximation  and  show  that  elongated  clunps  are  more  likely 
to  occur  in  systems  of  this  kind  than  are  square  or  circular-shaped 
clumps  with  their  minimum  circumference  to  area  ratio.  For  a  clunp  of 
length  n  we  therefore  have  the  following : 


Av  -  6A(1)  +  2A 


VH,=  3m(1)+M 


(2)  _ 

(3) 


=  10X 


(1) 


-  2m(2)  +  M (1) 


0  0  0  0 

0  X  X  ...  X  X  0 

0  0  0  0 


Ah  »  2(n  +  1)A(1) 

Mn  =  (n  -  2)m(2)  +  2m(1) 


(25) 
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ThtB  our  approximation  leads  us  to  the  fo  ll»dng  birth-death 


process  for  duip  size: 


4X*1’  6X'1’ 


X0X(1)  X.  =  2(n+l)X(1) 

.  n 


l/J  VJJ  I**/  ^  (2)  Q) 

--TTT^Sir- in+-7fi<2r - =  + 

2y(1)  2y(1,-»y(2)  2(y  1  +y'2' ) 


Vfe  will  simplify  yn  somewhat. 


Pjj  *  (n  **  2)m(2)  ♦  »,« 


-„U<2>  -  2*|  U t0)  +24u(0) 
-nw<2U„-'3> 

*  <n  +  1)m(2) 


Therefore,  to  simplify  the  solution  we  use  the  approximations 


\  ^  *  2(n  +  1)X' 

Un  -  (n‘+  l)y(2 


n  >  1 
n  >  2 


i 

Define  p  to  be  the  equilibrium  probability  of  n  blocked  nodes 
n 

the  clurp  and  by  elementary  queueing  theory  17] 


ni° 

D  _  n  X(0)  "Z1  2 (i  +  1)X(1)  1  , 

!.*  tt*a>  nil 


where  r  * 


(n  +  l)y 


T55" 


v 


-*o^(rV^7Mrhr))(28) 


p  is  obtained  in  the  usual  way: 
0  ' 


\  i 


V 

(-  r  -  log(l  -  r))j  , 


D  .  f,  +  »««  2 

po-  1  +  7oT7 


(25) 


and 


x<°> 


pn  “  po  ~w  rrr t"'1  nil 


where  r  -  ^ 


(1) 

^2T 


(30) 


Ebr  the  case  o  =  we  thus  have  the  equilibrium  fraction  of  blocked 

nodes  (Eq.  (28) )  and  the  distribution  of  cluip  size  (Eq.  (30) ) . 

\ 

I  .  ■  \ 
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In  Appendix  B  we  apply  the  chirping  analysis  to  the  8-neighbor 
race  (shown  in  Fig.  5)  to  obtain  the  equilibria  fraction  of  blocked 
nodes  for  this  network  configuration  when  o  =  .  The  results  ob¬ 
tained  from  the  chirping  analysis  are  good  for  the  case  o  =  in 

both  4-  and  8 -neighbor  configurations .  But  the  results  are  very  poor 
for  a  >  u(0),  and  this  is  probably  attributed  to  the  interaction  of 
clunps.  We  treat  this  case  next. 


(0) 

3.  Average  Cliitp  Size  for  a  >  u 

Although  we  have  not  arrived  at  a  method  for  determining  clinp 
size  distribution  in  the  more  heavily  blocked  cases  (i.e.  a  >  \i^) ,  we 
have  a  nethod  which  gives  a  crude  estimate  of  the  average  chirp  size  for 
these  cases.  It  is  based  on  the  idea  that  any  node  is  potentially  the 
"origin"  of  a  clwp.  Hie  network  model,  Eq.  (18) ,  gives  the  equilibria 
probability  p  of  a  blocked  node  for  oi  while  the  clmping 
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analysis,  Eg.  (30) ,  gives  the  distribution  of  cluip  size  front  an  isolated 
node  (only  for  o  =  y^) .  To  treat  the  case  a  >  \i^  we  must  oonbine 
these  ideas. 

tfe  assure  that  cluips  occur  as  overlaps  of  clurps  front  origin 
nodes  which  are  distributed  miformly  across  the  lattice  with  probability 
p.  It  would  seem  that  a  problem  of  conservation  of  blocked  nodes  might 
exist,  but  for  estimating  the  average  cluip  size  this  method  gives  good 
approximate  results. 

Let  ib  take  the  left-hand  extremity  of  a  cluip  as  its  "origin." 
Vfe  will  use  a  simplified  clurping  analysis  that  assumes  cluips  are  always 
linear  with  growth  or  death  occurring  at  the  ends,  This  is  generally  a 
poor  approximation,  but  it  has  the  advantage  that  the  length  of  a  cluip 
is  then  geometrically  distributed  and  analytic  results  are  possible.  In 
particular,  we  will  find  the  probability  that  a  point  is  neither  an 
origin  nor  is  "covered"  by  a  cluip  and  call  this  P (empty  "system"]. 

The  relationship  of  this  system  to  an  infinite  server  queueing 
system  (NVV“)  will  be  shown.  Using  arguments  similar  to  those  used  in 
(7]  to  get  the  average  length  of  a  busy  period  in  a  single  server  system, 
we  will  get  the  average  length  of  a  one-dimensicnal  cluip  for  the  case 
a  >  y(0).  Finally,  we  will  employ  three  different  topologies  for  the 
average  two-dimensional  cluip  and/or  different  interpretations  for  the 
average  one-dimensional  cluip  length  to  get  estimates  of  the  average 
cluip  size  for  the  two-dimensional  case  with  a  > 

Let  us  suppose  that  nodes  are  marked  (blocked)  with  probability 
p,  the  equilibrium  fraction  of  blocked  nodes  obtained  from  Eq.  (18)  by 
letting  t  ■*>  «.  Associated  with  each  marked  point  is  a  length,  geometri¬ 
cally  distributed,  which  extends  out  to  the  right  as  in  Figure  6. 
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Flaunt.  Ctumptnf  MotM  for  a  >nM 

p  (p  =  k)  **  (1  -  0)0*  ^  k  ■  1,2,... 

i 

where  k  =  1  corresponds  to  a  clurp  of  size  one 

00  i~l  k 

P(t  £  k)  -  1  -  P(A  >  k)  -  1  -  £  (1  ~  “  1  ”  a 

j=k+l 

If  a  point  is  not  covered  by  a  line  or  a  marie,  then  we  say  that  the  sys¬ 
tem  (i.e.  point)  is  "enpty." 

Plenpty  system]  i  pQ  -  (1  -  p)  f]  (d  -  P)  +  Pd  “  °k>) 

k-1 

00 

-  ]"[  d  -  pak)  <31> 

k=0 
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Using 


log(l  -  x) 


12  13  1  4 

~X  ‘  2  X  4  X  ‘  4  X 


we  have 


,  “  #  k  l,  k  2 

log  p  =  (-po  -  7(po  ) 

0  k*0  1 


1/ 

j(po  ) 


-JL-  4jL 

l-o  ?  1  -  a2 


1  P 
T  '  3 

J  1  - 


,  -  esq)  - j 

0  j=l  j(l-°j> 


(32) 


tfe  make  the  restriction  p  <  1/2  since,  analogous  to  the  conjectured 
exact  result  for  the  critical  probability  in  percolation  theory  (see 
Chapter  2,  "Related  Work") ,  the  probability  of  an  infinite  clutp  nay  be 
non-zero  for  the  case  p  2  1/2.  Approximating  pQ  by  the  first  term 
only,  we  have 


Po s 


(33) 


Let  us  aotpare  this  model  to  an  infinite  server  queueing  system 
dVH/°°) .  The  average  in  ter  arrival  time  of  customers  to  such  a  system  is 
1/A  seconds  and  a  customer  departs  after  receiving  an  average  of  1/u 
seconds  of  service.  For  this  case  we  have 

P[enpty  system,  i.e. ,  no  customers]  =  e  -  jj  =  PQ 


In  our  nodal  blocking  system  the  "average  interarrival  distance"  between 
marked  points  (in  nodes)  is 
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i  =  p  (1  +  2  (1  -  p)  +  3(1  -  p)2  +  ...  ) 


=  p  f)k{l  -  p) k_1 
k=l 


Let  X  =  1  -  p 


then 


(1  -  X)‘ 


1 

P 


The  "average"  service  distance"  (in  nodes)  is 


i-fkd- 0,^-1 

v  k*i 


1-0 

- 2 

(i  -  or 


and  P0  =  et‘e"  T^o  “P0 

Thus  the  system  we  are  considering  corresponds  approximately  to  an  in¬ 
finite  server  queueing  system. 

Our  systan  is  "enpty"  with  probability  pQ  and  "busy"  with  probabil¬ 
ity  1  -  pQ.  In  any  line  of  N  nodes  or  points  (N  »  1)  NpQ  will, 
on  the  average,  be  enpty  and  N (1  -  Pg)  will  be  busy  (see  Fig.  6) .  The 
average  length  of  an  enpty  string  is  the  average  interarrival  distance 
for  our  system  *  1/p  nodes.  Therefore,  the  NpQ  empty  nodes  will,  on 
the  average  conprise 

Np0 

I7p  =  ^0  p 


distinct  enpty  sets  or  strings.  Therefore,  the  average  length  of  a  busy 
string  is 


N(1  -  Pg) _  1  -  Pq 

Np0  p  Pq  P 


nodes 


(34) 
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which  is  analogous  to  the  average  length  of  a  busy  period  in  an  M/tVl 
queueing  system  [7] . 

We  irast  still  determine  o ,  the  parameter  in  the  geometric  length 
distribution  Wc  do  this  by  oiisidering  a  line  of  n  block?  d  nodes  and 
assuming  that  growth  or  death  can  only  occur  at  the  aids  of  the  string. 
Then  we  haw  the  following: 

OXX...XXO  X  =  blocked  node 

«  n  ■  — >  0  =  free  node 

This  iirplips  the  following  birth-death  process  for  chain  length: 


Define  =  P  [chain  is  of  length  n] ,  then 


\{1) 

Clearly,  we  should  take  a  =  — rjr  in  our  clurtp  model 

There  are  at  least  three  possible  approaches  to  the  determination 
of  the  average  clurtp  size  C: 

1  “  Pq 

I)  Assume  all  of  the  clunps  are  circles  of  radius  R,  and  i  *  p^  p 
is  the  average  length  of  the  intersection  of  a  random  line  with  a  circle 
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of  radivs  R.  Kendall  and  Moran  [221  give  the  average  length  of  the 
intersection  to  be  1/2 ttR.  Then  we  have 


and 


(35) 


where  pQ»e  -  p/(l  -  a)  and  p  is  the  equilibrium  fraction  of  blocked 
notes  obtained  from  the  network  model. 

II)  Assume  T  is  the  diameter  of  an  average  clunp  (assumed  circular) , 
then 


(36) 


III)  Assure  1  is  the  length  of  the  side  of  an  average  clunp  (assured 
square) ,  then 

C  “  I2  (37) 


Ibr  the  case  which  prompted  this  analysis  all  three  of  these  methods 
give  an  average  clunp  size  within  .9  of  the  value  observed  in  simulations 
(approximately  3.48) .  Method  I  overestimates  the  observed  value  by  .76, 
method  II  underestimates  it  by  .86,  and  method  III  underestimates  it  by 
.16. 


4.  Ftexinum  Clunp  Size 

A  model  which  predicts  the  size  of  the  largest  clunp  surpris¬ 
ingly  well  was  suggested  to  the  author  by  Mr.  Tom  Leavitt  of  the  UCLA 
Conputer  Science  Department.  In  previous  sections  we  assured  that 
"stringy"  clumps  are  more  common  than  round  or  square  ones  because 
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growth  in  a  probabilistic  system  occurs  by  shooting  out  projections  in 
random  directions.  These  random  projections  actually  "weaken"  the  clunp 
by  exposing  it  to  more  free  nodes.  We  expect  the  largest  clumps  to  show 
a  tendency  to  minimize  their  circumference  with  respect  to  their  area. 
Therefore,  in  modelling  the  largest  clumps  we  will  use  rectangular 
clunp  topologies.  We  assure  that  a  clunp  will  increase  in  size  until 
the  number  of  free  nodes  on  the  border  that  are  becoming  blocked  is 
equal  to  the  number  of  blocked  nodes  chi  the  border  that  are  becoming 
free.  This  equilibrium  point  corresponds  to  the  largest  clunp.  In 
order  to  perform  the  analysis  we  must  make  the  following  assumptions: 

1)  all  clumps  are  rectangular 

2)  blocked  nodes  not  on  the  border  will  remain  blocked 

3)  every  blocked  node  on  the  border  has  exactly  three  blocked 
neighbors 

4)  every  free  node  on  the  border  has  exactly  one  blocked  neighbor 
An  exanple  of  such  a  clunp  is  that  shown  in  Fig.  7. 


X- BLOCKED  NODE 
0  -  FREE  NODE 
III  ■  BLOCKED  BORDER 
^ ■ FREE  BORDER 


Fifum  7.  Maximum  Clump  Six*  Modal 
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He  see  that  the  lumber  of  blocked  nodes  on  the  border  is 

21  +  2(w  -  2)  *  2(1  +  w  -  2) 

and  the  number  of  free  nodes  cn  the  border  is 

21  +  2(w  +2)  *  2(1  +  w  +  2) 

At  equilibrium  we  have 

2(1  +  w  -  2)y(3)At  *  2(1  +  w  +  2)A(1)At 
or 


Far  a  fixed  border  size,  die  nurber  of  nodes  in  the  clunp  is  maximized 
for  l  ■  w,  or 


Therefore,  the  expected  maximum  clunp  size  is 


There  are  two  reasons  why  this  result  estimates  the  maximun  clunp  size 
and  not  the  average  clunp  size: 

1)  Clutps  do  not  grew  by  adding  entire  borders;  they  add  projec¬ 
tions  that  weaken  the  clunp. 

2)  The  model  assumes,  incorrectly,  that  blocked  nodes  within  the 
clunp  cannot  become  free;  but  they  do  and  this  further  weakens 
the  clunp. 

A  number  of  models  and  results  have  been  presented  to  characterize 
the  behavior  of  a  network  of  two-stage  Markovian  nodes.  The  efficacy 
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of  these  methods  will  be  shown  in  the  next  section  in  which  we  discuss 
the  network  s  inula  ticn. 
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CHAPTER  4 


MAFKOV  MDDEL  NETWORK  SIMULATION 


A.  Description 

Simulation  of  a  network  of  1024  nodes  enploying  the  Jfarkovian  in¬ 
ter-event  time  assumption  has  substantiated  the  analytical  approxima¬ 
tions  described  earlier.  The  two  different  programs  which  simulated 
this  network  are  listed  in  Appendix  C.  These  programs  run  cn  the  UCLA 
XDS  Sigma-7  ccnputer. 

The  first  program  simulates  a  network  arranged  in  a  square  grid 
32  x  32  and  simultaneously  displays  the  net  activity  on  a  Digital  Equip¬ 
ment  Corporation  340  Precision  Display  CRT  (Fig.  8) .  Each  node  is  con¬ 
nected  to  its  four  nearest  neighbors  (a  lattice)  except  in  the  case  of 
the  nodes  along  the  border,  which  have  only  three  nearest  neighbors  (or 
two  nearest  neighbors  in  the  case  of  the  four  comer  nodes) .  When  a 
node  changes  state,  new  event  times  are  chosen  for  it  and  for  all  of  its 
nearest  neighbors  based  on  the  new  nuntoer  of  blocked  neighbors.  The 
nenoryless  property  of  the  exponential  distribution  simplifies  the  cal¬ 
culations. 

The  second  program  simulates  a  randomly  connected  graph  in  which 
each  node  is  given  exactly  four  neighbors.  Due  to  rranory  size  limita¬ 
tions,  this  program  does  not  have  a  graphical  display. 

B.  Comparison  of  Observations  and  Predicted  Behavior 

1.  Fraction  of  Nodes  Blocked 

Comparison  of  the  network  model  (Eq.  (18))  and  the  simulation 
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Figure  8.  Network  Simulation  CRT  Display 
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results  for  the  lattice  ond  the  randan  graph  are  shewn  in  Figs.  9,  10, 
and  11  for  three  different  sets  of  system  parameters  a  and  each 

starting  both  from  completely  blocked  and  corrpletely  free  nets.  In  Fig. 

9  the  equilibrium  fraction  of  blocked  notes  is  obtained  from  Eq.  (28)  of 
the  clunping  analysis. 

At  any  point  in  time  the  network  model  (Eq.  (18) )  predicts  sane 
value  f  as  the  expected  fraction  of  blocked  notes.  Assuming  no  corre¬ 
lation  between  nodal  states  and  a  network  having  1024  nodes,  8,  the 
standard  deviation  of  the  measurement  of  the  fraction  blocked  is  [23] 

„  _  Vfc(l  -  f) 

0  - 55 — 

At  equilibriun  we  have  in 

Figure  9:  f  *  .07  8  ■  .00796 

Figure  10:  f  ■  .25  8  =  .0135 

Figure  11:  f  =  .833  8  *  .01165 

With  a  95%  confidence  limit  of  1.968  and  a  99.7%  confidence  limit  of  38, 
we  see  that  in  Fig.  9  the  as  sorption  of  independence  is  completely  un¬ 
acceptable.  Recalling  that  the  equilibrium  value  for  this  case  was  pre¬ 
dicted  fran  the  clumping  analysis  which  shows  a  high  degree  of  correla¬ 
tion,  the  deviations  ab served  in  the  equilibrium  value  in  Fig.  9  are  not 
surprising.  The  behavior  observed  in  Figs.  10  and  11  is  generally  with¬ 
in  the  99.7%  confidence  limit.  Occasional  excursions  outside  this  range 
show  the  effect  of  clurp  formation  and  dissolution. 

Figures  12,  13,  and  14  give  simulation  results  for  the  two-dimen¬ 
sional  integer  lattice  in  which  each  node  is  assumed  to  have  eight 
neighbors .  This  was  accomplished  by  extending  the  nearest  neighbor  defi- 


48 


S3CJ0N  03)13018  30  NOliOVUd 


49 


Figure  12.  Fraction  of  Block  ad  Nodes  IV 
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nition  to  include  nodes  which  are  diagonally  adjacent.  The  random  graph 
program,  because  of  oonputer  matory  limitations,  could  not  be  modified 
to  include  the  8-neighbor  case.  In  these  figures  comparison  is  made  to 
the  predicted  behavior  obtained  from  the  network  model  assuming  every  IMP 
has  exactly  nine  output  lines,  one  of  which  goes  to  the  HOST.  The  equil¬ 
ibrium  fraction  of  blocked  nodes  in  Fig.  12  is  obtained  frcm  the  clump¬ 
ing  analysis  given  in  Appendix  B. 

Figures  15,  16,  and  17  ocrpare  simulation  results  an  the  lattice  of 
degree  four,  when  a  free  node  with  k  blocked  neighbors  is  considered 
k-fourths  blocked,  to  the  predicted  behavior  based  on  a  non-linear  "par¬ 
tial  blocking"  model,  fliis  model  makes  two  assumptions : 

1.  The  disturbance  (i.e.,  blocking  propagation)  spreads  out  in  a 
wave-like  manner  from  blocked  nodes  and  can  be  characterized  as 
a  Poisson  growth  process  of  the  type  studied  by  Morgan  and 
Welsh  [15] .  In  particular,  we  assure  that  the  blocking  starts 
with  a  single  blocked  node  in  the  center  of  the  network  and 
that  blocking  is  limited  to  what  we  call  the  "disturbed  area" — 
those  nodes  which  are  within  a  distance  r(t)  of  the  center 
node. 

2.  If  the  number  of  blocked  nodes  within  the  disturbed  area  (com¬ 
prising  a  total  of  N(t)  nodes  is  n(t) ,  then  the  nuriber  of 
blocked  neighbors  k(t)  seen  by  an  average  node  within  the 
disturbed  area  is 


k(t)  =  4  SUf  +  4 


(1 


n(th 

NltP 


n(t) 

nTET 


where 
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Figure  IS.  Fraction  of  Blocked  Nodes  VII 
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Figure  16.  Fraction  of  Blocked  Nodes  VIII 
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n(t  +  At)  =  n(t)  -  n(t)y(k(t))At  +  (N(t)  -  n(t) )  X(k(t) }  At 
N(t)  is  found  by  assuming  that  a  free  node  on  the  edge  of  the  "disturb¬ 
ance  wave"  sees  on  the  average  1-1/2  blocked  neighbors  as  pictured  below: 
X  0 


These  equations  must  be  integrated  step  by  step.  The  results  are  gener¬ 
ally  poor  except  in  the  case  a  =  .02,  which  is  relatively  insensitive 
to  changes  from  the  basic  4-neighbor  network  model. 

2.  Distribution  of  Clump  Size  for  a  =  y^ 

Figure  18  compares  the  equilibriun  distribution  of  clump  size 
observed  in  the  4-neighbor  lattice  simulation  to  the  prediction  based  on 
Eq.  (30) .  Figure  19  gives  the  expected  clump  size  distribution  in  a 
lattice  when  the  blocked  nodes  are  placed  randomly  on  the  lattice  with 
an  average  fraction  blocked  of  .07  as  given  by  Roach  [16] .  We  compare 
this  to  the  clump  size  distribution  observed  in  the  random  graph  for  the 
pgcp  o  =  ,  n  =  50  by  formally  assigning  that  the  nodes  are  in  a 

lattice.  The  result  of  this  assumption  is  a  mapping  that  randomly  dis¬ 
perses  the  clumps.  The  agreement  is  excellent,  and  by  comparing  Figs. 

18  and  19  we  see  that  the  clumping  in  the  Markov  network  is  not  at  all 


random  (i.e. ,  uncorr elated) . 
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Figure  19. 


3.  A\ 


Clurp  Size  for  o>y 


i  i  ... 

Pear  the  case  0  =  .01067,  y  U  =  .01  an  average  clurp  size  of 

\ 

3.48  was  observed  in  the  simulation  after  equilibrium  was  attained.  Thdj 
three  different  methods  for  predicting  this  value  give  estimates  of  4.24, 
2.62,  and  3.32,  respectively,  straightforward  application  of  the  clump¬ 
ing  analysis  (Eq.  (28)) ,  which  is  valid  for  o  =  y^ ,  yields  a  value 
greater  than  6.  Hence  these  new  methods  offer  seme  improvement. 

4.  Maximum  ClUmp  Size 

- I  -  1  \ 

Figures  20  and  21  show  the  distribution  of  the  maximum  clurp 

size  observed  in  the  similation  for  two  different  sets  of  parameters 

after  equilibrium  is  reached.  Figure  21  shews  the  effect  of  "harmonics" 

of  the  expected  maximum  clunp  size  as  large  clurps  combined  for  short 

,  •  .  l 

times.  The  results  are  renerkably  good,  especially  considering  the  dis- 

i  I 

persion  in  the  distribution  in  Fig.  21. 

\ 

\ 

C.  "Hot  Spots"  -  Analysis  and  Results  \ 

In  this  section  vie  analyze  tile  effect  of  placing  a  small  number  of 
high  rate, of  blocking  (i.e.  o»y(0))  nodes  into  networks  of  predomi¬ 


nantly  low  rate  of  blocking  nodes  (o<y';) .  We  call  these  high  rate  of 
blocking  nodes  "hot  spots" .  The  simulation  of  a  single  hot  s^xjt  (with 
o  =  2y*°* )  placed  centrally  in  a!  32  x  32  network  of  nodes  with 
a  s  /2  revealed  that  such  low  rate  of  blocking  nodes  effectively 
prevent  blocking  propagation.  The  high  rate  of  blocking  node  was  the 
only  node  in  the  network  that  was  ever  observed  to, block.  Hence  in  the 


analysis  to  follow,  the  lew  rate  of  blocking  nodes  will  be  assumed  to 

have  a  =  y^ ,  N  =  50,  and  we :  will  approximate  the  hot  spots  as  being 

\ 

permanently  blocked. 
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Line  of  N  Permanently  Blocked  Nodes 


X  X  ...  X  X 

- - N 

Suppose  a  line  of  permanently  blocked  nodes  is  put  into  an  environ- 
nent  of  nodes  with  o  =  and  N  =  50.  For  a  network  consisting  en¬ 
tirely  of  this  latter  kind  of  node,  we  know  that  the  expected  maximum 
clunp  size  is  9  nodes.  This  leads  ui  to  expect  a  triple  row  of  N 
blocked  nodes,  including  those  permanently  blocked.  Therefore  in  a  net 
of  1024  nodes,  recalling  that  .07  is  the  expected  fraction  of  blocked 
nodes  in  the  absence  of  permanently  blocked  nodes,  we  should  have,  with 
our  blocked  line, 

E  [fraction  blocked]  =  (3N  +  .07(1024  -  3N))/1024 

«  .07  for  N  small  (39) 

For  N  =  32,  i.e.,  the  line  of  permanently  blocked  nodes  spanning 
the  network,  we  should  get 

E  [fraction  blocked]  -  (32  *  3  +  .07(1024  -  32  *  3))/1024 

=  .157  (40) 

For  isolated  permanently  blocked  nodes  we  must  again  consider  the 
growth  topologies  and  the  Markov  chain  structures. 

Let  indicate  a  permanently  blocked  node 

X  indicate  a  temporarily  blocked  node 

®  X1  =  4X(1) 
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We  see  that  from  a  single  permanently  blocked  node  growth  occurs  at  a 
rate  of  4X^ .  Let  us  look  at  a  clump  of  two  and  form  the  correspond¬ 
ing  Markov  chain: 


X  ® 


x2  =  6A(1) 


.  ,.d) 


y?  =  y 


=  ,,.(1) 


y,  =  3y 


The  death  rate  out  of  state  3  (i.e.,  a  clvirp  of  3)  assumes  that 
either  of  the  following  topologies 


X  ®  X 


or 


X 

X  ® 


is  much  more  likely  than 


XX® 

Already  we  have  been  forced  to  make  approximations.  The  topological 
problems  which  we  face  in  this  analysis  are  even  more  difficult  than 
those  faced,  before  in  analyzing  the  system  to  obtain  the  average  nurrber 
blocked  for  the  case  a  =  y^ .  At  that  time  we  found  it  useful  to  make 
the  approximation 

Xn  =  2(n  +  1)  (1)  n  >  1 

yn  =  (n  +  1)  {2).  n  >  2 

In  the  system  with  permanently  blocked  nodes  the  growth  rate  at 
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different  clutp  sizes  should  be  the  same  as  those  above.  Hcwever,  the 

pemanently  blocked  node  cannot,  by  definition,  become  free.  Assume  that 

it  is  well  within  the  clutp  at  larger  clutp  sizes .  Then  we  should  use 

the  y  qiven  above  diminished  by  where  km  is  the  highest 

n 

muter  appearing  and  the  expression  for  v^.  Hence  we  will  assume  the 
following  growth  and  death  rates: 


yn  =  2{n  +  1)X' 


n  >  1 


yn-nw 


n  >  2 


n  >  1 


"fr1  2(i  +  1)A(1) 


n  /u  + 


PjTn-1  n  ^  1 


where  r  = 


— 1ST 

y' 


00  p 

r'  n  _  *1 


£pn  -  1  -  Pi  Er n  =  rT7 
n=l  n  1  n=0  1 


Therefore 


Pn  =  1  -  r 


=  (1  -  r)r  n  >  1 


00  00  1 

in  system]  =  £np  =  d  “  r)  111:11 

n=l  n=l 

1  2X  ^ 

=  —  v^rer  =  -75y 
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For 


o  —  .01  —  p(0),  N  =  50 

Ad)  .  0  .  p(l)  _  0  .  ,(0)  +  M 


(1) 


=  .002 


H<2)  »  p<°>  -  |y(0)  =  .006 


Therefore 


.  .  2X(1)  2 

r  T55"5 


Therefore 


E[#  in  system]  =  j-i-j  =  3 

3 


Thus  an  isolated  permanently  blocked  node  should,  on  the  average,  cause 
a  clunp  of  size  3  to  be  produced,  i.e. ,  itself  plus  two  temporarily 
blocked  nodes.  If  there  are  N  isolated  permanently  blocked  nodes  and 
N  is  less  than,  say,  100  we  should  have 


E [fraction  blocked]  =  (3N  +  .07(1024  -  3N))/1024  (41) 

If  N  is  large,  i.e.,  greated  than  100,  we  must  iterate  to  a  solu¬ 
tion  as  in  the  following  exanple.  Consider  a  lattice  of  256  permanently 
blocked  nodes  superimposed  on  the  1024  node  network: 


2 

0 
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The  nunbei  =ide  a  node  indicate  hew  many  permanently  blocked  nodes 
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that  rode  has  as  neighbors.  It  is  easy  to  see  that  one-third  of  the 
non-permanently  blocked  nodes  sire  of  the  0  type,  and  the  other  two- 
thirds  are  of  the  2  type. 

Front  the  amount  of  time  spent  in  the  blocked  state  and  the  free 
state,  we  know  that  a  rode  with  X  blocked  neighbors  is  blocked  with 
probability 

i(jc)  X(x)  _  X(x) 

f*  '  £(x)  +  i(x)  =  *lX>  +  =  ~ 


Therefore,  we  have  as  a  first  step  in  the  solution 


E[#  blocked]  =  256  +  512  f2  +  256  fQ 


with 


x(2) 

f2  =  ~ 


o  -y(0)  +  5m(_^_  =  2  ^  f  al 50  =  _Q2 

o  3  0  O 


Let  =  average  #  of  blocked  neighbors  for  a  type  2  node 
kQ  =  average  #  of  blocked  neighbors  for  a  type  0  node 
then  our  iteration  proceeds  as  follows: 

[k2  =  2  +  2*f0  =2  +  2  ( .02)  -  2.04  «  2 

: =  0  +  4  *  f 2  -  4( .4)  =  1.6 


X(k2}  *2 

o  5 


.32 

o  5 


k2  =  2+  2*f0=2  +  2(.32)  =  2.64 

kft  =  0  +  4  *  f,  =  4(  .4)  =  1.6 
VO  ^ 
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In  the  limit  the  E  [fraction  blocked]  =  .66176.  (42) 

The  last  case  which  we  will  consider  is  that  of  an  R  X  R  clunp  of 
permanently  blocked  nodes,  with  R  >  2.  Modelling  the  border  of  this 
clump  as  a  line  of  permanently  blocked  nodes  formed  into  a  square,  we 
should  expect  the  cluip  to  increase  to  (R  +  1)  X  (R  +  1) . 
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Therefore, 


Elfraction  blocked]  =  ((R  +  1)^  +  .07(1024  -  (R  +  1)  ))/1024  (43) 


Table  1  lists  the  results  observed  in  the  simulation  of  hot  spots  on 


the  32  x  32  grid  for  the  following  cases: 

1.  TWo  hot  spots  side  by  side 

2.  TWo  hot  spots  separated  by  one  lew  rate  of  blocking  node 

3.  Three  hot  spots  in  a  connected  straight  line 

4 .  32  hot  spots  in  a  line  (one  whole  row  of  the  network) 

5.  A  lattice  of  64  hot  spots  spread  evenly  over  the  32  x  32  grid 

6.  A  lattice  of  256  hot  spots  spread  evenly  over  the  grid 


7.  Pour  hot  spots  in  a  2  x  2  clurp 

8.  Nine  hot  spots  in  a  3  x  3  clump 


9.  25  hot  spots  i 


Time 


Case 

%  Blocked 
Hiqh 

1 

10.0 

1327 

2 

8.4 

953 

3 

8.6 

1901 

4 

16.8 

2412 

5 

25.4 

838 

6 

64.3 

632 

7 

9.6 

1897 

8 

10.7* 

2237 

9 

10.7 

1731 

i  a  5  x  5  clurp 

Total 

%  Blocked  Observation 


Average 

Time 

7.5 

1636 

7 

1331 

7 

1985 

13.5 

3581 

24 

1265 

63 

758 

7 

2060 

8.4* 

2380 

9.2 

2098 

%  Blocked  Pertinent 
(Prediction)  Equation 


7 

39 

7 

39 

7 

39 

15.7 

40 

24.4 

41 

66.2 

42 

7 

43 

7 

43 

10.2 

43 

HOT  SPOTS  RESULTS 
TABUS  1 

♦The  high  value  and' "the  overall  greater  average  were  due  to  the  forma¬ 
tion  of  a  large  clurp  that  was  not  connected  to  the  3  x  3  cluip. 
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These  models  have  clearly  proven  their  applicability.  This  completes 
our  analysis  of  hot  spots. 

So  far  we  have  permitted  ourselves  the  strong  assurrption  of  two- 
state  Markovian  nodes.  In  the  next  section  we  treat  the  application  of 
these  results  to  a  simulated  caiputer-cximunication  network  of  64  nodes 
which  has  many  real  world  properties. 


CHAPTER  5 


SBMATICN  OF  A  NETWORK  WITH  MJSSAQ3  TRANSFER 
A.  Description 

A  program  which  simulates  a  s tore-and- forward  ccmruni cation  network 
of  64  nodes  was  run  on  the  CCLA  XDS  Sigma-7  conputer  (see  Appendix  C  for 
a  listing  of  this  program) .  In  this  network  messages  are  sent  from  ori¬ 
gin  to  destination  nodes  under  nearly  fixed  routing  strategies.  The  es¬ 
sential  characteristics  of  this  sinula tion  network  are  the  following: 

1*  Nodes  are  arranged  in  an  8  x  8  grid  and  are  numbered  consecu¬ 
tively  from  1  to  64  by  rows.  Any  node  i  is  connected  to  nodes  i  +  1, 
i  +  8  modulo  64.  The  result  is  a  "twisted  torus,"  which  shows  complete 
symmetry  for  each  node.  (A  torus  network  prevents  the  center  of  the  net 
from  becoming  a  bottleneck,  and  a  "twisted  torus"  is  conveniently  pro¬ 
grammed.) 

2.  Message  lengths  are  exponentially  distributed  with  an  average 
of  5/u^  units. 

3.  Every  node  has  storage  for  exactly  N  messages  (1  <  N  <  50) . 

4.  The  arrival  rate  of  requests  for  inputs  to  the  IMP  from  the 
HOST  is  (a  -  4n(0)/5) . 

5.  When  a  blocked  node  becomes  free,  each  of  its  neighbors  who  has 
a  message  for  it  makes  a  request  to  send  that  message  to  it  at  a  rate  of 
0  FETRSf  (or  just  a  RE) . 

6.  Routing  is  fixed.  The  routing  algorithm,  after  being  queried 
by  a  node,  relays  to  that  node  the  "best"  next  node  and  the  "second  best" 

Preceding  page  blank 
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next  node  for  that  message  based  on  its  final,  destination.  However, 
every  queue  within  a  node  far  an  output  line  from  that  node  is  limited 
in  length  to  N/4  +  1.  This  avoids  the  "deadly  embrace"  that  could  re¬ 
sult  if  two  adjacent  nodes  should  fill  up  with  messages  for  the  other 
and  thus  both  become  permanently  blocked. 

7.  Messages  are  sent  to  and  from  the  HOST  on  lines  equal  in  capa¬ 
city  to  an  IMP-IMP  line. 

8.  Message  destinations  are  chosen  within  a  node  from  a  uniform 
distribution  on  the  ranaining  63  nodes. 

With  these  assumptions  the  network  was  simulated  with  y^  =  .01, 

N  *  50,  and  various  values  of  a  and  a  FE. 

B.  Observations 

The  surprising  result  of  these  simulations  was  that  eventually,  the 
network  blocked  conpletely  in  every  case  observed  for  a  >  y^  .  The 
network  in  the  case  a  *  y^  ,  did  show  a  degree  of  stability,  however, 
requiring  an  extremely  long  time  to  block  conpletely.  After  the  network 
had  blocked  conpletely,  an  inspection  of  the  contents  of  the  nodes  showed 
that  each  was  filled  with  messages  destined  for  the  other  IMPs,  i.e., 
they  contained  no  HOST  messages.  An  explanation  and  model  for  this  be¬ 
havior  and  the  complete  blocking  of  the  network  is  given  in  the  next 
section. 

C.  Derivation  of  the  Modified  Network  Model 

The  basic  reason  that  the  IM’s  become  conpletely  filled  with  mes¬ 
sages  for  the  other  IMPS  can  be  stated  very  siitply.  In  a  non-blocking 
network  an  equilibriun  exists  between  the  input-output  rates  (and  the 
average  storage  required)  for  both  HOST  and  ncn-HDST  traffic.  Blocking 
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causes  a  decrease  in  the  output  rate  of  non-HOST  messages  while  the  in¬ 
put  of  such  messages  remains  constant.  Cn  the  other  hand,  blocking  has 
no  effect  cn  either  the  input  or  the  output  rate  of  HOST  traffic.  The 
loss  of  equilibrium  between  the  input  and  output  rates  fear  ncn-HOST  traf¬ 
fic  causes  a  gradual  increase  in  the  storage  required  for  such  traffic. 
Eventually,  the  storage  is  conpletely  taken  over  by  ncn-HOST  traffic, 
and  thus  the  rate  at  which  the  network  delivers  messages  to  destinations 
(HOSTS)  goes  to  zero. 

Vfe  now  present  a  mathematical  model  for  this  phenomenon.  Consider 
once  more  the  simplified  network  model  shown  in  Fig.  4. 

The  rate  at  which  the  system  becomes  free  is  p^/5,  which  is 
equal  to  the  average  rate  of  message  transmission  into  the  HOST.  Simi¬ 
larly,  the  rate  at  which  the  system  becomes  blocked  is  o  -  which 

is  the  excess  of  the  arrival  rate  over  the  total  service  rate.  This 
model  assumes  that  there  is  always  a  message  in  the  IMP  that  is  destined 
for  the  HOST.  In  real  networks  such  may  not  be  the  case. 

Let  P(t)  =  P  [there  is  a  message  in  the  IMP  destined  for  the  HOST 
at  time  t] .  Then  the  average  rate  of  transmission  into  the  HOST  is 
p(t)i/°Vs  and  a  better  network  model  would  be  that  shown  in  Fig.  22. 


Figure  22.  Modified  Network  Model 
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\ 

This  mocfel  yields  the  following  system  equation: 

=  Ffe(t)  to  p(0)  +  5  }p(t))  f  o  -  y(0) 

1 

where  pg(t)  =  P [system  is  in  state  B  (blocked)]. 

Before  solving  this  equation  we  must  derive  an  expression  for 
p(t) ,  which  we  do  by  employing  the  Ehrenfest  model  of  diffusion  [24] . 

We  will  mate  the  optimistic  assumption  that  the  IMP  is  conpletely  filled 
with  messages  (optimistic  because  it  increases  the  change  of  finding  a 
message  in  the  IMP  tfet  is  destined  for  the  HOST)  and  neglect  the  fact 
that  this  means  it  is  blocked!.  We  will  use  the  modified  network  model 
(Fig.  22)  to  get  the  fraction  of  blocked  nodes  given  p(t) ,  and  p(t) 

will  be!  determined  at  the  same  time  by  means  of  the  blocking  history. 

\ 

Wfe  will  then i solve  this  systan  of  equations. 

Suppose  that  we  haVe  two  barrels  labeled  HOST  (H)  and  Store-and- 
forward  (SF) .  Distributed  between  these  two  barrels  are  N  marbles 
(rres sages) .  At  random  tines*  one  or  the  other  of  these  | barrels  is  chosen 
according  to  sone  probability  law,  and  a  marble  is  taken  from  that  barrel 
(if  it  has  a  marble) .  With  some  probabiUry  the  marble  is  put  into  the 

SF  barrel  and  witll  the  complementary  probability  it  is  put^  in  the  H 

I 

barrel. 

The  state  of  the  systan  is  the  number  of  marbles  in  the  SF  barrel 

at  time  t,  or  equivalently,  the  nunber  of  storage  cells  required  for 

I  ,  i 

store-and-forward  traffic.  In  particular,  we  want  to  knew 

! 

=  p  [barrel  SF  contains  all  N  marbles] 

*rhe  interval  between  these  times  is  presumed  to  have  an  average  value 
equal  to  the  average  time  required  for  a  transmission  plus  an  arrival 
given  the  condition  of  the  network,  i.e.,  the  fraction  of  blocked  nodes. 

! 

I 


i 
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which  corresponds  to  the  case  of  a  node  with  no  traffic  deliverable  to 
its  HOST.  Then  it  follows  that 

p{t)  =  1  -  E^(t)  \ 

Choosing  barrel  H  and  withdrawing  a  marble  from  it  represents  the 
transmission  of  a  message  to  the  HOST.  If  there  is  a  message  to  be 
transnitted,  the  transnission  rate  is  y^/5  (more  generally  it  is 
M,y^/M,  if  there  are  M  output  lines  of  which  M^  go  to  the  HOST) . 
Choosing  barrel  SF  and  taking  a  marble  from  it  represents  the  transmis¬ 
sion  of  a  store-and-forward  message.  If  a  fraction  f  (t)  of  the  nodes 
are  blocked  at  time  t,  then  the  average  output  rate  for  store-and- 
forward  traffic  is  4/5y(0)  (1  -  f(t))  assuming  that  there  are  at  least 
four  stare-and-forward  messages  in  the  IM5  and  all  of  the  output  lines 
to  other  IMPs  cure  being  utilized.  The  total  output  rate  from  the  IMP 


is  thus 


i^  +  V0)a-£(t))  -  y#)«t) 


The  probability  of  choosing  barrel  H  given  that  there  is  a  marble  in  H 

x 

arri  at  least  four  marbles  in  SF  is  thus 

1 


arri  the  probability  of  choosing  barrel  SF  under  the  same  conditions  is 

i.(°) 

5^  (1  -  f(t)) 

1  y™  -|y(W£<« 


For  the  case  of  j  SF  messages  in  the  IMP  with  j  <  4,  the  output 


i 
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rate  for  SF  traffic  is  j/4  times  the  output  rate  for  j  =  4.  The  total 
output  rate  and  the  probability  of  choosing  a  barrel  nust  then  be  adjust¬ 
ed. 

Assure  that  the  average  path  length  in  the  network  is  L.  Then,  on 
the  average,  a  message  visits  L  +  1  IM>s  in  making  its  way  through 
the  network.  If  the  time  spent  in  any  segment  of  the  path  is  approxi¬ 
mately  the  same  for  all  segments,  then  the  probability  of  a  message  be¬ 
ing  in  any  particular  segment  of  its  path  is  1/(L  +  1) .  In  particular, 
the  probability  that  a  message  is  in  its  fined  path  segment  is  1/(L  +  1) . 

Placing  a  marble  into  barrel  H  represents  the  arrival  of  an  IMP  of 
a  message  that  is  destined  for  the  HOST,  which  occurs  with  probability 
1/(L  +  1) .  Similarly,  placing  a  marble  into  barrel  SF  represents  the 
arrival  of  a  s tore-and- forward  type  message,  and  this  event  occurs  with 
probability  I/(L  +  U  • 

Let  us  define 

P„.(t)  =  P  [HOST  type  message  arrival]  =  P  [placing  a  marble  in 
barrel  H] 

P  (t)  =  P  [ store-and- forward  type  message  arrival] 

On 

=  P  [placing  a  marble  in  barrel  SF] 

POT(t)  =  P [message  transmission  to  HOST]  =  P [taking  a  marble  from 
nr 

barrel  H] 

p__  (t)  =  P  [message  transmission  to  another  IMP] 

ST 

=  p  [taking  a  rrarble  from  barrel  SF] 

E.  =  event  that  there  are  j  marbles  in  barrel.  SF 
1 

aij(t)  =  P[going  from  to  Ej  in  one  step,  i.e.,  one  message 
transmission  plus  one  message  arrival] 
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If  there  are  no  store-and-farward  messages  in  the  IMP,  then  the  proba¬ 
bility  of  a  transmission  to  the  HOST  is  one,  and  if  the  IMP  is  ccnpletely 
filled  with  store-and- forward  messages,  the  probability  of  a  store-and- 
forward  transmission  is  one.  Analogously,  we  are  not  allowed  to  choose 
an  empty  barrel  from  which  to  withdraw  a  marble.  As  a  result  we  get 
the  following: 

a01(t>  =  p»(t) 

a00^  * 

aj^(t)  =  +  Pg^(t)Pgp(t) 

ajj-i(t)  *  Par^PHA10 
ajj+i<t>  "  Pjfr(t>PsA(t) 
aNN-l(t)  "  %A(t) 
aKN(t)  ’  PsA(t) 

where,  for  sinplicity,  we  have  not  listed  all  of  the  cases  for 

i  or  j  <  4. 

Let 

Mt)  =  [aij  <  t)  ] 

Pj  (t)  =  P  [E j  at  time  t) 
and 

P(t)  =  [pQ(t)  ,p1(t)  ,p2(t)  ... 

then 

P(t  +  At)  =  P  (t)A(t) 


We  have  assured  that  the  IMP  is  ccnpletely  filled  with  messages;  there¬ 
fore,  we  must  have  a  message  departure  before  we  can  a llcw  a  message 
arrival.  Thus,  given  the  fraction  of  blocked  nodes  in  the  network 
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f (t) ,  and  a  message  arrival  rate  to  the  IM5  of  a  message/ sec.,  we  have 
that  the  average  time  required  for  one  step  (1  departure  +  1  arrival)  is 


At 


f(t) 


(44) 


The  other  equations  catprising  the  system  of  equations  that  must  be 
solved  to  get  p(t)  are  the  following: 

df  (t)  ■  -  f(t)  (a  -  y(0)  +  }p(t))  +  a  -  y(0)  (45) 

“3t  5 

P(t  +  At)  *  P(t)A(t) 

p(t.  +  At)  *  1  -  Pjj(t  +  At)  (46) 

lb  actually  calculate  the  solution  to  this  system  of  equations  we 
must  be  given  the  initial  values  pi(0)  and  f  (0) .  Equation  (45)  is 
integrated  step  by  step  using  a  value  of  p(t)  that  remains  constant 
for  a  length  of  time  At  given  by  Eq.  (44)  whereupon  it  is  recalculated 
using  Eq.  (46)  with  the  new  values  of  ai;.  (t) . 

The  solution  of  this  set  of  equations  shews  that  the  fraction  of 
blocked  nodes  changes  very  slowly  and  the  final  value  is  hi^ier  than 
that  preducted  by  the  unmodified  network  model  (Fig.  4) .  If  state  N  is 

an  absorbing  state,  i.e. ,  once  the  IM5  becomes  filled  with  store- 
and- forward  messages  it  remains  in  that  state,  the  model  predicts  that 
the  network  blocks  completely  for  the  case  a  >  y(0)  with  probability 
one. 

D.  Comparison  to  Simulation 

The  modified  network  model  predicts  that  the  case  a  -  y^  should 
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be  stable,  i.e.  should  not  block  oorrpletely.  This  is  a  weakness  in  the 
model.  By  the  clumping  analysis  we  know  that  the  equilibrium  fraction 
of  blocked  nodes  in  a  network  with  these  parameters  and  N  =  50  should 
be  about  .07.  Any  amount  of  blocking  will  result  in  a  loss  of  equili¬ 
brium  between  the  input  and  output  rates  of  store-and- forward  traffic 
and  thus  we  expect  ocnplete  network  blocking  to  be  the  final  result. 

In  one  simulation  run  with  c  =  y^  and  N  =  50  the  network 
stayed  in  the  range  3.1%  to  12.5%  blocked  for  98,000  time  units.  This 
nay  be  catpared  with  a  time  of  2,000  units  which  was  the  time  required 
bo  reach  equilibrium  in  the  unmodified  network  model  (Eq.  (18))  for  this 
set  of  parameters.  This  message  transfer  simulation  required  a  net  time 
of  250,000  time  units  to  block  oonpletely,  the  net  time  being  the  amount 
of  time  from  first  observed  blocking  until  the  entire  net  is  blockod. 

A  subsequent  run  required  130,000  time  units  to  block  oonpletely.  Both 
of  these  runs  used  a  value  of  1,000  for  a^,.  The  effect  of  this  value 
was  almost  to  insure  that  when  an  IMP  becomes  free  its  empty  spot  gets 
filled  with  a  message  frcm  another  IMP,  which  may  be  a  HOST  message. 

This  tends  to  free  the  net.  When  was  decreased  to  a  value  of 

.002,  a  rate  comparable  to  that  at  which  messages  are  arriving  frcm  the 
HOST,  the  net  time  to  total  blocking  dropped  to  91,000  time  units.  A 
further  decrease  of  to  10-®  caused  this  net  time  to  drop  to 

66,000  time  units. 

A  simulation  run  with  o  =  .01067,  y^  =  .01,  and  =  1,000 

(Fig.  23)  again  showed  seme  stability.  The  net  time  to  complete  blocking 
was  observed  to  be  118,000  time  units.  Reduction  of  to  ,002 

(Fig.  24)  and  then  to  10_®  (Fig.  25)  caused  the  net  time  to  drop  to 
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64,000  and  52,000  time  units  respectively.  The  predicted  time  to  com¬ 
plete  blocking  from  the  modified  network  model  with  initial  conditions 
P4q(0)  =  1  and  f (0)  =0  is  77,000  time  units.  In  Fig.  25  the  predic¬ 
tion  from  the  modified  network  model  assumes  P  [  store-and-f  orward  message 
arrival]  =  1,  and  P[H0ST  arrival]  *  0  and  the  same  initial  conditions 
as  before.  Okie  of  the  reasons  that  the  fit  between  the  simulations  and 
the  predicted  trajectories  is  not  better  is  the  difficulty  of  achieving 
uniform  initial  conditions  for  the  simulated  network,  which  are  assumed 
in  the  modified  network  model. 

Simulation  results  for  the  network  with  a  =  .02,  y  ^  =  .01,  and 
Opg  equal  successively  to  1,000,  .002,  and  10~®  yielded  net  times  to 
total  blocking  of  46,000;  22,000;  and  24,000  time  units  respectively. 

The  prediction  from  the  modified  network  model  is  18,000  time  units. 

We  see  that  this  model  is  far  from  being  perfect,  but  it  does  pro- 

f 

vide  nearly  quantitative  and  certainly  qualitative  understanding  of  the 
behavior  of  these  simulated  networks. 
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CHAPTER  6 


OONCLUSICNS 

A  nuntoer  of  new  models  that  have  application  to  store-and-forward 
ccnnunication  networks  have  been  presented. 

First,  we  have  the  probabilistic  model  for  nodal  blocking  due  to 
finite  storage  space  (Fig.  2  and  Eqs.  (3-7)) .  The  model  is  applicable 
when  the  average  message  arrival  rate  o  equals  or  exceeds  the  average 
message  service  rate  y(0) .  The  model  shows  that  the  blocking  behavior 
of  an  IMP  is  approximately  a  two-state  Markov  process. 

Our  second  nodel  is  for  the  fraction  of  blocked  nodes  in  a  network 
of  such  nodes  and  also  has  a  two-state  Markov  process  representation 
(Fig.  4  and  Eq.  (18) ) .  The  result  appears  valid  for  both  randomly  con¬ 
nected  and  lattice  networks  and  for  a  variety  of  system  parameters  (Figs. 
9-14) .  However,  the  model  for  the  fraction  of  blocked  nodes  in  a  "par¬ 
tial  blocking”  network  (Figs.  15-17)  needs  to  be  greatly  improved. 

Various  clumping  models  have  been  presented  and  shown  useful  for 

(0) 

such  a  network.  The  clump  size  distribution  for  the  case  a  =  y 
(Eq.  (30)  and  Fig.  18)  and  the  rtaxintun  clunp  size  model  (Eq.  (38)  and 
Figs.  20  and  21)  appear  adequate  to  describe  these  cases.  The  average 
clump  size  for  the  case  a  >  y(0)  (Eqs.  (35-37))  is  a  fair  apprcximation 
and  needs  further  work. 

The  modified  network  model  (Fig.  22)  provides  a  clue  to  the  funda¬ 
mental  behavior  of  store-and-forward  ccmunication  networks  that  are 
subject  to  overutilization.  The  model  treats  the  case  a  >  y^  fairly 

Preceding  page  blank 
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well  (Figs.  23-25)  but  does  not  appear  applicable  in  the  marginal  car*-! 


Further  work  on  models  of  this  type  appears  justified.  An  effort 
should  be  made  to  improve  the  modified  network  model  for  the  case 
o  =  ,  and  investigations  should  be  made  into  the  transient  clunping 

behavior  in  ocnpletely  blocking  networks.  Also,  the  variance  of  the 
measurement  of  the  fraction  blocked  in  such  networks,  and  the  time  depen¬ 
dent  connectivity  requires  investigation. 

Questions  regarding  the  behavior  of  networks  with  selective  blocking, 
as  in  the  ARPA  network,  remain  unanswered;  nor  have  we  introduced  the 
effect  of  multi-packet  messages.  These  would  be  important  (and  difficult) 
areas  for  research. 

The  whole  subject  of  blocking  in  networks  of  this  type  appears  to 
be  absent  from  the  literature.  We  believe  that  this  field  contains  many 
jvMtHrmal  challenging  research  areas. 
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APPENDIX 


A.  Solution  of  P  =  AP  +  C  for  sane  special  cases 
In  this  section  we  solve  the  network  equation 

P  =  AP  +  C 

for  some  special,  network  topologies.  Recall  that  if  there  are  m  nodes 

til 

in  the  net,  then  P(t)  is  the  m  x  1  matrix  whose  i  ccnponent  is 
the  probability  that  node  i  is  blocked  at  time  t.  A  is  an  m  x  m 
constant  matrix  and  C  is  an  m  x  1  constant  matrix.  The  solution  is 

P(t)  =  eAtP(0)  +  A-"1  (eAt  -  DC 

'Thus  our  problem  is  to  find  the  exponential  and  the  inverse  of  the 
matrix  A. 

1.  Lattice 

2 

Consider  a  network  consisting  of  m  =  n  nodes  arranged  in  an 

n  x  n  grid  with  4-neighbor  connections .  For  this  system  the  matrix  A 
2  2 

is  n  x  n  and  takes  the  following  form: 


A  = 


D  A 
ADA 


O 


O 


ADA 


A  D 
A  D 
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b  a  b 

where  D  =  b  a  b 


b  a  b 


b  a  n  x  n 


A  =  bl 


where 


a  =  -a,  b  =  —  ,  and  1^  is  the  n  x  n  identity  matrix. 

We  must  first  find  the  eigenvalues  yv  of  D  which  are  the  solu¬ 
tions  of  | D  —  yl |  =  0.  Let  a  stand  for  a  -  y  in  D;  we  wish  to  find 
the  zeros  of  the  determinant  of  D.  Expanding  by  the  elements  of  the 
top  row,  we  note  the  following  recurrence  relation  for  the  determinant 
A^  of  the  n  x  n  matrix  D: 

m  aVl  '  b\i-2 

with  initial  conditions  A^  =  a,  AQ  =  1,  A^  =  0.  Following  Grenander 
and  Szego  [25]  we  substitute  a  =  2b  cos  0,  assune  a  solution  of  the 
form  Ah  =  pn,  and  solve  the  resulting  quadratic  in  p.  After  satisfy¬ 
ing  the  initial  conditions  the  result  is  sirrply 

,  _  ,_n  sin(n  +  1)6 

*n  =  b  — sGTT” 


which  vanishes  for 


6  =  V7r/n+l  v  =  1,2,  n 


Therefore,  the  eigenvalues  of  D  are 
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a  -  2b  cos  n  +  ^  v  =  1,2,  n 
which  are  all  distinct.  The  eigenvectors  are  the  solutions  of 


a  b 

*vl 

V 

b  a  b  C3 

Xv2 

Xv2 

b  a  b 

\3 

= 

*v3 

1 

o 

• 

• 

• 

or 

•  •  • 

-V 

•  •  • 

-V 

It  is  easy  to  verify  that  the  normalized  solutions  are 

kvrr 

V  -  - •  sm - 


n  +  1 


so  that  the  (i,j)  element  of  e 


D  _  ^  JV  y  y 

i'j  ‘  5l  ^  ^ 


and 


■  2  «\TX  \i  Vi 


,-l 


where 


w 


yv  =  a  -  2b  cos  j-j 


and 


.  (-Dn'k  J22L. 

n  +  1 


V  =  J"^=  sin 


Vn  +  1 
2 


Similarly,  it  is  easy  to  show  that  the  transformation  R*AR  (where  R* 
is  the  transpose  of  R)  where 
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XllIn  *'•  Vn 


.  X  ,1 
nl  n 


X12In  •“  Xv2In 


X,  I  ...  XI 
In  n  vn  n 


•  Xn2Xn 


•  Xnr>In 

nn  n 


with  X^.  as  given 


above  reduces  A  to  the  quasi-diagonal  form 


»!  O 

«2 

— ^  •  •  • 
O  "n 


where 


=  D  -  2b  oos  In 


Since  M  is  equal  to  D  with  a  change  of  the  diagonal  element,  we  have 

A 

that  the  (k,l)  element  of  the  (i,j)  block  of  e  is 


ei,j;k,l  -  I.Vvj  co3m  -»  “^'Vpl 

*LW  *  l/vi^  £<•-»  - 213  “^r>"Vpi 


where 


V  -  <-Dn~k  ;in 

**  cn  n  + 1 


In  our  system  a  =  -c  and  b  *  u^/5  so  the  time  constants,  i.e., 
the  arguments  in  each  of  the  exponentials  appearing  in  the  solution  for 


e  are  of  the  form 
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4  (0)  it  .. 

exp  -  (a  -  5  y  °°s  n  + 

Hie  nuntoer  n  is  the  square  root  of  the  nurrber  of  nodes  in  the  square 
lattice.  Hus  result  shows  that  as  n  +  “  the  system  attains  its 
steady  state  at  a  rate 

exp  -  (o  -  y  )  t 

which  agrees  with  sinulation  results  for  n  =  32  (see  Eq.  (18)  and  Figs. 
9-11) . 


2.  Torus 

Again  we  consider  a  network  of  m  =  n2  nodes  arranged  in  an  n  x  n 
grid  with  4-neighbor  connections,  but  this  time  we  assume  that  opposite 
sides  of  the  grid  are  connected  together.  The  result  is  a  torus,  and 
for  this  case  the  matrix  A,  again  n2  x  n2,  takes  the  follcwi.ig  form: 
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b  a  b 


where  D  = 


b  a  b 


b  alnxn 


A  =  bl 


where  a  =  -a,  b  =  y^/5,  and  I  is  the  n  x  n  identity  matrix.  The 
solution  follows  from  the  fact  that  A  is  a  block  circulant  matrix  [26] . 

It  is  easy  to  verify  that  the  transformation  R*AR  (where  R*  is 
the  transpose  of  R)  where 

IWn  ¥n  Vn  1 


R  E  xllXn  *  *  *  XilIn  ’ • *  XnlIn 


Xln-lIn  Xin-lIn  Xnn-lIn 


xik= 

V^si"T 


k  —  0,1,  • • ,  n-1 


v  even,  ^  0,  <  n 
k  —  0,1,  •*,  n-l 


j — l—  if  n  even;  k  =  0,1,  n-1 

/Aj"'  /a 


reduces  A  to  the  quasi-diagonal  form 
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o 


where 


D  +  2A 


•V-Vl  *  D  +  2A  “s  5T 

M  =  D  -  2A  if  n  even 


v  even,  ^  0,  <  n 


Therefore,  the  (i,  j)  block  of  e  is 


n  M 

£ w s 


with  N  and  X  ,  as  given  above. 

t>  bK 

We  observe  that  the  matrix  D  is  simply  matrix  A  wherein  each  block  is 
of  dimension  one;  therefore,  the  (k,l)  element  of  eD  is 


_  n  m 

=  £W  p 


=  J^VWV 
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where 


m^  =  a  +  2b 

=  n\v+1  =  a  +  2b  cos  ~  v  even,  ?  0,  <  n 
mn  =  a  -  2b  if  n  even 

Since  A  =  bl  ,  each  of  the  matrices  M  is  equal  to  D  with  a  change 
n  s 

of  the  diagonal  element.  Hence  the  (k,l)  element  of  the  (i,j)  block 
of  eA  is 


■"is  -  as  +  25 

\s  -  Vis  ■  as  *  21  003  r  V  even,  f  0,  <  n 
=  ag  -  2b  if  n  even 


a1  =  a  +  2b 

ar  =  ar+1  -  a  +  2b  oos  ~  r  even,  7*  0,  <  n 
a  =  a  -  2b  if  n  even 


and  with  as  given  before. 

3.  Twisted  Torus 

2 

Cnee  more  we  consider  a  network  of  m  =  n  nodes  arranged  in  an 


n  x  n  grid  with  4 -neighbor  connections.  We  assure  that  nodes  are  nun- 
bered  from  1  to  n2  by  rows  and  that  node  i  is  connected  to  nodes 
i  ±  1,  i  ±  n  modulo  n2.  A  "twisted  torus"  results  for  which  the  connec¬ 
tion  me  xix  A  is  as  follows: 


The  matrix  is  a  circulant,  i.e. ,  any  row  is  a  cyclic  shift  of  the  pre¬ 
vious  row.  Following  Bellman  [26]  we  find  that  the  eigenvalues  are 

/  2Tiki  2iTtau  2TT|(n2_n)i 

=  a  +  b\e  n+en+en  +en  / 

=  a  +  2b (cos  ^  +  cos  k  =  0,1,..,  n2-l 

n 


where  i  =  /-I,  and  with  associated  eigenvectors 
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1 


2rrk  , 

“71 


^(n2-l)i 

n 


Ihe  eigenvalues  occur  in  pairs  in  all  but  the  extrene  cases.  This  can 
be  shown  easily: 


.  .  _ 2rr(n2  -  k)  _ _ 2Tr(n2  -  k)x 

y  =  a  +  2b  (cos - —5 - +  00s - — - ) 

n-k  n 


n 


=  a  +  2b  (cos  ^  +  cos 

n2  n 


=  Yt 


Thus  the  eigenvalues  are 


Y^  =  a  +  4b 


Xi^iA 


V  -  ft  lw» 

*vk  -  *n  sln  ~Z 
n 


„  _ kVTT. 

Vlk  =  n2 


v  even,  ^  0,  <  n 
k  =  0,  1,  n2-l 


/-l)*  2 

X  ~  =  v  '  n  even;  k  =  0,  1,  n  -1 
n\  n 


therefore,  for  the  twisted  torus 
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3-2)  0  0  0  0  0 


0  X  X  X  0 


0  0  0  0  0 


0  0  0 


0  0  0  X0 


0  X  X  0  0 


0  0  0  0 


0  0  0 


0  0X0 


0  0X00 


0X00 


0  0  0 


ooooo 


0X0X0 


0  0X00 


0  0  0 


0  0  0  0 


0  X  X  0 


0  X  X  0 


0  0  0  0 


4-2)  0  0  0  0  0  0 


0  X  X  X  X  0 


0  0  0  0  0  0 


,<3>  +  4X(2)  +  6X(1) 


y3  =  y(2)+*(1) 


X3  =  X(3)  +  4X(2)  +  9X(1) 


=  U<2>  +  2y(1) 


X3  =  4X(2)  +  12X(1) 

y3  -  y(2)  +  2y<1> 


X3  *  X(3)  +  3X(2)  +  12X (1) 

y3  =  y(2)  +  2y(1> 


X4  =  8X(2)  4X(1) 

y4  -  4p<3> 


X4  =  4X(3)  +  4X(2)  +  6X(1) 

Va  -  2y(2,.+  2y(1) 
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\ 


\  • 


4-8) 


0 

0 

0 


0  0  0 

i  • 

X  0  0 

O  X  X 

I 

0  0  0 


0 

X 

0 

0 


0 

0 

0 


X4  =  2A{3)  +  4A(2)  +  12A(1) 
y.  =  2y(2)  +  2y(1)  . 


4-9) 


0 

0 

O 


0  0 
X  0 
O  X 
0  0 


0 

X 

0 

0 


0 

0 

X 

0 


0 

0 

0 


X4  *  2A(3)  +  4A(2)  +  12A(1) 
y4  =  2y<2>  +  2y<» 


4-10) 


O 

O 

O 


O  0 
O  X 
O  X 
X  o 

o  o 


0 

0 

o 

X 

o 


0 

0 

0 


A.  =  3A(3)  +  3A(2)  +  11A(1) 
4 


4 


i<3>  +  3U<1> 


4-11) 


0  0  0 

0X0 
O  0  X 

0  X 
0  0 


o 

X 

0 

0 


0 

0 

0 


A4  =  3A(3)  +  2A(2)  +  llA^a) 
y,  =  y<3>  +  2m«>  +  y(1) 


The  approximation  far  large  n  is  again  a  straight)  line  topology.  In 
simulations  on  an  8-neighbor  lattice  round  clumps  are  observed  mare 
frequently  than  stringy  ones,  but  for  clunp  topologies  up  to  4,  the 
straight  line  clutp  yields  better  approximate  growth  and  death  rates 

l 

than  does,  say,  a  square  clunp.  The  approximation  gives  ' 
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1 


\ 

1 


A  =  2(n  -  2)A(3)  +  4A(2)  +  6X(1)  =  (6n  +  2)X(1) 
n 


-  7n  +  2,  (0) 

9  M 

Wte  f  ind  that  a  better  approximation  for  n  =  4  (and  thus  we  will  assume 


for  larger  clumps  as  well)  is 


Xn  =  era*1’ 


u  _  7h  tfO) 
n  9 


For  the  stationary  probability  of  a  clunp  of  size  n  we  then  get 


where 


,(°)  n-1 

pn'p0^ir~  "-1 


r.  X(0)  JnU  -  r)  1 

r?- ~  J 


„  54  X*11 

Z-T-J5T 


X(°)  PQ 

E(#  in  system]  =  ^  _  ^.y 


For  the  case  a  =  ,  N  =  50  this  yields 


E[#  in  system]  =  .1338 


The  result  observed  in  simulations  is  about  .10. 
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The  following  programs  performed  the  simulations  described  earlier. 
They  are  written  in  Fortran  IVH  and  run  on  the  XTS  Sigma-7  conputer  at 


U.C.L.A. 


1.  Lattice  (with  graphical  display) 


2.  Random  Graph 

3.  Message  Transfer  Network 
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Lattice 


1 

2 

3 

4 

5 
4 

7 

8 
9 

10 

11 

12 

13 

14 

15 
14 
17 
II 

19 

20 
21 
22 

23 

24 

25 
24 
27 
21 

29 

30 

31 

32 

33 

34 

35 
34 
37 
31 

39 

40 

41 

42 

43 

44 

45 
44 

47 

48 

49 
5C 


C  MAIN  PROGRAM 

CRMM8N  NSETI22<3P<KI<  PARAM ( 10. ID) <  JRATE C 10) <  T-f  I N C 1 0 1  •  ATRIBI8I# 

1  Kh  ITE<  JRUN<  KRUNT (10)#  ISEEO,  TNBW<  NBD\<  NRSECl lh) j  PCSl C(16>< 

2  NfiANDl  12AI#  I.RAND#  MRANO#  KRANu#  KriflBO<  KBAD<  PARAM?! 1 j< ID  1 
nif'MON  r.rt»0F«5000).  JAR|J2#32),  1nT1(?)<  INT2I2) 

0 1  MENS  1  BN  JINT1(2)«  JINT2<2) 

IN (EGER  ATRIb#  TNR*<  TE  IN 
INTEGER  GCBOE 
DATA  IGBBJ#  I BAD/1  *»»8'/ 

KCittBO*  I  3M6D 
«HA0«1BAQ 

CALL  DS6PEM GCBOE #5000) 

DATA  wINTKl  )<J1NT1  IPI/'OKOO'dlOJA'/ 

DATA  wlNTPI 1 )< JlNTP(P)/' JaOO*  <  MOOF •/ 

INTI  I l)aJlNTl(l) 

INTI (2)aJlNTl(2l 
I NT2( 1 1* JINT21 1 1 
INTr<2)aJlNT2(2> 

CALL  IfEGPlCl ’PICTURE* ) 

C  BEGIN  PICTURE  OFF  INI T IAN 
»X*l53 
17*950 

Ob  300  I  at <32 
Jal 

K* INTENSI ? ) 

CALL  OOSTP(IFT) 

C  IPT  IS  THE  CURRENT  STACK  POINTER 
r.  INTENSITY  a  2  IS  INTENSITY  ME  GBRD  NBDFN 

C  INTENSITY  *  7  IS  INTENSITY  »F  BAD  NUDES 

JAM  I#  Jl  ■  IPT 

C  JAR ( 1 1 J ) ■ ABS  ADDRESS  BE  INTENSITY  INSTR.  EAR  PRINT  t*J> 

KaPCilNTl IX< IV) 

OB  200  J*2<32 
KaINTENS(2> 

CALL  ODSTPIIPT) 

JAM  I#  Jl  ■  IPT 
KaRPfilNT ( 24<0 I 
200  CBNTlNUE 
IX* 153 
IYaIY-24 
300  CBNTINUF 
CALL  ENRPIC 

CALL  UDPAOl 'PICTURE •< IBP) 

OB  31  I  at  #  3? 

OB  31  Jal*  3? 

31  JAM !<J|aJAR( l«J)-TMP4l 

C  JAR(I<J)*REL  ADDRESS  RE  INTENSITY  INSTR*  FRR  PRINT  (IfJI 
CALL  0 ISP  1C 
C  DISPLAY  PICTURE 

RE  ADI  IDS# 143)  LRANO<MRAN0<KRAND< NRANOI 1 ) 
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51  1*3  FNRMAT (419) 

58  DA  14*  JaPs 12A 

53  1**  NRAND(J)aNRANO(J.lU?aj 

5*  JRUM1 

55  5  CALL  DATAN 

54  CALL  GASP 

57  ON  3  l«1* J2 

58  Oft  3  w»  1 »  J2 

59  3  CALL  DDfi£P( •PtCTUHF * # JAR( I# J) # INTI  I 

6C  IF  { JKUN  *GE  t  10)  Qfl  TS  57 

41  C  LAST  RUN  IS  NB  10*  THIS  CARO  CAN  BE  ALTERED  TA  ALLBm  AW  NUMBER 

42  C  OF  RUNS  UP  TA  10  PITH  DIFFERENT  PARAMETERS  FBR  EACH  RUN* 

43  ,)RUNa«jRUN*l 

4*  GA  TB  5 

45  57  ST6P 

44  ENL 

47  r 

48  C 

49  SUbKBUT  INS  ORANO  (RNUM ) 

70  COMMON  NSET ( 32* i?i A ) •  PARAMI 10* 10) » JRATE ( 10 > »  TFlNiiO)*  ATRIBI8)* 

71  1  Kk I TE«  JHUN,  KAUNTMOI*  ISFEDt  TRBUt  NBDN,  N3SEC<14).  PfSEC(16)A 

72  9  NKaNOI 128) i  I.RAND*  MRAnD*  KRANJ 

73  LRANCaLRAN0a45ST9 

7*  MNAN0aMRAI*0a33&5**31 

75  Jalal aBSI lRaND) /14777214 

74  RNUna  . SaFLBAT (NRAnD( J) +LRAND+MRANU) a  •232A3064E-9 

77  KRAN0aKRAN0a342«360A9 

78  NRANUlJIaKRANO 

79  RFTURN 

80  FNC 

81  C 

82  C 

83  SUBROUTINE  BNDER  (M,N) 

8*  ent'KBN  NSET ( 3?» 3P* A )*  PAKAMI 10* 1 0) » JRATE ( 10 ) *  TFtN(lO)*  ATRIB(S)* 

85  I  KKITt#  JKUN*  NAUVT  (10)#  (SEED*  TNBW*  NBDN*  NRSEC  ( 14 )  •  PCSECC14) 

84  INTEGER  ATS*  ATfc*  AT 7*  ATS 

87  InTEGnH  ATRIB*  TNB*.  tfin 

88  C  AHDFN  ADDS  NflDC  ( M«N )  TS  BRUEREO  LIST  BF  EVENT  TIMES  FAR  DaTAN 

89  Mhao 

90  M7«0 

91  C  M5  AND  M7  ARE  FLAMS*  taHEN  BbTMai  wE'VE  FOUND  RIGHT  SPBT  FAR  ( M»N) 

92  fal 

93  >1*1 

9*  C  KNOW  THAT  NODE  (1*1)  HAS  BEEN  ORDERED 

95  IE  (KRITE  *E0*  A)  OA  TB  1* 

94  laATRIB(S) 

97  JaATRIBI 4) 

9M  IE  (J  »NF *  7777)  GA  TB  U 

99  | aATRlRI 7) 

IOC  JaATRIB(M) 
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101 

102 

103 

10* 

105 

106 
107 
106 
109 
UC 
til 
112 
113 
11* 
lib 
116 

117 

118 

119 

120 
121 
122 
123 
12* 

125 

126 

127 

128 

129 

130 

131 

132 

133 
13* 

135 

136 

137 
136 
139 
1*0 
1*1 
1*2 
1*3 
1** 
1*5 
1*6 
1*7 
1*8 
1*9 
15C 


1*  IF  |ATWIBil)-NSFTf l,j*l))  h,6*7 
C  «LT«  0  MEANS  SB  TB  PREDECESSOR*  «GT*  0  SB  TO  SUCCESSOR 

5  lAaNSET ( 1 « J«7) 

IF  |M7  ■  •£  J*  1)  SB  TO  9 
IF  (1*  .F:.  9099)  GO  TO  9 
C  SEE  IF  IT  HAS  A  PWEDFCFSSOR 
.laNSETI  I «  j*  o ) 

1  a  I A 
Mhal 

Gb  Tt)  1* 

7  UaNSETl f«J«5) 

IF  (Mb  *E0*  1)  SO  TO  6 
IF  I IA  •Eh*  77771  GO  TO  6 
C  SEF  IF  IT  HAS  A  SlICCFSSttR 
JaNSET I !« J«6) 
la|  A 
M7«  1 

no  to  i* 

6  ATn IB  15) aNSET  1 1  * J*5 ) 

C  (H,N)  SUCfFUES  II..J)  »RV  CONVENTION  IF  TFiEV  HAVE  a  TIMFS) 

ATR IB  16 ) aNSET 1 1 *U* A ) 

ATkIBI7)aI 
ATk IB (ft) aj 
NSETI  I  * t J*  b )  ah 
NSE  T(  I  *.l«  ii )  aN 
ATbaATR ld<5 ) 

IF  IAT5  •E0«  7777)  GO  TO  90 

C  TFST  FAILS  IF  (I*. i)  HAS  A  SUCCESSOR.  AND  THUS  MUST  UpDATf  HIM 
AT6aATRIB(6) 

NSl T(ATS«AT6«/)aM 
N8E  T ( ATS*  AT6*  8 ) aN 
r.n  TO  SA 
9  ATHIBl 5) a  I 

C  <M,m)  PRECEDES  (I,J) 

A1RIBI6) aj 

ATGIb I  7 )  aNSET ( I*  J* 7 ) 

ATk IB ( A ) aNSET ( I « J*  A ) 

NSETI I»J»7)aM 
NSETI  I  « .)«  H  )  aN 
AT7aATR!8I7) 

IF  IAT7  »EQ»  9999)  GO  TO  90 

C  n*'°  W*S  A  RECFRSbR  I mHICH  MUST  RE  UPDATFD) 

ATN Inis) 

NSE  T I  AT  7*  ATO*  5  )  aF* 

NSET IAT7*AT0«6)aN 
96  DO  99  Ka5*0 
99  NSET  (M*k*K)aATNIRtiK) 

RETURN 

END 
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151 

152 

153 
IS* 

155 

156 

157 

158 

159 
16C 
161 
162 
163 
16* 

165 

166 
167 
161 

169 

170 

171 

172 

173 
17* 

175 

176 

177 
176 
179 

160 
161 
182 
183 
J8* 

185 

186 

187 

188 

189 

190 

191 

192 

193 
19* 
19b 

196 

197 

198 

199 
?0C 


c 


r. 

c 


L 

C 


SUbNBUT  INE  RACK  I  J*K ) 

CftKMBN  NSETl3r*9?*M)»  PARAMC1C»10>*JRA1EI10)*  TF1NI10)*  ATRIBIB)* 
1  KR 1TE  *  J*UN*  kAIINTMO)*  ISCED*  TN9K*  KBON,  NBSECI16)*  Pr!>EC<16)* 
?  NRANPI  1281*  I  HAND#  MHAnD*  KWANU*  KGBBD*  KBAD*  PARAM2C 10*  10) 
CHhltOK  GChOEI aOOOJ*  J7RI32* 321*  1NT1I2)*  INTPI2) 

INTEGER  ATKIb*  TN9a»  TFIN 
INTtGER  GCBDE 

(J*K)  IS  AN  INITIAILV  DAO  NtlDF 
NSt TCd*K*?)»l 

CAU  DDRFPI  'PICTURE  »*jAR< J*K)»  TNT2) 

Will.  GPL'ATE  KB  RF  P*D  NSDES  IN  SFCTJR*  AND  PCT  PAD  IN  SfcCTAR 
NSECT  aNSFT ( J« K«  3 ) 

NHSECINSECT)«NBSFCINSECT)*1 
RSEC-KBSFC INSECT) 

PCSECINRFCT ) «RSFC/6*»0 
IF  IJ-3?)  2*3*? 

2  NSET ( d*1 *K«  * ) "NSFT I J+1 *K*+ ) +1 

3  IF  (J-1)  **5» * 

*  NSE Tlo-1*<**)«NRFT( J-I*<**)+1 

5  IF  IK-3?)  6*7*6 

6  NSET I d»«*1 *  * ) «NSFT IJ*K+l«*)+t 

7  IF  IK-1)  «*9*rt 
H  NSETIw*K-1«*)>NSFTIJ«K-1**)*1 

9  RULRN 
FNt 


SOcRBgT INE  FIND  I NRBR*NCBL»  NCPPF ) 

rfll.BBN  NSET I3?*3?*rt)«  PAWAMC  10*1 0 >* JRATE C 10) .  TFINI10)*  ATRIBIB)* 
1  KRlTE*  JRUN*  KBUNT(IO)*  1SFED*  TNHM*  NBDN*  NUSECI 16 )«  PCrECI 16) 
INTEGER  ATR lb*  TNfiw*  TFIN*  ALT 

TF  ||NRr,M.LT.1.BR*NRBK.3T*3?).f)R*lNCftL.LT*1*NH*NCBL.bT.3?)) 

1  OB  T6  6 
JRBMNHftrt 
JCBL+NCM. 

ATK Ibl 1 ) »NSET IJRBW*  ,JC3L*t ) 

ATR IBI?) »nSET I JHflK* JCBU«2) 

A  TR IB I  * ) »NSET I JH«R» JC9L*  * ) 

,ITIIiE»ATRIRI  1 1 
JNEhBaATR 131*1+1 

IF  INCBCF  +Ew*  0)  ATRIBI*)«ATRIBI*)+1 
IF  INCBPF  .EC..  1)  ATRIfll*)«ATR|BU)-1 
NRETIdRRR* JCBL**)«ATRI9I+) 

knerb»athibi*)+i 

JSGb«5» ATR IHI r  I+JNFH8 
KSOfaab*  ATR  IHI 2  >+KN'FHB 
CALL  CAM jR6K*jrftL*AI  AMD) 

CALL  uHANDIRNuH) 

At  TaTNB.-INTI ALAMO*A(  BGIRNjM ) -0.5) 


108 


701 

702 
203 

704 

705 
7  06 
707 

?o« 
709 
P  1C 

711 

712 

713 
214 

715 

716 

717 

718 

719 
??C 

721 

722 

723 

724 

725 

726 

727 
72  fc 
725 

730 

731 

732 

733 

734 
73b 

736 

737 
736 
739 
74C 
741 
?4it 

743 

744 
?4S 
746 
74  7 

748 

749 
75C 


3  ATkiatlltALT 
(ib  T6  6 

4  IF ( ATs 1«( 1 ) "ALT )  b«6«7 

5  t  F  ( NCOOPt  «NE  aATWIE*I711  ATR IP ( 1 1 aALT 
(irt  TS  6 

7  IP  iNLBOFaPGaATHIt  I?))  ATRIP  11 JaALT 

6  CSMlNUP 

NSETI uR6*a JCBl# 1 laATRtblll 
IP  I  jTIhEaEbaATwlDM  >>  lib  TB  in 
r.ALL  H'lfiVEURttWaJCflLl 
CALL  ttROEHI  JPIrtW.JCBLl 
10  C.HMINUE 
U  «f  TUK.m 

fnd 

c 

c 

RilbHBUT  INt  HMHVF  (JRBWa  JCHU 

f.fthnttK  KSET ( 37« 37<N  I  j  PARAM <  1  Ot  1 0 » *  JRATE ( 10 ),  TP  I N(  1 0  )<  ATRIBI Alx 
1  kHITEj  JRLlNi  KBUMUblx  IRFED*  TI4«R#  NBDN#  Nh»ECI161#  PrsEC<16> 
IKTE6P.H  ATHIbi  TNBR.  TP  IN 

on  io  Kai«a 

10  ATMBlKla  NSE T I ,IR9W ( JC BL « A ) 

I SbChaATR 16 15 ) 

I SLCCaATH IB ( 6 1 
IPRfcKaATHlB ( 7  I 
IPRLOATRtBlbl 

IF  (  IPREW  aEQa  49991  uB  Tft  15 
NSETI iPKFRa IPkECablaliUCR 
NStTlIPt'EKi  IpHtCablafSUCC 
IF  ilSbCR  .EO.  77771  GB  TB  25 
lb  NSL I ( ISUfHj I SbCC«7 leTPRER 
HRETI ISlJCRx ISuCCaB 1 aIPKEC 
2b  RPToKN 
FND 
C 
c 

Rl.bkllbTlNP.  HASP 

f  rihhtJN  NSET  (37#  37*^  l«  PARaHl  ic#  1 0>  *  JHATEI 101 1  TFIMlOla  ATRIBI81# 
1  KMTEa  J*bM#  ABUNT  (1  Ola  ISP  ED#  T\BW*  NBDN#  NHSECI16I#  PrsECIl6> 
IKTLbl  k  ATKIb#  TNt.W,  TF  IN*  SUCR*  SUCC 
M  XTPial 
NP.ATLal 
JRITtaO 

4  IP  (NafcT  <NEXTR«NFXTC»7>  afcO.  99991  lid  TB  5 
C  TFST  AbtJVf  FAILS  IF  IAFxTR#HEXTCl  IS  NBT  THE  HFAI)  BF  THF  LIST 
NPXaNSP  T(  ('if.XTH*  NFATCa  7 1 
NFxTCaKSF.T<NE*TH.PFxTC#91 
NPXTPiaMX 

C  SFF  IF  HIS  PHP.Of  CFSSBR  IS  HPAO  MF  LIST 
Gh  Tt  4 


251 

s 

f.ALL  ftMeVE(NEkTR.NFxTC) 

25c 

ITEST.ISITBN(B) 

253 

IFI ITEST.nE.S)  OR  T8  59S 

PS* 

597 

ITEgT.ISITOMi*) 

?  56 

IFI 1TEST.E0.0)  RETURN 

256 

GH  Tb  59/ 

257 

598 

ChNTlNUE 

258 

TNB»«ATMH(1) 

259 

IF  (TNOw  tGE.  TFIMJRJN))  G8  Tft  17 

26G 

SUCR.ATRJ«(5) 

261 

SUCCaATRlt>(6) 

P62 

JWlTEaJH)TE*l 

263 

HKIU*0 

P64 

IF  IJRITE  .IT.  JRATF I JRUN 1 1  GB  TO  7 

P65 

r. 

IF  AbbVF.  TFST  IS  MFT.  MILL  NOT  PRINT  RFSULTS  AT  THIS  EVFNT  TIME 

266 

JRlTiaO 

267 

KR I TL»1 

268 

r 

WILL  PRINT  RESULTS  AT  THIS  EVENT  TIME 

269 

7 

CALL  EVNTSINEkTR.NF  XTC I 

27C 

NFkTR.SI.CH 

271 

NFKTC.Sl.CC 

272 

GH  Tb  4 

273 

17 

KRITE-I 

274 

URlTt  (108.211 

97b 

21 

FORMAT  (lH0.52k.26H*.  FINAL  MEPHRI  FOLLOWS  \ 

276 

CALL  EVNTS  (NEXTM.NFXTC) 

977 

RETURN 

97 a 

FNU 

979 

c 

PH  0 

c 

281 

SUbRbUT  |.Nt  OAT  AN 

282 

CbHMBN  NSFT(32« 92.81*  PARAM (10. 10>. JR ATE (10).  TFIN(IO).  ATRIBl 

9M3 

1  HR  I  Tt*  JRUN.  KAUNT  ( 10) «  ISEEO.  TNOW.  NBDN.  NHSECM6I.  PCSECIl 

9  »ts 

2  NKkNC ( 128 ) .  1  RANO*  MRAND,  KRANU.  KGBOO.  KHAO.  PAHAM? ( 1j, 10 ) 

99  5 

OIMEN&IAN  MAP (64 I .  BAD(2j) 

Pit 

INTEGER  ATRIB*  TNOW,  TFIN.  BAD 

997 

real  lambda,  mu 

999 

IF  (  jHUN  .NE.  1)  tlA  TO  29 

299 

READ  (105.31)  KRITF 

290 

31 

FARMAT  (ID 

291 

IF  (kRITE)  32.32.33 

292 

3c 

op  19  j.i.in 

293 

REAU  (105.21)  LAMBDA,  MU.  N 

294 

WRITE  (108*10)  1  AMHOA.  Mu.  N 

29b 

IF  (LaMBOA-MU)  6.7.8 

296 

6 

WHITE  (108*9) 

297 

9 

FhWMAT  (50H0ERHAR-  LAhBOA  MUST  RE  GREATER  THAN  OR  FOUAL  Ts  MU) 

999 

J9LN.20 

299 

RFTuRN 

300 

/ 

P.N 

110 


.101 

102 

303 

104 

305 
106 
307 

306 
309 
31C 
311 
31? 
313 

114 

115 

316 

317 

116 
319 
32C 
121 
322 
123 

324 

325 

326 

327 
326 
329 
33C 

331 

332 

333 
.134 

335 

336 
.137 
336 
339 
34C 

341 

342 

343 

344 

345 

346 
34  7 

348 

349 
35C 


PAKaM( l,J)aP/MU 
c»f»  T8  11 

6  PARAM1 1 *J)*1*C/M  AMHDA-MU) 

11  PARAM(6« J) ■ 1 • C/MU 
DH.TAaMC/5«0 
08  19  I  a?#  5 
MlJaMU-OFLTA 

PAKAM(I*J)  •  1  •  C/  ( L  AMflDA"fiU  ) 

19  PAR AH( 1*5* J ) a  1  *0/Mli 
tu  FORMAT  ( MHOLArlBDA**  F7*  4*  5»*  .IHMUa 
1  22MS1/E  BA  UllEllETNG  Rt>8M**  15 

21  FORMAT  <?(F7.C)*tb) 

3a  OH  34  1*1*10 

34  RF  AO  (105*55)  (PAHAM21 J* I )* J*l, 10) 

55  FORMAT  ( 10E7*2) 

DO  37  I  a  1  *  1 0 
Oh  37  Jal* 10 

37  PARAM2I  I*  j ) a  1  ,/PARAM? (  1*  J) 

35  FORMAT ( 1UF7*0) 

36  RA  AD  (105*22)  ( JRATF( 1 >« 1. 1* 10) 

22  FORMAT  ( 1017) 

Rt  AO  ( 1  Otj*  23 )  ( TF I N ( I  ) «  I  •  1*  1 0 ) 

23  FiRmAT  (1017) 

08  4  wa 1  *  10 

WRITE.  (104*3)  .)*  JRATE(J)*  TFIN(J) 

4  WRITE  (108*5)  (PARAhf I*J), 1*1,10) 

5  FtihhAT  <1H  *115PARAHFTtRSa*l0(F7.4*5X) ) 

3  FhRMAT  (1H0*7hRun  N(|.*I2*10X*l?HRtPoRT  RATF  a*  1 7*  10X* 

1  12HFIMSH  TIMF a*  |  7) 

29  RFAti  (105*25)  MUM 
tSttOao 
KHlTE*6 

25  FORMAT  ( T 4 ) 

C  NIR\  IS  Nh  »F  IMTTAILV  BAD  NBDFS 
NBCNaNIRM 

RA  AO  (105*36)  KOUNTM) 

36  FhRMAT  < II ) 

WRITE  ( 1  OR*  26  )  .JRUN*  Ml  BN*  KBbNT  ( 1  )*  .IRATE  ( JRLiM*  TFIN(JRljN)* 

1  LHANO*MRAND* KHAND* NoAND ( 1  ) 

2h  FhRMAT  (lrtl*3CX*7HRUN  M)*»  1 2*  10X*2/MN0.  OF  IMTlALLV  RAD  MODES** 

1  14«10**12HF INO  BPT10Na*I1/lH0«4X»12HREP0RT  HATE** 1 7* 10X* 

2  12RA IRISH  TIMF**T7«10X*1 3HRANDBM  SEE DS.*4( 19* IX) ) 

WRITE  ( 1  OH*  5 )  ( PAH AM ( I  * JR JR ) *  1*1* ]0) 

C  WllL  NO*  SEt  SECTION  NOS  AND  RESFT  OTHER  ELEMENTS  BF  NSrT 
l  \»0 

00  14  M|.  a  1 «  25*6 
ML7»Mi.*7 
00  14  1*1*25*6 
1 7» 1*7 
l  N«LN*1 


111 


351 

M2 

363 

354 

355 

356 

357 
350 

359 

360 

361 

362 

363 

364 

365 

366 

367 
366 
369 
37C 
371 
37? 

373 

374 

375 

376 

377 
376 
379 
360 
381 

362 

363 

384 

385 
366 
387 

386 
369 
39C 

391 

392 
39  J 

394 

395 

396 

397 
396 
399 
40C 


r.  ln  will  be  the  sect i an  \e 

Dfl  14  NaML/HL7 
D0  14  Ja I / 1 7 
06  13  Kat/? 

13  NSETI  N/.lz  K  )  a0 
NSE  T ( N/.J/  3 ) aL\ 

OH  14  Ma4/6 

14  NSET(N/J/M)aO 
OH  17  I  al / 16 
NbSECUjaO 

17  PCSECUJaO 

HF.'.O  ( 105*47)  I  DATA 
47  FARHAT  (1.1) 

IF  (ICATA  #£&*  1)  GA  T6  94 
C  tOATAal  JF  USING  ALTERNATE  FARM  AF  INPUT  BF  INITIALLY  BAD  NODES 
IF  (NEON  .EO.  0)  GA  T0  44 
WRITE  ( 108*53 ) 

53  FhNKAT  ( IHO/ 19Xz 33HL I  ST  BF  INITIAL  BAD  NBOES  FBLLBwS) 

DC  43  Ial/NHDN/10 

READ  (105/45)  <4A0(.J)»  jal#?0) 

45  FORMAT  (?0(12)i 

WRITE  (106/54)  ( LAf)(  J )  /  Ja  1«  20  ) 

54  FHkhAT( 1HJ/I0(?H  ( / IP/ lH/z I ?/ PH)  )) 

Ob  43  Hat/10 

JabAO(?aH-l> 

KabAO(?aM) 

0  BAD  NbLE  IS  (J/“) 

IF  (J  iFdi  44)  GA  Tft  *4 

f.  J/44  HAWKS  END  BF  LIST  5F  INITIAL  HAD  NODES 
CALL  RACK  (J/K) 

43  CHMINUF 

44  Oil  77  Hal/ 3? 

Ou  77  Nal/3? 

.JFVNTaNSFHh/Nz?) 

Nf  H(jC  aNf  F  I ( Mz  N/ 4 ) 

.iSliba5a.iEvNTaNEHH0al 
CALL  CAP(H/N/ALAHD) 

call  oranlxrnuM) 

NStT(M/Nzl)a«lNT(AI  AMD* ALbG ( HMJM ) *0*5 ) 

77  CONTI NUF 

IF  ( NSF  T ( 1 / 1 / 1 ) >NSFT( 1 /2/ 1 ) )  76/76/79 
76  NSET(l/1/5)al 
NSET()/1/«t)a2 
NSfc  Kl/1/7)  a9599 
N5ET(l/1zM) a 9999 
NSfcT(  t/?/5)a7777 
NSfc  T(l/?/6)a7777 
NSET ( 1/ ?/ 7  )  a  1 
NSfc I(t/?/6)al 
r.h  TO  63 
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401 

402 

403 

404 
40b 

406 

407 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 
4  3C 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 


7 9  NSLT ( I*1»5)a7777 
NSU(l*1»6)«7/77 
NSLTI 1* 1  *  7 ) al 
NSLT(1«1*a)*2 
NSH(1*2»&)«1 
NSfcT (l«2*6)al 
NSt  T( 1 #2#  7 ) «9999 
\Sl  T  ( 1  *  2*  8 ) *9999 

83  ri.»o 

Oh  80  1«1»32 
1 1«  1 

Oh  M  j*l*3? 

JJ«j 
Ml.»hl>l 

IF  (ML  .Lt.  21  Ift  T8  80 
*TKibm»Nsrrm»jj»i) 

CALL  HRRF.*(  I  1#  J.l) 

80  C(tM  INUF 
KMlU-0 
RETURN 

99  hrtlTfc.  (108*101) 

101  Ff*RMAT(lHd/lHC*16X*33HlNlTIAL  GRID  (M'S  MARK  HAD  N8DFR)/1H0/1M0) 
DO  127  1al*32*2 

HEAD  (105*121)  (MAP( J) * J»1 *  64 ) 

121  FhfiMAT  (6411) 

r  TWO  KdwS  OF  tMfc  GRID  COMPRlaF  BNF  LINE  IF  DATA 
hUllt  (108*122)  (MAP(J)* J»1 » 32 ) 

122  FhKmAT  (in  *32(11*1X1) 

C  HAD  NhOES  •  8*  GBOO  NROCS  ■  1 

WHlTt ( 108*122)  (MAP(J)*J«33*64) 

Dh  123  Hal  *  32 

IF  (MAP(M)  .Ew»  1)  8*  T#  123 

jal 

Kah 

CALL  HACK  (J*K) 

123  COMIMJF 

Ob  127  Ha33* 64 

IF  (HAP(H)  •  £$•  1)  08  T9  127 

.!•  I  ♦! 

MM-32 

Call  rack  ( j*<) 

127  CHNT1NUF 
fid  TO  44 

C  J81NS  PROGRAM  RHfcRF  REGULAR  FORM  BF  INPUT  F.NDED 
END 
C 
C 

SUBROUTINE  CNVRT(N» IR9W*IC8L) 

Oft  10  I  a  1 »  32 
IF  ( N«32 )  5*5*10 


113 


451 

452 

453 

454 

455 

456 

457 
456 
459 
4*0 
4*1 
462 
4*3 
4*4 
4*5 
4*6 
4*7 
46* 
4*9 

470 

471 

472 

473 

474 
471 
47* 

477 

478 

479 
4HC 
481 
4H? 

483 

484 

485 
48* 
487 

486 

489 

490 

491 

492 

493 

494 
49b 
4*6 
497 
496 
499 
*»OC 


lii  Nal>»32 
5  INtrfaN 
lf.tL»  I 
RFTu*N 
FNlJ 

r. 

C 

RljbhbUT  INC  FVNTR  (t.J) 

r.MMMBN  N.HET ( 32*  3?«8 1 *  PARAM C  1C*  10 I *  JRAU  ( 10 I *  TFINC10I*  ATRIB(H)* 
1  KhlTt*  JRUN*  KAUNT (10)#  ISFED*  TNUW*  NRDN*  NbSEC(lM*  PCSECIl*)* 
?  nRAND ( 1  ?M )  j  I  RAND*  HRA«D*  KHANO*  nr, 090*  KHAD*  PARAP2I  III*  101 

CPMbfc  CiCuDEI 5000) *  JAR|i**32»*  INTI ( 2 )*  INT?<?) 

DIPtNSIPN  HAPI3PI 
INTEGER  ATR lb*  TN0W*  TUN 
INTEGER  GCODt 
I  Ml 
•JJa  J 

JEVNT«ATRIBI2) 
imjEwnt 

IF  (JtVNT  .EL.  0)  ATK  ld( 2 )  ■ 1 
IF  (JEVNT  .FCl.  II  ATWI9I2M0 
NRLT(ll»JJ*?MATRlbC?) 

KFVNT  a ATR I B( 2  I 
NdbRaa  Ml 
NMfcCLaJ 

TALL  H Ml  <Nh-jR<.*NHhCL*  I2l 
NHbRWaMl 

CALL  FIND  ( MluRdj NMHCL* 12 ) 

NMBkwat 
NMfcCLaJ*! 

Call  find  (NKoRw.NM.5CL*  1 2 1 

NHcCLaJ«l 

CaLl  FIND  (NMrsHw*NMHCL*I2> 

DM  1  Llal«M 

1  ATKIbCLI MNSET(|J*JJ*Ll> 

NFhkDaATHItf(4Ml 
NRECTaATRIBO) 

JSLbanajFVNMNEMBD 
KSLnab**E  VNMNEmHU 
CALL  CAP! 1 1 *JJ« ALAND) 

C  ALA-lL  lb  THE  MEAN  INTFRARRIVAl  TIME  I.F*  1/LaMRDA  0*  1/*ll 
CALL  ckanDirnlm) 

ATK  lb  1 1  MTNttN>lNT(Al  AMDaALbf,  C  KNUM I -0.S  I 
C  JFVNTbO*  WEVNMi  IF  NrlDE  CHAN6EL  FRtIM  UrfcttO  T9  BAD*  AND  VlCF  VfRSA 
ATr>lb(3)aNSECT 
AT*»lb«4laNEHBJ-1 
if  «*>fcv?.r  tFo.  n  or  Td  10 
CALL  .JURE?! •PICTIJKF 1  * JAR| II* JJI* INT2) 

C  UPLATtb  INTENSITY  HF  P01NT(I*J)  IN  DISPLAY 
C  JaR ( I  * J I  ■REL  AD0RFSS  0F  INTENSITY  INSTR.  FAR  PHINT  C!»J) 
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1 


I 


1 


Nt»DKaN3f>N*l 

,NrtON  16  M  BF  BAD  NBDFS  IN  SR  ID 
I  WiStC(NfiFCT)sNhSFCCNRECT)*1  , 

NRiEC  IS  Mi  BF  bAO  NBOFS  IN  SECERN 
nn  Tb  i? 

V,  N-iUNaNBONM 

nmseccnrfct i»nbsfc«nsect)»i 

call  ddrepi »picturf • # jar« 1 1« jji« inti ; 

2  MSEC»NSSEC«NStCTI 

PC&EC ( NSFCT ) ■HSFC/AA*U 
OB  99  If  1  r » 

9  NSIT( II»JJ»IlaATRIBII 

IF  IKRITE  cEO*  01  r.B  TB  S7 
BStCaNBDN 
PCTM»»NSFC/10if*.0 
hRITEl 1PA*S?)  TNBW#  NhON*  RCTNB 
2  FtiAHAT  I t HO/lNO#BX#SMT INE»»  I5*SX» 17MNB*  BF  BAD  NBOCRa, 
1  U.S*»9HPCT#  RADa«F*»«l 
0(1  IB  MRel«32 

OB  li  MCal,32  \ 

IF  (NSET(NR«MC«Pl)  14«IS«1* 

«  MAP(MCla«f)AD 
SB  TB  It 

$  HaPINCURGBBD 

I  CBMINMC 

MUTE  (10**21 1  (MAP (HA I »HMa| » 3P l 
•  CBNTlNliF 

X  FBRMAT  (1H  «3?( At* 1  Ml t 
1 7  CALL  BRnERI  lltj.ll 
9  NEToBn  1 

FNO 


bubrbutine  capj i * j* alahdi  , 

CBMMbN  NRET ( 32# 1?* Alt  PARANI 10. tOMJRATE I 10 ) *  TFIM10I*  ATRIBI i)» 

1  KRITE*  JRUN*  KBUNTItuI*  ISFED*  TNBN*  NBDN,  N3SECI1AI*  PCSECdBI*  \ 

2  NRANOI 12*)«  I  RAND#  MRANO*  KRANO#  KGBBD*  KBAO*  PARAM2I 10* tOI 
INTEGER  ATRIB*  TNBW»  TF IN 

NCBbEaSaNSCT ( I*. 1*2 laNSETI I* J«*l«t  ! 

IF  ( l*EQ*lA*AND«J*Ffl*16)  GB  TB  30 
IF  ( I •  EQ«15*AND*J*Ff3«l6)  GB  TB  30 
IF  ( I*EQ*lA*AN0«J«F3*ltl  GB  TB  30 
ALAMO*PARAMP(NCBOE*.JRUN| 

RETURN 
30  CKNTINUF 

AL  AMD  aP  AR  AN ( NCBDF  » JRUN I 

RETURN  I 

END 

SYMBBL  , 

CtF  ISITBN 


I 


I 


\ 


551 

JSPflN  C*t2#l 

C 

S52 

kb*0 

c 

M3 

sTcr 

i 

55* 

1 

MS 

SCS#9 

• 

j»S* 

\  L**.4 

13 

557 

Lh*n 

1#* 

MM 

ANU*3 

•13 

559 

a 

It* 

StC 

fc*D 

1 


I. 


\ 


\ 


\ 
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Random  Grac 


MAIN  FHMiRAM 

rt*KMbN/7/«SETt3?.J?»Iif ) »F *N AM (1 O* 1 U ) *  jSATtl  10 1  * TF I N ( 1 0 1  *  ATHIBlBl* 
l  aMTI#  JHUN*  KmiNT(IO)*  ISFEO*  TNow#  NBDN*  NBSECdM*  PCSEC(16>» 

9  nRaNO ( 1 2A )  «  I  RANO*  PRA»D*  KPANJ 

OIMNSIPN  NR6*  (  A I «  NC9L  ( 8 ) 

INTEGER  ATRIb*  TN9W,  TF  IN 

RF A0(1G5*  1A3I  LWAi»0*MRAND*<CPAN0»  NR  AND!  1  > 

1  A3  FHNMATUfSI  . 

Or.  1 A  A  .1*2*  1KB 
1 A  A  NWANl)<J1*NRAMX.I-l  W?*J 
CALL  GRAPH 
NF  AO  ( 103*21  I  IliRAF 

21  F  OhMAT  (Ml  A 

IF  ( IfaRAF»tO»£.I  Gh  T9  9b 
WRITE!  lf)M*  37 ) 

3/  Frt**MAT  ( 1H1»?(  Aj(«  AHN90E*  lb»»9HNF  I  GHbBRP*  23X  II 
ON  9A  1*1*3? 

Oft  9A  0*1*3?*? 

JF  «*J*1 

Oh  a  k*9*12 

MF*NStT(t«J»fc) 

CALL  LNVHKME*  jnn*»  JC.9L1 
NUt-«(k«AI*JRO* 

NC.t. L I  k»M  1  ■  JCttL 
Mf  *NSET( I«JFX*<I 
CALL  LNVRT  (  ME*  JRH*«  .)C.hL  1 
NWbw(»*AI*JR6«i 
A  NCbLI N"A I ■ JCBL 

wWlTEdOA»171  I*.i«SiRf»*.<n*MC0L(1  I  «  NNPk ( 2 ) *  NCPL (  ?  I  «NRBM  3 )  t  \CBL  (  3) * 
1  NN6*( A I  •  NCUL !  A 1  # 

?  I.M*NR«w(5»*NC0L(f»)*NR(»  ,  ( A  1»NCBL  (  6  1 .  NRR*  (  7 1  #  MCBL  (  7  1  * 

3  NHbMA)*NCBL(M| 

17  FhRMAtllH  *2(9(3H  <  ,  I?»  lH»  *  I?»lHn*  10*  >  1 

9*  (ONTIMF 
9b  ChNTlNUF 
JMLN*1 

b  CALL  l)ATA'» 

CALL  GASP 

IF  (jRLiN  »UE*  If. I  M  T9  97  „  _ 

l*st  nun  js  np  io,  this  Card  can  be  alieRfD  t»  alls-/  anv  nu^ee* 
tF  runs  up  tp  io  rith  oifff_rfni  paaametfRr  fbh  laCh  Ri.n. 

JhUN*uRU4»l 

r.t*  TO  s 
b7  RTcl 
F\L 


loi  an  T9  u7  .  „  „ 

o.  this  Caro  can  be  alfeRfD  t»  allo-i  anv  nu^e* 

Trt  10  R|TH  UlFFcRFNI  PaAaMETFRR  FbH  LACh  Run. 


Sl.fckblT  I  ti.  OR  AND  IRMJM) 

Ci*N  MCN///NSET  (3?«'1?,1?)  »PaFi  AM  (10<10 )  *  jRATE  (10)*TFtN(1C)»ATRIBC8I» 
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51 

52 

53 
5« 
55 
51 
57 
5k 
59 
to 

41 

42 

43 

44 

45 
44 
47 
4« 
49 
7C 

71 

72 

73 

74 

75 
74 
77 
71 
79 

ac 

ai 

12 

S3 

14 

as 

14 

17 

11 

19 

90 

91 

92 

93 

94 

95 
94 
97 
91 
99 

100 


1  UNITE*  JRUN#  MOUNT (10)#  ISCED*  TN8W*  NBDN,  NHSECI16)*  Pr.iiECIll), 

2  NRANCI 12H)»  I. RAND*  MRAND*  MRANO 
L  RAND«LRAND#45539 
MRANbaMHAND*335544.33 
JalalAHSILR AND 1/147772 14 

RNl/H*  .SaFLBATINRANOlJ)4LMANDaMRANO)a  .232A3044E-9 

KNANDaKNANDa3624340A9 

NRANDlJIaMRANO 

RETURN 

FNO 


C 

C 


c 


c 


c 


c 


c 


r. 


c 


SUBROUTINE  ORDER  IM.N)  , 

COhHON///NSET  I3P»3?»12)»PARAMI  10.'  101* JRATEllO)*TFINI10)» ATRIBIB)* 
1  aRITE*  JHUN*  MOUNT  1 10 ) «  I SEED*  TnOW*  NBDN*  NOSEC  1 14) •  PCSEClll) 
INTEGER  ATS*  AT4 «  AT7 «  ATI 
INTEGER  ATRIb.  TNOw*  TFIN 

OROFR  ADDS  NODE  IM.N)  TO  8RUERE0  LIST  OF  EVENT  TINF.S  FOR  DAT  AN 


MMO 

M!»ManS  M7  ARE  FLAGS.  WHEN  8»TH«I  WE'VE  FOUND  RIGHT  SPOT  FOR  (M»N) 
1  al 


Jal 

KNRri  THAT  NOOE  11.1)  HAS  BEEN  OROERED 
IF  (MRITE  .EG*  R>  SR  TO  14 
I aATRIB(S) 

J.ATRIBI4) 

IF  U  »NE.  7777)  GO  TO  14 
I aATRlBI 7) 

JaATRlBIR) 

14  IF  { ATHIBI 1 )«NSFT ( I *  J« 1) )  5*4*7 

.LT.  0  MEANS  GO  TO  PREDECESSOR,  ,GT.  0  GO  TO  SUCCESSOR 
5  I AaNSETl I  * J«7 ) 

IF  |M7  .EG*  1)  GO  TO  9 
IF  (IA  .Ed*  9999)  Go  TO  9 
SFF  IF  IT  HAS  A  PRFDFCFSSBR 
JaNSET | I«J«S) 
tall 


Mh«l 

GO  TO  14 

7  I AaNSETl !«J*5> 

IF  IMS  *EU*  1)  SO  TO  4 
IF  ( IA  .Ed*  7777)  SO  TO  4 
SFF  IF  IT  HAS  A  SUCCFSROR 
JaNSET I I.J*4> 

1  a  I A 
M7«l 

GO  TO  14 

4  ATRIB 1 5) aNSET 1 1 • J«S)  . 

|M,N)  SUCCEOES  ll*J)  IBV  CONVENTION  IF  THEY  HAVE  •  TIMES) 
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101 
102 
tOi 
104 
1 0h 
106 
107 
10* 
109 
UC 
111 
112 
113 
1U 
11b 
116 

117 

118 
119 
12C 
121 
122 
123 
12* 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 
146 
149 
15C 


ATK  Id  ( f>)  a.»SET  ( I #J# A  ) 

ATMb(/)a| 

ATK | U ( * ) ■ J 
NsETIltJ#*)"* 

NSETI  l#.l»6)aK 
ATbaATRTH(b) 

11  (Alt.  .tU.  77771  .18  TO  9* 

f  TEST  MILS  II  (t#.l)  MAS  A  SJCCESSAR,  AND  THUS  MUST  UPDATE  HIM 
ATfcaAtR 1b(6) 

NSL  t ( Alb# AT6# 7 1»M 
NSf T ( ATK#  ATA#  m ) aN 
r,n  Tb  9A 
9  ATK  lb (b  I  ■ I 
r.  (M#M  PRECEDES  <1  < J) 

ATK |B ( 6 1 ■ J 
ATKlb(7)a*SET(l#J#7) 

ATK lb(tt ) sNSET ( 1 « J«8 1 
NSETI I #J*  7) »M 
NK4  T(  I#.l#«)aN 
AT7aATRtB(7) 

IE  (AT7  .tQ»  99491  08  TB  9h 

C  TEST  KAILS  IF  llt.il  HAS  A  PREDECESSOR  (-HICH  MUST  BE  UPDATFDl 
ATMaATRIRIB) 

NSETI AT7#AT8#b>*M 
NSET  ( AT7# ATR«6)aN 
96  Dti  99  KabtR 
99  NSET(M#N»<leATRI8(K) 

RETURN 

end 

c 

c 

SUtiRbUT  TNL  EVNTS  (la.il 

r.BMr,6N/7/NSET(3Pt3P#l2)#PARAMM0#10l«  J8ATE(10I#TEIN(1o)«ATRIB(8I# 

1  aHTE#  JRUN#  KBUNT(IO)#  I  SEED#  TN8M#  NBON#  NBSECl 161 #  PrsECIlBI# 

2  NNANOI 128) *  1  RAND#  HRANO#  KRAN'J 

DUENSIAN  I  STM  (1024)#  1 01  ST  ( 10241#  IMAXI1024)#  ISTKI1024) 

INTEGER  ATRIB#  TNBW#  TF IN 

I  1*1 

JJ»«# 

J(VnT*ATRIB(2I 
I /■ jEvnT 

CALL  ALTER!  1 1  •  J.i#  1 7  I 
D»  1  Llal«8 

1  ATS  IB (Ll  )  aKSET (IT#  JJtLl I 
NEHdL*ATR|B(4U1 
NSECTaATRIBO) 

JSLH'ba.lE  VNT  AnEhBD 
IE  ( JE.VKT  .ECl.  01  ATR  18(2) a  1 
IE  (JEVNT  •F.O#  1)  ATRlS(2)aO 
KFVNTaATNlB(2) 


r 


1S1 
1  52 
153 
15'. 
155 
1  !>fc 
1117 
I5t< 

159 

160 
161 
162 

163 

164 

165 

166 
167 
I6e 
169 
1  70 

171 

172 

173 

174 
1  7b 

176 

177 
1  76 

179 

1 80 
1M 
162 
1  M3 

164 

165 

166 
167 
16b 
169 
19C 

191 

192 

193 

194 

195 

196 

197 

198 

199 
20C 


KSLh«5»<F  VNT4NEH6C 
Al  AHUaHAHAM(KSUH» JRUN) 

C.  ALAMO  IS.  THE  Mt AN  TNTFRARRtVAl  TIME  I.E.  1 /LAMBDA  Ok  1/MU 
CALL  OHANU(HNUH) 

ATkTb(  D.rNeb- INK  AL AMCaALBG ( &NUM  1-0.5) 
c  JF.VNTbO,  <K VNT»1  IF  NflOt  CHANGED  FROM  (ilOD  T»  BAD*  AND  VICE  VERSA 
ATMbOtaNSELT 
A1kIb(4)»NCHh')»1 
KNliNT ( KSUrt ) aKliUNT  <  <SUb ) 41 
KtlUNT  I  JS.UU )  abSUNT  ( JSUH  >  «1 
TF  ( JEVNT  .EL.  It  r.6  TO  10 

NbLN4(.6RN4l 

C.  AMOK  IS  N6  OF  dAD  NODES  IN  GRID 
NrlSEC  ( NFFCT )  aivBSFC  (  NSECT  )  *1 
r  NiARF  C  TS  NH  OF  MAO  N6DFS  IN  SECT 
bn  TO  12 
1C  NhLN»>i9RNa1 

Nr-S)  C  ( NSFCT ) aNbSFC ( NSECT 1»1 
It  RSECaNBEFC INSECT  I 

PCSEC ( NFFCT )akSFC/A4.G 
05  99  I al « 2 

99  NSE  T( I T< JJ* I )  .ATRIH ( I  ) 

TF  (<H T TE  .EL.  0)  36  TO  a7 

IF  (  TNOk.LT <1600 )  06  TO  3()5 
6  LVL1*C 

D»  216  Mai,  1024 
21a  ISTMiritaa 

DIj  IK  |Hal,32 
Otj  11C  TCal,3? 

TF  (Nat T(TK,1C«2).F!J.U>  60  TO  110 
Nxa  |fc.32a( |C*»1 ) 
l VL 1 "L VI  141 
TsTkCLVI 1 taNK 

IK  CONTINUE 
MP»lVl141 
05  111  Jl "HP* 1024 
111  TSTbKDaO 
SuMaLvLI 
Nl  bTSaO 
NOL'ESaO 

TF  (LvLI  )  33*33,30 
29  TF  (LVLI)  32,32,30 
3C  CALL  CNVRTI 1STKCLVI  1 ),NR*NC) 

TSTkILVLI t»0 
l VL1»LVI 1«1 
LHTa  1 
LyL2alii24 
C.fl  TO  67 

114  TF  (LVL2.E3.K24)  06  TO  lb 
LVL2»LVL24l 
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201 

CALL  CNvWTl  16TMI  VI  ?).M»NCl 

702 

IF  (NSl T (NO.NC* a  1 .1  T«*l  JKNT a JKNTal 

703 

ISTkILVI ?)*0 

9 0 » 

LPl-LOT-l 

70b 

*7 

IF  CNfc-.VI  21*2*2 

706 

21 

MMlaM.1  /v. 

? 07 

NC1*KC 

LINfl 

NVaNhla3?.(NCl-H 

IF  ILVL1.ME.0I  fill  T«  117 

IF  (LVL2.E0* 102*1  -3  TB  16 

LhT*Lt)T.l02*-LVI  2 

20A 

709 

77 

?1C 

711 

212 

713 

l  VL J*LV|.  2*1  ^60 

7i* 

ON  302  I6*LVL3*102* 

2  IS 

CALL  CKvHlf  !i>Tkt!S).MH«»C) 

91* 

IF  (NbET(MR.KC«*)«LT.*l  JKNT  aJKNT.l 

717 

3o? 

CONTINUE 

719 

SB  TB  If 

919 

117 

Db  20  INal.LVLl 

920 

IF  (ISU(lh)-WI  20.119*20 

92 1 

119 

|STMLVl2t*NV 

722 

L vL2«LVI  2-1 

723 

r.  LVL7  IS  AM  OPEN  SPOT 

22* 

ISTK( IH| alSTKlLvLl 1 

72b 

1 RT  AC  L  VI.  11*0 

92 1 

l  VLl«LVI.l*t 

921 

GN  TB  121 

222 

20 

CONTINUE 

229 

121 

GO  IB  I2.*«6«ll*|*l  TME 

230 

2 

IF  IMN-Ii  *»*»3 

231 

3 

NtilaNFt-1 

232 

MCiaNC 

233 

LlNfe»2 

23* 

G»  16  77 

23S 

* 

IF  (NC-321  5*6*6 

23* 

9 

NMlaNfc 

237 

MC1*MC*1 

232 

l  I  ML ■  3 

239 

Gb  TB  77 

2*0 

6 

IF  lMC-11  11**11*. 7 

2*1 

7 

NKlaNn 

2*2 

^laNC-1 

2*3 

l  IKE** 

2** 

GO  TB  77 

2*S 

c  havf  poller  out  all  bf  bad  mshk  bf  i*m*nC)  *no  placed  tmfh 

2*2 

c 

AT  bBTTNM  BF  STACK 

2*7 

GN  TB  11* 

2*2 

16 

CONTINUE 

2*9 

ISTIMLBT1*ISTN1LNTU1 

2SC 

Ml.bThaKI  NTSal 

121 


251 

P52 

32 

?53 

33 

P54 

PSS 

PS* 

PS7 

311 

PSS 

30* 

?59 

P*0 

P*l 

P*2 

?*3 

P*4 

301 

?*S 

P6* 

?*7 

PSS 

?*9 

P70 

P71 

P72 

303 

P73 

?74 

P75 

P76 

277 

304 

27S 

P79 

PSO 

?S1 

PS2 

PS3 

307 

PS4 

2*5 

305 

PS* 

306 

PS7 

310 

PSS 

309 

?*9 

312 

29C 

57 

291 

19 

P92 

293 

r 

P94 

r. 

P95 

?9* 

P97 

29* 

299 

300 

c 

Gh  TB  Z<* 

PLBTfcaNIBTS 
CONTINUE 

IF  (ISEEO.GTfC)  GO  TB  303 
OB  311  J»lf10?4 
IDIST ( J) aO 
IMAXI  jUO 
CBMInL'F 
iSElDatSEEOal 
On  3b 1  Jal«10«4 
IF  (  ISTMC J)  •LJ*0)  OB  TO 
lOISTlJla  ID  ISTI.il  ♦ISTNIJ1 
MAXaJ 
CBMINUF 

IMAX(MAX)a;MAX(MAX).»l 

IF  (ISEFn.LT.2001  GB  TB  JOS 
UNITE! 10*f 312)  TNfiU 
KRL'MaO 
KSUMaO 

Dt>  303  Jal#  1024 
xSUrta  aSIjM*  1 01  ST  (J) 

KNLMauROM*  I  MAX  ( .)> 

OB  304  Ja  1  *  1 0?4 

IF  C  IlilSTCJl.EL.O)  G9  TB  304 

STF  MaFLBAT ( ICISTl J) I /FLOAT (KSUM) 

UN  I TE ( 1  OB*  306 1  .if  STFM 

CHNTINUF 

UNITE  ( 10A#309) 

OB  307  Jalfl0<!4 
IF  (IMAXUlfEO.O)  GB  TB  3(»7 
STEMa  FLBAT(IMAX|J11/FL0AT<KRUM) 

UNITE  <1  OK# 31 Cl  J#  STEM 

r.nNTiNut 

HFTliKN 

CBNTlNUF 

FHNI1AT  (1H  fUMCIUMP  S  !2Ea#  1 4  jSX#  lOHFREOUENCYa#  F7.5 » 
FBMAT  11"  fllHCLUMP  SIZE.#  I4#5X#13HFHE0.  AS  MAXa#F7#51 
FBNMAT ( 1H01 

FORMAT ( lHJ# 5X#5hT IMFa# 14) 

CALL  bNRENC I IfJJi 

RETURN 

FNO 


SUBROUTINE  SACK  I  Jt K ) 

CHMnBN/7/NSET(3?#32#12>#PAMAM(10#1U1#.jRATEC101#TFINI101#ATRIBI81« 
1  XRITE#  JNUN*  KBUNT ( 1 0) f  ISEEOt  TNBU#  NBDN#  NBSECIU1#  PCSEC(1*1# 
?  NMANOIlZaif  I RAsDf  MRA*D#  KRAnJ 
INTEGER  ATRlBf  TNBUf  TF IN 
<Jf<)  IS  AN  INITIALLY  BAD  NBDE 


122 


.101 
102 
'0.1 
10* 
10b 
106 
107 
106 
109 
11C 
111 
.112 
113 
11* 
11b 
116 
.117 
11 8 
.119 
120 
121 
122 
123 
12* 

325 

326 
127 
.126 
329 
13C 
131 
112 
.131 
13* 
13b 
336 
.137 
136 
139 
.1*0 
1*1 
1*2 
3*3 
1** 
3*b 
1*6 
3*7 
1*6 
.1*9 
15C 


r. 


r. 

c 

c 

r: 


NSi  t  ( >Jt  K#  i  )  *  1 

kill  UPDATE  K«  hF  HAO  NjDES  IN  SFCTaM*  AND  PfT  HAD  IN  RFCThk 
NSfcLTaNRET  ( J»»»  1 ) 

K«S{  C (KSFCT  )  »KbSFC  ( NStCT  )  *1 
fcSIC'NbfiFCtKSkCT) 

PCStClNSFCT) aWSFC/6*.0 
OB  10  I «1j 1? 
t  !■  I 

nMNUtlJiKiII) 

IF  (Ht.HJfO)  RFTUSn 

call  ckvhtimejjrbw.jcili 
KHfcwa*»Rb«* 
nciiL'-ch. 

lu  NSt  T  I  KKC8*  NCtJt  #  *  1  »NSFT(  NRSW*  KC.BL*  *>♦! 

9  M  TUAN 

F  NO 


RUt-RbuT  INE  ALTER  ( J# < j \C8D(.. ) 

ALT)  N  UPDATES  NefcS  SF  tJtOt  WHICH  MAS  JUST  CHANGED  STATE 
THIS  COHRf SPBNDb  TH  THE  *8Ra  BF  FIND  IK  THF  NGN-RANDBM  6RAPh 

r.BK;SK/7/KSETl3P*3P,l?)#PAKAmiO#10)j  JSATEMOIj  TFIKMOlf  AHIBCB1# 
1  kK|TE,  JRUN#  WB4IKT ( 10 ) <  ISFEOj  TvBW.  KBDK*  NHSECI16I.  PC ;iEC{16) 
TrtTEGtR  ATRIb#  TN3 A<  TFIN#  ALT 
KK  tl  n  ■ 

KCfcL* A 

Ort  10  I »9» 1 P 
T  T- I 

MF  aNSt  T ( NHBh«\CH|  #|l) 

IF  (HE,EC.»0)  RETURN 
CALL  CKVHT I HE« jJR2h( JutoBL ) 

JRtl*a».JWrtA 

JCKL'-JCrtL 

AT  a  IBM  )  ■  VSET  ( JRBW«  JCSL#1 ) 

ATS Ib( ?) aKSE  T ( JRHaj JCSL»2 ) 

ATh  tb ( * ) aNSET ( JRHWj JC8L# *  I 
jTIMLaATHlbl 1 1 
JNEHBaATMIB(«)*1 

IF  ( NCOCF  •£&•  0)  ATRIB<*)«AT»fHI*)*1 

IF  ( KCttf.'F  *EU«  1)  ATR!6I*>aATRI8C«0«l 

KNthoaATR|S(*)Al 

jSUijAbA  ATh  1ft  1 2 )«  JKFHR 

KSUbab*  ATh  Ifi  I  ?  J  4KKFHB 

A BUNT  t  <SUn ) aKoUNT ( KSUB 1*1 

KBUNT « JEUrt ) aAHUKT ( JSUo ) -l 

ALAr.OaPAHAMIWjUHj  JRUKl 

CALL  URA.mO(RNvjH) 

Al.  T  aT\6a»  I  NT  I  aL  AMD#  Al  SGtH  WM ) -0.5 ) 

3  ATHblllaALT 
Of*  TB  6 


351 
35  c 
353 
351 
35b 
356 
157 
356 

359 

360 

361 

362 

363 
366 

365 

366 

367 
366 
369 
37C 

371 

372 

373 

374 

375 

376 

377 
376 
379 
SRC 

361 

362 

363 
384 
365 

386 

387 
386 
369 
39C 
391 
39c 

393 

394 

395 

396 

397 
396 
399 
«0C 


4 

5 

7 

6 


K 

8 


C 

C 


lu 


15 

25 


IFI ATKlft(lI-ALT)  5#6c7 
If  (NCftCFcNEcATHIHIP) )  ATn 16 1 1 > »AlT 
QH  Tb  6 

if  (KCPrK»fco#*THib(?) i  atsiriii.alT 
CftETINUf 

NSLT(wHPHiJCbLil)»ATRlB<l  ) 

NStT  ( uNSH* JCBL*  4 ) • A  TRI B ( 4 ) 

IF  (jMkK»KU»aThIB(  1 ) )  Sd  TB 
CALL  HMf* Vc  ( JhdM*  JCflL  ) 

CALL  BfiOFM  jKoWcJCBL) 

CtiKTIMJf 
RETURN 
f  \C 


SubMlLTlNE  KKBVF  JCBM 

f.HMMbN/7 /i»SE  T ( 3?  «  3?  c 1 2 1  *PaRAH  ( ’  3 1 1 0 )  •  JK  ATE  tlO)*TFINMO>#AT*<IBC8># 
l  kKITE#  JxUNj  KfliAT (10)c  ISFEC*  TNBiaa  NBDN;#  NBSECJ16)*  PC:»EC(16> 
INTEGER  AIR  lb#  TNfid*  TF  IN 
OH  10  6«ti8 

ATKIbtO*  NSfcT  ( jRuWc  JC6L«  a  ) 

I RtCH*ATRIb(5) 

ISUCOATRIBC6I 

|PRtH»ATRIB(7) 

|PKEC»ATRIB(B> 

If  (IPHM  *£0*  99991  OB  Ts  15 
NSEll  |PKf>j  |Phtr»b)«lsOCR 
Nfif.TMPKF-ljlPHEr#6l-ISUCC 
IF  (IRUCH  «F«.  77771  GO  19  25 
KSI  T(  ISlifwc  ISoCr#7l«IPR£R 
NSl  T  ( ISl,r«#  ISuCr.«8l«|PREC 
RfTohN 
f  NO 


c 


4 

r 


t: 


rtirKDN/7/i<SFT  I3P#  jp«12  IaPARAN C I o#  10  > c.iKATE  C 10  >»  TF  IN  ( 1  r,  >  #  AT*  IB(  8  >  # 
l  aKITEj  jhliNj  MMP.T(lt))#  ISFEO,  TNB4*  NHDN#  NHSEC(U)»  Pr$»EC(l6> 
jNTtGLR  aTRjB*  TN9W*  TF |Nf  RUtR#  jUCC 
MaTA-1 
M aTC«1 
J*  ITL»0 

IF  |NSI T ( NEXT**  NFaT  C»  7 )  »KJ»  99991 
TF31  AdoVf  FAlLa  IF  I \FXTR#.\fcXTC )  IS 


fVi 

N8T 


T9  S 

THE  hF  AO 


BF  THf  Lts>T 


»  »  NaL  T  («EXT>i»\FXTC»/l 
M  XTC«N6ETIStxTH«f.F<TC*8) 

\F*Th»\f  X 

SFf  If  HJ-  PHFOKCFSiiHR  IS  He  AO  BF  LIST 
G‘<  T6  4 


124 


401 

*02 

403 

404 

405 

406 

407 
406 

409 

410 

411 

412 

413 

414 
4 1  £> 

416 

417 
416 
419 
4?C 

421 

422 

423 

424 

425 

426 

427 
426 

429 

430 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 


5  CALL  *-MfcVL<N£XTR»NF  XTC) 

ITLSTatStTOMA) 

IF!  1UST.  \t*fc»  liH  Til  f.»i 
597  I TtST a  If  1  TOM  »> 

1F| ITfcST.tO.O)  Hfc  TUBS 
f.h  TO  597 
556  r.  iKTl  MUl 

TNUwaATKttil  1  I 
It-  <TN6«  auF  •  TF  IM.IRUN)  I  SO  TH  17 
SUCRaATFlMS) 

StjCt»ATAJ-M6l 

J11Tt»jKlTt4l 

XRlTEaO 

IF  URITF.  »LT.  JRATF I  jHUNI  >  OH  TB  7 
r.  IF  AcliVi  TFST  la  MfT*  WILL  NOT  PRINT  HfSULTS  AT  THIS  EVFNT  TIME 
JWITtaO 
KHITtal 

t  hill  PRINT  RESULTS  AT  THIS  KVf NT  TIRE 
7  CALL  EVNTSINExTH*NF xTC I 
NFXlKaSlJCH 
MXTCaSlCC 
Or  TB  4 
17  KWITtal 

RKITE  I10H.P1) 

21  FORMAT  (1hO#5?X.P6m»#  FINAL  RF PilRT  FOLLOWS  •• ) 

CALL  fctfNTS  (NtATH.NExTC) 

RETURN 

FNL 


C 

C 

SUBROUTINE  DATAN 

rBPMtlN/7/'*SETi3P,3P,12>«PARAMno*lC>»jBATUiCiI*THNncIaATRlBI8)« 
t  KRtTE,  JRUN*  KAI1NT ( 1U>#  ISEEU.  TNBW*  NRDN*  NBSECI 16)*  PrsECIl6>* 
?  NRANM  12H)*  IRANI)*  MRAND*  XRANJ 
DIMENSION  MAPI A4 I «  BA0I2OI 
iNTElith  ATRIb.  TNBw*  TUN*  PAD 
RFAL  LAMHOA*  MU 
IF  I  >LN  •  NE  •  II  Qfl  TB  ?9 
READ  <105*311  KHITF 

31  FHRMAT  111) 

IF  IKR1TC)  32*32*33 

32  05  19  J.l*10 

READ  1105*21)  I  AMPOA*  MU*  N 
RRITL  1108*10)  LAMBDA*  MU*  N 
IF  <LAMHOA»MU>  A* 7*A 
h  WRITE  1106*9) 

9  FORMAT  <50Hc£RHflR-  LAMBDA  MUST  BE  GREATER  THAN  BR  EQUAL  To  MU) 
JHLNafcO 
RETURN 
7  PaN 


125 


r 

r 


4b  1 

*52 

453 

8 

*54 

11 

45b 

456 

457 

456 

459 

ly 

460 

it 

461 

462 

21 

463 

464 

33 

465 

34 

466 

3b 

467 

Jfc 

468 

22 

469 

470 

23 

471 

472 

473 

4 

474 

3 

475 

476 

5 

*77 

29 

478 

*79 

2b 

480 

C 

481 

482 

483 

3H 

484 

485 

486 

26 

487 

488 

489 

490 

491 

3 

492 

493 

494 

495 

C 

496 

497 

498 

499 

500 

PARAM(1*J)*P/NU 
GH  TB  11 

PAR  AN  ( 1  #  J )  ■  1  *G/(I.AMBOAaMU) 

PARAM  ( 6*o )  ■  1  •  C/Mil 
DFuTA*M(j/b*0 
ca  19  i*2*5 

MIj*ML»DFLTA 

PARAMJ I,J)*l.L/(  LAMBDA -NO) 

PAK AM ( I+5*J)*1*0/MU 
FORMAT  MfiOLAMBDA*»F7.4*b»*3HMU**F  7.4*5X, 

1  22HRT7E  8F  QUEUE  IMG  Rc»BM**T5) 

FORMAT  (2(F7*L)«tb) 

GH  te  36 
Oil  34  1*1*10 

RF AC  (105*35)  (PARAM(J*I)*J*1*10) 

FORMAT  ( (Of 7*0 ) 

RF AD  ( 105*2?)  (JRATF( l  )*I*1*10) 

FORMAT  ( 1017) 

RFAO  (105*23)  (TFIN(T)*I*1*10) 

FHRHAT  (1017) 

Ofl  4  u* 1  *  10 

.RITE  ( 108#  3)  .)»  JRATE(  J)j  TFIN(J) 

mRITE  (108*5)  ( PAR  AM  (I*J)*T*1*10) 

FhRMAT  (1H0»7HRuN  Nfl.,  I2*10X*1?HREP8RT  RATE** 17# 10X* 

)  12HFIMSH  TIMF./I7J 

FHRMAT  (Irt  *  11HPA’<AMETERS**10(F7*2*5X)  ) 

RFAO  (105*25)  NtfiVl 

KRlTEaR 

FHRMAT  ( I 4 ) 

NI8N  IS  N6  JF  IMTtAI.LV  6AD  NHDF  S 
NnLMMRN 

RFAO  (105*38)  ISFEn 
FORMAT!  M) 

USITt  ( 108*26)  .(RUN,  MlBN*  tSEFO*  jWATE(jRUN)*  TF IN (jRLN  1 , 

1  LRAKn»MRAN0»KRARD»NRANO( 1 i 

FORMAT  ( 1H1*22X.13HRAN08M  GRAPH  * 

1  7HRUN  N8. »  I ?*  1 0X*27HNH*  HF  INITIALLY  RAD  N.9CES** 

2  14* 10X* 12HF INI  HPTT&N-# i 1/1M0»4X* IPhREPPRT  RaTE*« |7* 10X* 

3  ) ?HF IRISH  TINE**  t  7* lOx* 13MRAND8M  SEEDS** 4 ( 19* 1 X ) ) 

RRlTk  ( 108*30 )  ( PAR AN ( I* JRUN) « I *1* lb) 

FORMAT  ( 1H0*9HPARAH(1)*»F7*P*10X*9HPARAM(2)*,F7*2* iOX«9HPARAM(3)*, 
1  F  7* ci 10X*9hPARAM( 4  )  *»F7t2* 10X*9HPARAMI5) **F7*2/1hO*  9HPARAM( A ) ■* 

?  F  7*e* 1  OX* “HP ARAM ( 7 ) **F7*2* 1CX*  9HPAHAM ( 8 ) **F7*  2* 10X* 9HRAR AM( 9) ■* 

3  F7*2*10X» 10hPA«AM(1O)*»F7*2) 

RILL  N8r  SET  SECTION  NOS  ANU  FESfT  BTHER  ELEMENTS  BF  NSFT 
LN*0 

OH  14  Ml*l«?5«8 
ML7*ML*7 
DM  14  I .1*25*8 
I7*I*7 
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<501 

502 

*03 

504 

5Ca 

*06 

50/ 

soe 

50V 

510 

511 

512 

513 
515 

515 

516 

51 7 
516 
519 
52C 

521 

522 

523 
5?4 

525 

526 

527 
526 
529 
53C 

531 

532 

533 
536 

535 

536 

537 
536 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 
55C 


r. 


r. 


r. 


c 


i  .■ln*i 

lim  aTll  fir  THt  section  no 
Oil  14  NsrtLaML/ 

04  14  Ja  I  a  I  7 
Oh  13  Ka1*2 
1J  NSE  T  ( Na.i*  <  )  *G 
NSE1(N«.i*3)aL\ 

Oh  14  I“la4*6 
la  NS)T(N».l».1)aO 
Oh  37  la  la  10 
37  KULhl < I ) aO 
Of*  17  |  a  1  a  1 6 
NhbtC {  I  )  aj 
1/  PrbtL(I)aO 

fif  AO  <105*47)  I  OAT  A 
47  FHhMAT  (111 

IF  (IjATA  .EG.  1)  G«1  Tb  94 
lOATAal  IF  USING  Al  TFRNATE  FORM  HF  INPUT  6F  INITIALLY  RAC  NftOfS 
IF  (NcDh  »LQ.  0)  On  TO  44 
WHITE  (109*53) 

53  FhhMAT  I l«0al9K*3JHLIST  SI  INITIAL  BAG  NODES  FOLLOWS) 

DR  93  Ial«KBDNa10 

RFAj  <105*45)  (RAD ( J ) *  J* 1*  PO) 

4b  FhhMAT  (PCI  12)1 

WHITE  ( 108*54  )  ( *5  A  i  J  ( J  )a  Jal  *20) 

54  FffkMAT ( ImC* 10(2M  (  . IP* 1M, , 1 ?, ?h )  )) 

DU  43  M. la  10 

JabAD ( ?aM* 1 ) 

KabAO<?*M) 

BAD  NODE  IS  (J*K) 

IF  (J  •((!•  44)  06  TO  44 
j.44  tfARrtfi  ENC  SF  LIST  OF  INITIAL  BAD  NODES 
CALL  HACK  (J**) 

45  CSMlMJF 
44  Dh  77  Ha  1*3? 

DH  77  Ka 1* 3? 

.if  VNTaNSFT  (M*N*P) 

NF>£iD*NSF  T I Ma<a«  4) 

JSLb*b»JFYNT4NEHrtD*l 
Al AMUbPaHAHI JSUft.JRUN) 

TALL  DRAnO(RNUF') 

NSEKKaK*  1  )  a- I  NT  (  AL  AMDaALOu  (RNUM  )-0.5  ) 

77  KBliNT  ( JS03 )  aAOUNT  ( JSlib )  *1 

WnITE  (106*64)  ( KOUNT ( 1 ) a  I ■ 1 « 10) 

6a  FnhnAT  (lM0*7HKNT(llaaI4a2Xa7HKNT{2)a«I4a2Xa7HKNT(3)a*I4*2Xa 
1  7HKNT(4)Bal4aPX*7M<NT(5)a,I4a?Xa7HKNT(6)aaI4aPXa7HKNT(7)aal4a 
?  2 A*  7HKNT  ( 6  J  a*  1 4a  P  A*  7HKNT  ( 9 )  aa  1 4a  2X*  8HKNT  ( 10 )  a*  1 4*  ?A ) 

IF  (NSET (1*1*1) »NSFT (1*2*1))  76*78*79 
76  NSET I l*1«5)al 
NSET( 1* 1*6) *2 
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j)»ilia4aNMi<NNfBAiM>66ta^lB4iaaaMateaaM«4aa6lfc4ai6a6AaaMi<4Mi*4fa 


551 

552 

553 

550 

555 

556 

557 

551 

559 
S6C 

561 

562 

563 

560 

565 

566 

567 

568 

569 
S7C 

571 

572 

573 

570 

575 

576 

577 

578 

579 

580 

581 

582 

583 
580 

585 

586 

587 

588 

589 

590 

591 

592 

593 
590 

595 

596 

597 

598 

599 
60C 


NSETl i« 1*  7)a9999 
NSE  T  ( 1  *  1  •  H  ) *9999 
NSET ( 1*2*5) a7777 
NSE T(l#?#6)»7777 
NSE  T  ( 1  *  P*  7 )  a  1 
NSETl 1*2*81 al 
6a  TB  83 

79  NSETl l*1»3)a7777 
NSE  T(l*1*6)a7777 
NSET(1«1* /)«1 
NSETU*1»H|a2 
NSET(l«2*B)al 
NSET( 1 «  2*  6 ) al 
NSETl 1*2*  7! a9S99 
NSET ( 1 «  2*  8 ) a9999 

83  ML*0 

Off  fcl)  la  1 «  32 
11*1 

OB  80  Ja 1 #3? 

JJaJ 
ML»MLoI 
IF  (ML  »LE»  2)  GB  TB  80 
AT*Ib(l)aNSET(!f»JJ«l> 

CALL  BROCK ( 1 1  a JJ1 

80  CONTINUE 
KRITEafl 
RETURN 

99  WRITE  (103*1011 

101  FflfiMAT ( 1 HO/ INC* 16X«  33HIN1 TI AL  GRiO  (M'S  MARK  HAD  NBDFS 1/lMO/lHO! 
OB  127  |al*32*2 

READ  (105*1211  (MAPI Jl* J«1 *  AO) 

121  FflRMAT  (60|11 

C  TWfl  NOWS  PF  THE  GRID  CBmPRUE  BNC  lINE  6F  DATA 
WRITE  (108*1221  (MAPI Jl* Jal*32) 

122  FORMAT  (1H  *32(11*1X1) 

c  bad  mjdes  a  i,  ubbo  nodes  a  1 

WRITE (108* 122)  (MAPI J) * J*33«60) 

DH  123  Mai *32 

IF  (MAP(M)  aEOt  1)  GS  TB  123 

Ja  I 

KaK 

CALL  HACK  (J«K) 

123  CONTINUE 

Dm  127  Ha33*60 

IF  (MAP(M)  t£Oa  1)  38  TB  127 

J*10l 

KaM-32 

CALL  RACK  (J«K) 

127  CrtNTINUF 
Gfl  TB  00 
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j„l*S  PH B MAH  hHERF  pFBULAR  FRRH  »F  INPUT  ENDED 
F\C 


CHHMBN/7/NSETO?»3?*l2)»PARAMt  10*10  J»J*«ATEI  10 1*  TFI N( 101# ATRIBI • 
t  kKITE*  JHUN*  KBIJNT (101*  1SFED*  TNOA*  NBDN*  NBSECI14I*  PCSECM6) 
N  IS  HL'NaUP*  K  HAS  <*  NgRS  H«NF  eF  hH|CH  IS  N 
NDS4"0 
NMN 

CALL  CNjVRT  INN»NP*NCl 
K««K 

HANGUP  HAS^UCCUHRf^BECAUSE  BKLV  AVaILARLE  nrdfs  ARE  ALRFAOr 
i«HS  BF  IHR.NC)*  RR  ELSE  TmERE  ARE  NO  AVAILABLE  NODES 
bKlTLIlOM*?)  N*  K 

FfiRMAT  ( 1H0#?H1  ■  *  I  A*  ?HJ“*  1  ♦)  '  ***Z?^-^ 

Drt  10  IP.9*l2  C 

'  l-kset«m#nc#iP)  i 

IF  (L  *MJ.  0)  G«  T«  10  J||h 

CALL  CNVRT(L*LR«LC» 

if  «nset(lR»lc*si  *lt«  ai  afl  tb  i* 

o  CBKTINUF 

ALL  OF  HIS  *BRS  ARF  FU» L 

L  IS  A°MR  OF  N  hHB  H»  fEHF.R  THAN  a  NfiRS  |IF  AT  STATEMENT  14) 

A  KX*  12 

3  | aKsET (KR*KC*KX) 

IF  (I  *EG*  OJ  GB  TB  P7 
IF  cl.tn.LI  GB  TB  27 
CALL  CNVRT  «  I*IH.ID 

on  20  kb. 9* 12  .  ,, 

IF  <NSETIIR*IC*KGI  *E0.  L»  SB  Tn  27 
>W  CHKTINUF 
GO  TB  ?f. 

}7  KA.KX-1 

IF  ( Kk  .GE.\  91  GR  TB  19 


EHF.R  THAN  *  NB«S  (IF  AT  STATEMENT  14) 


IF  (Kk  *GE*\  91  GR  TB  1 
IF  (K"IN«1)1  a7»*8»AR 


47  K.K41 

GO  TB  49 

48  K«K»1 

49  CALL|CNvRT(K*KK*KC> 

DR  30  KP.9* 12 

IF  |NSET(kR»KC»KMI  *F0*  N1  GB  TB  46 
30  CltNTlNUF 
GR  TO  3 

CAN  exchange  branched  if  at  statement  ? 
as  IF  (NSET(NR*NC«3)  Ieq*  3)  ndsa.ndsam 
IF  (NSET(L«*LC*3)  *FQ*  3)  NDS4.NDS441 
CALL  REHROI KR*KC* IP* IC ) 


AT  STATEMENT  ?5 
3)  NPS4.N0S4.1 


1 


AS! 

NSETI*R#XC»3)-NSFT!XR#KC»3)-1 

*52 

NSET! IR. IC#3) -NSFT ( IR. IC.31-1 

*53 

CALL  PLACE !NR»NC#XH#XC> 

*5* 

CALL  PLACE! IR# IC.LR.LC ) 

ASS 

RETURN 

ASA 

75 

X  X- 1 2 

A87 

7A 

I -NSET !KR#KC#KX ) 

ASS 

WRITEI100#*! 

V 

ASS 

* 

FBRMATI5M0ST7*) 

1 

AAO 

fF  (I  »E0#  6)  GS  TB  7S  \ 

AA1 

CALL  CNVRT < I* JR# ICl 

AA2 

OB  77  XG-9#12 

AA3 

IF  INSET! IR#  IC#KG)  #20*  NN)  GO  TB  7S 

AAA  < 

77 

CrtNTIhUF 

AAS 

CALL  REORQIXR#XC»IH»IC) 

AAA 

NSETIKR#XC#3)-NS€T{XR#KC#3I-1 

\ 

AA7 

NSET! IR# IC#31-NSf T ( IR# IC#3)-1 

\ 

AAS 

CALL  FLACEINR«NC#XR«XC) 

1 

AAS 

CALL  PLACE  {NM«NC«1R#IC} 

A70 

IF  |NSETINR#NC#3)  .EO.  *)  ,  NDS*-NUS*-1 

*71 

GB  TB  7S 

A72 

78 

XX-KX-1 

*73 

IF  {XX  .GE.  S>  30  TB  7A 

*7* 

GB  TB  «* 

A7S 

7S 

RETURN 

*7* 

FRO 

*77 

C 

A7S 

c 

A7S 

SUBROUTINE  REBRDiNR.NCaXR.XC) 

AIO 

C0MMBN/7/NSET!3P#3?#12) 

AS1 

WRITE (100# 11 

*•2 

1 

format (13H0REBRD  Calleoi 

*•3 

WWITE!10M#5)  NR.  NC 

AS* 

5 

FORMAT ( ahonbde • . 1 2. ?X# I 2 I 

ASS 

WHITE! 100. 6)  INSET (NR.NC# 1 )» 1-3.12) 

ASA  , 

A 

FORMAT! 11H0ATH 103-1 2-# 10! 3X. U) ) 

AS7  \ 

WRITEIIOO.5)  «R*  Kt 

1 

ASS 

WRITE! 100# A)  INSET !KR#KC# 1 1# !-3#li 

r> 

ASS 

R-NR-32- (NC-1 1 

ASC 

K-KR-32- (XC-1 ) 

AS! 

C  RFBRC  PUTS  N  AND  K  AT  THE  END  OF  THE 

6TH 

AS2 

27 

M-NSET(NH.NC.3> 

AS3 

GO  TB  I0#2#3#*l#  NL 

AS* 

2 

IF  INSETINR#NC«10)  .P0»  XI  GB  TB  8 

ASS 

NSET {NR.NC#  S ) -NSC  f (NR#NC# 10 ) 

ASA 

GB  TB  8 

*S7 

3 

IF  (NSET(NR.NC«11)  *F0*  K|  GB  TB  8 

ASS 

IF  INSET INR#»C# 10) -X)  9# 10. 9 

ASS 

10 

NsET! NR.NC# 10 )*N£ET (NR.NC. 1 1 1 

700 

GO  TB  8 

NBRS 


701 

70S 

703 

70* 

705 

70* 

707 

7oa 

709 

710 

711 

712 

713 
71* 
715 
71* 
717 
71* 

719 

720 

721 

722 

723 
72* 
725 
72* 
727 
72i 
729 
73C 

731 

732 

733 
73* 
735 
73* 
737 
731 
739 
7*0 
7*1 
7*2 
7*3 
7** 
7*5 
7** 
7*7 
7*a 
7*9 
75C 


9  NSLT1NR.NC.9UNSFT (NR# NC«  1 1 1 
GO  TO  • 

*  IF  (NSETlNft.NC.12)  .CO*  K>  GO  TO  it 
IF  CNSETlNR.NC.il  UK)  12.13.12 
13  NSET (NR.NC. 1 1 UNSET (Nft.NC. 121 

r.n  to  h 

12  IF  1 NSE T I  NR.  NC«  1 0  UK )  15.J6.15 
1*  NSET (NR.NC. 10  UNSET (Nft.NC. 17) 

GO  TO  0 

15  NSE  T 1  NR.NC.  9 ) .NCFT 1 NR. NC« 1 2 1 
0  IF  (K  tl ii*  N)  GO  TO  19 
K.N 
HR.NN 
HT.NC 
NR.KR 
NC*KC 
GO  TO  27 
19  NM*HR 
NOMC 
RF.TUNn 
end 
c 
c 

SUbRbUT INE  CNvftT ( N. IRON. ICOL I 
00  10  U1.32 
IF  (N«32>  5.5.10 

10  N«N«32 
5  I  RB<*«\ 

ICOL- I 
RF  TORN 
ENO 
C 

c 

SUBROUTINE  PLACE ( NR. NC.KR.KC 1 
COhfiOK/7/NSET  ( 32. 32.  12 )  .PAR  An  1 10.  1  o  I «  JRATE I 
1  UNITE.  J»UN,  KOUNT  (  1i3).  I  SEED.  TNOrt.  NOON. 
LRaNft 
LCNC 
IR*aR 

ic*kc 

IF  ( (LR.E3. IR1  .AND.  (LC.E0.IC11  GO  TO  2* 
IF  1 (NSET1LR.LC. 3) «flf •* ) .OR. ( NSET ( IR, IC.3I. 
27  LEVEL.NSFT1LR.LC.3U1 
C  LEVFL  IS  Ml  OF  BRANCHES  ♦  1 
K.IR*32.(1C-1) 

GO  TO  M.2.3.*). LEVFL 

1  NSET1LR.LC.9I.K 
GO  TO  19 

2  NSETCLR.LC.10Uk 
GO  TO  19 


IOI.TFINi 10). ATRIOlO). 
NBSECll*!.  PrSECIlAI 


GF.*))  GO  TO  7% 
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781 

78* 

783 

784 
788 
78* 

787 

788 

789 
7*0 
7*1 
7*< 
7*3 
7*4 
7*8 
7«* 
7*7 
7*a 
7*9 

770 

771 
77* 

773 

774 
778 
77* 

777 

778 

779 

780 

781 
78* 

783 

784 
788 
78* 

787 

788 

789 

790 

791 
79* 

793 

794 
798 
79* 

797 

798 

799 

800 


3  N8fcTCLR#LC#ll )•* 

88  T8  14 

4  N8ET  C  LR#LC# 12 1 *8 

19  N*ETUR#LC#3I»LEVEI. 

17  (LR#EQ»KR»ANO»LC»EQ»RCl  88  TO  *4 

LR»*R 

lc»kc 

I  RaNR 
!C*NC 
68  T8  *7 
*4  RETURN 
END 
C 
C 

?ii!IS21i3J^fi?Ttl*.a8.1*»4R*RAH«10#10»#JRATEC10»«TFINC10)#ATRI8C8)# 

M0S4*U 

C  N0S4*N0  BE  N80C8  THAT  HAVE  4  N8R6 
08  1  l>l<3* 

08  1  J*li32 
N8ETCI«J«3»*0 
08  1  K«9«l* 

1  NSETCl#J#*l"0 
08  IS  NCal«3* 

OR  IS  NRal#32 

C  NC  IS  C8L  N8#  NR  18  RBW  N8 
N.NR43**|NC-1I 

C  NSET(NR#NC*3?,TELLr"h8H  HANV  NBM  N80E  * , 1 1 jT” ' S  ^ 
c  AT  END  BF  8UBR8UTINE  N8ETC If J# 3)»4  FOR  EVERY  N80E  CI#JI 

8  NRtpa4-N4ET(NR«NC#3l 
08  1*  Kf.al#NRER 

*  SS^INTcTFLBAlMloas-NOaAn'RHW^Ot**  SEARCH 

c  RICH  UR  NUH-TH  AVAILABLE  N80E  IN  NSCT  AnD  BEGIN  SEARCH 
MAaO 

*7  08  14  <CaNC*32 

?FI(8C#E0»Nf  I  .AND.  (8R»E0»NR  II  flB  T8  14 
C  AB8VE  TEST  PREVENTS  IT  FR8H  PLACING  BRANCH  8N  IT8ELF 
IF  (NSET(8R»KCf3la3l  7#7#14 
7  HAaNAal 

IF  JHA-NUHI  14# 4*9 

9  08  13  LNa9#12 

c  TCIT  CMMCTIN.  T*~ 

13  C8NTINUE 
NNR  a  NR 
NNCaNC 
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HOI 
1102 
MO  3 
MO* 
MOS 
MO* 
M07 
MOM 
MOS 
M10 
Mil 
M12 
M13 
Ml* 
MIS 
Ml* 
M17 
MIS 
MIS 
M20 
Mtl 
M22 
M23 
M2* 
*25 
M2* 
M27 
MSS 
MSS 
M30 

sai 

*32 

N33 

S3* 

S3S 

M3* 

N37 

N3S 

*3* 

M*0 

**1 

*** 

N*3 

a** 

M*B 

S»* 

**7 

N*S 

s*s 

NSO 


KKRaKR 

KKCaKC 

CALL  PLACECNNR.NNC.KKR.KKC) 

C  PLACE  UPDATES  MET  WITH  THE  NEW  BRANCH 
IF  (NSETINNR.NNC.ai.FQ.*)  NOS*aNOS**l 
IF  (NSETCKKR.KKC.3I.EQ.AI  NOSAaNOS**! 

GB  70  IP 
1*  CRNTINUE 

IF  <MA  .GT*  1S00I  GB  TO  100 

HA«1*00 

GB  TS  27 

C  ABOVE*  FAILED  TB  ADO  A  BRANCH*  MUST  TRY  AGAIN  AMONG  *|RST  NUM 
C  AVAlLABLF  NGOES 
12  C6NTINUF 
*3  CONTINUE 
GO  T6  IS 

C  MUST  PERFORM  BRANCH  FXCHANGE  IF  AT  STATEMENT  100 

100  KRITEaNStTCNR.NC.3l 

IF  (NSETCNR.NC.3i  .FO.  *>  GB  TB  IS 
CALL  DRANi)IRNUM| 

JalMT (Float IN )*RNUM) 

IF  (J  'Ed*  01  GB  TB  100 
IF  ( J  .EQ.  N)  GB  TB  100 
OB  101  KaS.12 

IF  CNSETCnR.NC.KI  .EQ.  J)  GB  Tm  100 
C  I.F..  IF  THEY  ARE  ALREADY  NORM*  TRY  AGAIN 

101  CONTINUE 

c  any  hangups  involve  *  br  fewer  nboes  ts  is  host  likely  nb.i 

WRITEC10R.2)  n.  j 

2  FORMAT ( 1H0*2HNa« l*«?HJa. I A ) 

CALL  XBRAN  (N.J.INCR) 

NOS*aNOM*«|NCR 

IF  (N8ET(NR.NC.3)  .EQ.  *)  GB  TB  IS 
IF  (KRIT£aNSCT(NR«NC«3>>  100.1*. 1* 

1*  WRITE! 10M.3S) 

35  FORMAT (MHONB  HOPE) 

IS  CBNTINUF 

WRITEH0M.37)  NOR* 

37  FORMAT  (1H  .2*HNB  OF  NB0E6  WITH  OEGREE  *a.l*| 

3  OB  G  1*1.32 
OB  M  w»1.32 

IF  (N8ET( I.J.3) <NE.*>  WRITE  CIO*. 10)  I.J.NSETc I. J.3) 

OB  S  KaS. 12 

IF ( NSETC t.J.KI.EQ.OI  MR ITE I 10B. 10 )  I.J.K 
S  CONTINUE 

10  FORMAT (7H0NB0E  C . 12. IN* . 12. HH)  HAS  ONLY  .11. SM  NORM) 

*  re’uRn 
eno 
SYMBOL 

OEF  ISITBN 
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MSI 

ISIT0N  CAlt*l 

C 

IU 

R0*0 

c 

111 

srcr 

3 

MS* 

CAL*«1 

1 

MSS 

tCtil 

* 

MSS 

LM«* 

13 

MS7 

LM«13 

lit 

MSM 

AW)«a 

•13 

MSS 

M 

2i* 

MSO 

CNO 
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Message  Transfer  HSfetwork 


1 

2 

3 
a 
S 

4 
7 
II 
9 

10 

11 

12 

13 

14 
16 
14 
17 
14 

19 

20 
21 
22 

23 

24 

25 
24 
27 
21 

29 

30 

31 

32 

33 

34 

35 
34 
37 
34 

39 

40 

41 

42 

43 

44 

45 
44 
47 
44 

49 

50 


C  MATS  PROGRAM 

COMMON  N*  IfttxT*  mult (4)#  STOMA.  SIGR.  RMU.  NST44).  LR(44*S>. 

1  NEXUi.2).  NXNC(AS).  HESSU).  LSM(44).  NRAN0( All).  LRASO. 

2  MftANC.  KM AND.  NBS.  JRATE.  *0(44.61.3).  MC*Q(A4.44> 

COMMON  101.  |j2»  163.  NSIOE.  NSQ*  ISEEC 

RFAO  (105.1)  LRANO.  MRANO.  KRANO.  NRASD(l) 

ON  20  TaP.124 

20  NRASD I t J • SR ASP ' T •! ) aP* I 

1  RFAO  (105.2)  SIGMA.  RMU.  SIGR.  N.  (MULT(K).K*1*4>*  JRaTE.  NSIOE. 

1  TSEFO 

2  FORMAT  (JF7.D.7I2. T5.I2.Il) 

NStiaNStOEaNSIDE 

ON  3  I *1.44 
NS( I )  «N 
LRM(I )«0 
NXKDdlaO 
ON  4  ja1.A4 

4  MCSGC  T . J) *0 
ON  5 

5  LS(  I* j) aO 
NFX(I.l)aO 
NFX( T.2)a0 
ON  4  tjal.bl 
ON  4  K«1 . 3 

4  NO(I.J.K)aO 

3  CONTINUE 
NXN0C45) aO 
Off  7  !»1.4 

7  MISSIIIaO 
NEXI45.1 )aO 
NFX(45.2)aO 
NffSaO 

I  FORMAT  (419) 

WRITE  (104.10)  SIGMA.  RMU.  RlQR.  N«  (MULT(K).Kal.4).  JRATE. LRANO. 
1  NSC.  ISEEO 
STGMAal. /SIGMA 
STOk*l./SIGR 
RMUal./RMU 
CALL  GASP 
GN  TO  1 

10  FORMAT  (INI.  AHST3MAa.F7.ff. 3X.3HMl|a. F  7.5. 3X.4HSIGRF a.F7.5.3X. 

1  PHNa. I2.3X.8HMUI  T( I ) a. A ( T2«2X).5HRATEa. T5.3X.5HRAN0a. T9. 

P  3X.4HNS0a.r2.3X.SHRREFa.il  ) 

STBp 

END 

C 

C 

SUBROUTINE  ORAND(RNUM) 

COMMON  N.  (NEXT.  MULT(A).  SIGMA.  SIGR.  RMU.  NS(44>.  LS(44.5). 

1  NEX(45.2).  NxnD(ar).  MESS(a).  LSM(ff«).  NRAN0(12R).  LRAsO. 
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SI 

SI 

CRM MSN  101#  I0»#  103#  NflOC#  NSQ 

S3 

LRAN0*LRANO«8B8.** 

St 

HRAM>bMRAND#3S5SAA33 

ss 

Jalsl ASS ILRANO 1/18777118 

St 

RNUMt  #S«ELSAT INRANO I J 1 «LRANO«MRANO 1 • *238030888- 9 

S7 

KRAN0aKRAN0a388s38089 

SS 

NS ANOlJI •SHAKO 

S9 

RETURN 

*0 

CnO 

*1 

c 

St 

c 

S3 

susseurtNC  cnvrti j.  iron# t col i 

SS 

COMMON  N#  INCUT#  MULTI SI#  StOMA#  SI OS#  SHU#  NSISSl#  LS(S«<5)« 

SS 

1  NEXtSB#tl#  NXNOISSI#  NCSSISl#  LSMISSl#  NRANOIItSI# 

LRAND# 

ss 

2  MRANO#  KRAND#  NRH#  JSATE#  NDISS«S1«3I#  MEsOISS#SS) 

S7 

COMMON  tot#  lat#  103#  NS10E#  N80 

SS 

MaJ 

89 

DO  ID  Iat#NS10£ 

70 

IE  IJ-NStDEl  5<S# 10 

71 

10 

JaJ.NBIOC 

71 

s 

ICOL* J 

73 

ISON* I 

7S 

JaM 

78 

RETURN 

7S 

CNO 

77 

c 

78 

c 

79 

SUBROUTINE  ROUTE  |t*tDE8T«IC) 

SO 

COMMON  N*  INCUT#  MULTlSI#  SIOMA#  SIBR#  RMU#  NSISSl#  LSISs<5l# 

St 

1  NEXI AS*tl#  NxnO|S§)«  NCSSISl#  LSMISSl#  NRANDl ItSI# 

LRAND* 

St 

2  MRANO#  KRAm fit  NOW#  JRATE#  N0lSS#St#3l#  MEoOISS«SS| 

S3 

COMMON  lot#  I0.«#  103#  NS1DC#  nrq 

SS 

taiaD 

SS 

CALL  CNVRTI 1# IRON# ICOL I 

ss 

CALL  CNVRTI IDEST# JROM# JCOL 1 

S7 

IE  IJROW-IROMI  1 #3*2 

SS 

1 

IE  1 I ROW- JR0N-NB I OC /2 1  S«S«B 

S9 

s 

IC-I-NSIDE 

90 

OR  TO  S 

91 

s 

IC# ISNSIOC 

9S 

GO  TO  S 

93 

t 

IE  1 JROw- IRON-NS I DE /2 1  S#5#S 

9S 

3 

IE  IJC0L-IC0L1  7*S*9 

95 

7 

IE  IlCOL-JCOL-NStDC/t)  ID# 10# 11 

9S 

10 

IC-I-1 

97 

00  TO  S 

9S 

11 

ICMst 

99 

OR  TO  S 

100 

9 

IE  | JCOL- I COL-NR I OF /2 1  UtlitlO 
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101 

• 

NR1TE  ( 10*# 12 )  |#  fOfST 

102 

12 

FORMAT  ( INI, 10X, 11MERR0R*R0UTE,  POX, (211 

ioa 

RETURN 

104 

* 

IF  (IC.OT.NSG)  ttatC-NSQ 

108 

IF  ( !C«LT*1)  IC.TONRQ 

10* 

IF  (151)  15*16,15 

107 

1* 

1GUIC 

10* 

IF  I JR0W*EO,IR8W)  RETURN 

'  09 

IF  I JC0L*IC0L)  3,1*, 9 

110 

15 

CALL  DRAN0  (RNUM) 

111 

IF  (RNUM,aT*0*S)  RFTURN 

112 

R»IC 

113 

1C* 161 

114 

IQlaK 

118 

It 

RETURN 

11* 

END 

117 

C 

111 

r 

119 

SUBROUTINE  INVERT 

120 

COMMON  N,  (NEXT,  MULTI*)#  SIQMa,  SIOR,  RMU,  NS(*4),  LKiAt.jti. 

121 

1  NEX(A5#2)«  NXND( *5) ,  MESSI4),  L8MC64),  NRANDI12*), 

LRANO, 

122 

*  MRANO,  KRAnO,  NOW,  JRATE,  N0(*4,Bl,3),  ME*G(*4,*4) 

123 

COMMON  |Q1,  Ig2,  lG3«  NSIDE,  NSO,  ISFED 

124 

lalGl 

128 

Ka|Q2 

12* 

IF  (K,EQ.l)  IfaUl 

127 

IF  (K«E0,2)  ICa |*1 

12* 

IF  (K,E0«3)  ICa laNRIDE 

129 

IF  (K*E0*4 )  ICaJ. NSIDE 

130 

IF  ( 1C,QT,NS0)  ICalC-NSO 

131 

IF  (  I C *LT a  1 )  |Ca |C*NgQ 

132 

IfllalC 

133 

RETURN 

134 

END 

138 

c 

13* 

c 

137 

SUBROUTINE  GASP 

13* 

COMMON  N,  INEXT,  MULTI*),  SIOMA,  SIGR,  RMU,  NSI64),  I.RfAnSJ, 

139 

NtX  I  A5,P ) ,  NXND(AS),  MESSU),  LSM(64),  NR  A  NO  ( 12S), 

LRANO, 

140 

2  MRANO,  KRAnO,  now,  JRATE,  nD| *4,51,3),  MEs6(*4,*4) 

1*1 

COMMON  lei,  lG2*  Ig3,  NSIDE,  NRO 

142 

INEXTaO 

143 

00  S  I » 1 , NSO 

144 

call  oranoirnuM) 

14S 

NO!  l,51*l)a*lNT (RIGMAaALOG(RNUM) )♦! 

14* 

NEXI |,l)aNO( 1,51,1 ) 

1*7 

NEXl I,2)a5l 

14* 

IGlal 

149 

CALL  OROERN 

ISO 

5 

CONTINUE 

191 

192 

1 

193 

4 

194 

199 

3 

196 

197 

192 

30 

199 

160 

161 

14 

162 

163 

164 

169 

166 

167 

168 

169 

32 

170 

33 

171 

34 

172 

35 

173 

36 

174 

37 

179 

38 

176 

177 

10 

178 

179 

180 

23 

181 

182 

183 

24 

184 

26 

185 

186 

26 

187 

27 

188 

28 

189 

29 

190 

22 

191 

192 

193 

194 

6 

195 

196 

7 

197 

198 

199 

200 

S 

JTIMEaO 
CALL  EVNT8 
jTIMEaJTIMEal 

IF  ( JTIHC.LT. JRATE)  08  TB  10 

nbdn*o 

08  ao  I *1»N8Q 

IF  (NS(I)*E0jC)  NBDNaNBONal 
CONTINUE 

OCTaFLBAT (  NBON 1 /FLOAT <  NSQ 1 
WRITE  (108*14)  NAN*  NBON*  ACT 
FORMAT  (lM0///4X.5MTIMEa*I6*3X*5HNBDNa 
IF  (NSI0E*E0*2)  WRITE! 108*38)  (NSC  I »# I 
IF  (N8I0E*E0*3)  WRITE! 108*33)  (NS! I )» I 
IF  (NSI0E*E0*4)  WRITE! 108*34)  (N8(1)«I 
IF  (NSIOE.EQ.S)  WRITE (108* 38)  (NS(!)*I 
IF  (NSI0E.EQ.6)  WRITE! 108*36)  (N8(I)«1 
IF  (NSIDE*E0*7)  WRITE! 108* 37)  (N8! I )« I 
IF  (N8!0C*£0*8)  WRITE! 108*38)  (N8(I)«I 
FORMAT  (2(2(4X«!2)7) ) 

FORMAT  (3!3(4X*lf)/)> 

FORMAT  ( 4 ( 4(4X* 18) ✓) ) 

FORMAT  (S(S(4X« 18)7) ) 

FORMAT  (6(6(4X* 18)7) ) 

FORMAT  (7(7(4X«l8)/)> 

FORMAT  (R(fl(4X* 18)7) ) 

JTIMEaO 

CONTINUE 

ITEST»ISIT8N|1) 

IF  (ITE8T*NC*1)  00  TR  22 
ITESTalSITONIl) 

IF  (ITE8T.NE*0)  00  TO  23 

IF  (MULT(l)-l)  84*88*24 

OR  26  Ial*5 

MULT!  Dal 

00  TO  28 

00  27  !al*S 

MULTdlalO 

WRITE  (108*29)  MULT ( 1 ) 

FORMAT  ( 1H0»1SX* 10MMULT( 1*5) ■* 14 ) 

CONTINUE 

I  TEST  a  1 SI TON ( 4 ) 

IF  | !TE8T*NE*4)  00  TB  803 
RFAO  (101*6)  I 
FORMAT  (12) 

WRITE  (100*7)  I 

FORMAT  UM0»17HC0NTENTS  OF  NBDF  *12) 
OR  800  llal*51 

WRITE  (108*801)  II*  NO  11*11*1) 

NZaNOl 1*11*2) 

LH*4  NZ 


I2*3X.4nRCT*,F5.3//> 

1*4) 

1*9) 

1*16) 

1*20) 

1*36) 

1*49) 

1.64) 
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901  >  ITU#*  ILENTH 

902  S  LI* 1  2 

903  >  LI#*  NM 

90*  S  ITU#*  IDENT 

90S  S  LI# 1  3 

90*  S  l.t  #*  NX#  1 

907  t  STU#*  IC 

90S  UNITE  ( 101# *09)  ILFNTH#  IOEIT#  IC 

909  N/aNOl I#I1#3) 

910  S  LH#*  NZ 

911  S  STU#*  I IUCC 

919  •  LI# 1  1 

913  S  LH#*  NX# 1 

91*  S  STh#*  I QUEUE 

91S  UNITE  ( 101# *09)  tIUCC#  IQUEUE 

91*  100  CONTINUE 

917  101  FORMAT  I 1H0#2X# I9#3X# If) 

911  109  EONMAT  C1H  #7X#3I9) 

919  103  CONTINUE 

990  ITEST*ISIT0M2> 

991  IF  ( ITEIT (NE*9)  SO  TO  11 

999  19  I TEST* I SI TONI 9) 

993  IF  (ITESTtNEtO)  SO  TO  19 

92*  RCTUAN 

925  11  CONTINUE 

92*  I TEST* I SI TONI I) 

927  IF  lITESTiNE#*)  00  TO  1 

921  2  ITEIT* I SI TONI!) 

929  IF  <!TEST.NE*C)  00  TO  2 

230  UNITE  1 101# 90) 

931  UNITE  1 101# 21 )  (IMFRQ(I#J)#J*l#*«)#!al#**) 

232  RETURN 

233  20  FORMAT  I 1M1,S*X#22HHE8S  OEN  <8R10!N,DEIT) ) 

93*  11  FORMAT  I 1H1«5*X#22HL!NE  BLOCKING  DISCARDS) 

93S  21  FORMAT  I  *1391 IX# |3)/)//) 

93*  1*  FORMAT  I  lil*f3x#!S)/>//> 

237  IS  FORMAT  I lMl#S*X.23MN00AL  BLOCKING  DISCARDS) 

23*  ENO 

939  C 

2*0  C 

2*1  SUBROUTINE  EVNTR 

942  COMMON  N(  I NEXT,  MULTI*)#  SIOHA#  SION#  RMU#  NSI**)«  LSI**#S)« 

2^3  1  NEXI AS«2)«  NXNDIAR)#  MESSI*)#  LSMI**)«  NRANDI12D#  LRANDf 

9**  2  MRAnD#  KRAnDi  NAN#  JRATE#  nDI**#S1«3)«  MEftU (**«**) 

9*S  COMMON  tfll#  1q2<  I Q3<  NSIOE#  NSO 

?**  I ■ I NEXT 

2*7  NAUaNEXI 1# 1 ) 

2*1  INExTaNxNDII) 

2*9  NXNDIDaO 

2S0  JaNEX 1 1 #2) 
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r 


251 

c 

THIS  It  HCS2A0C  NO  IN  OUCUCIII  WHICH 

252 

IF  <J*£Q*0)  MITE  ( 102*999) 

HI 

999  FORMAT  (IHla  7HERREVNT) 

254 

c 

TIC  U2  LIST 

255 

NZaNDl 1 aja3l 

rat 

2 

LH«4  NZ 

a 

STM#  4  I  SliCC 

2sa 

2 

Llal  1 

>55 

2 

LHa4  NZal 

MO 

a 

8THa4  {QUEUE 

Ml 

NEX( 1*2) • ISUCC 

MS 

NO ( I a  Ja  3 ) ■ I QUEUE 

M3 

IF  ( ISUCCtNCtO)  OB  TO  5 

2*4 

NCXIIaltaO 

2*5 

00  TO  1 

Mi 

5 

NCXI la  1 1 aNO( la iSUCCal 1 

2*7 

1 

CONTINUE 

n 

IF  ( JtNCaSl 1  00  TO  12 

2*9 

lOlal 

S70 

CALL  HCSSAQ 

S71 

IQSaSl 

272 

CALL  0ROCRQ 

273 

IF  l(Ntin*CQ.OI*00.(L2Hm.NC«on 

274 

00  TO  2*8 

S7S 

12  CONTINUE 

27* 

NZaN0(IaJa2l 

277 

2 

•  LH#4  NZ 

27* 

a 

STWa4  ilentm 

279 

s 

Llal  2 

220 

a 

LOa  4  NZal 

221 

a 

8TW«4  I0C8T 

222 

a 

Llal  3 

223 

a 

LB# 4  NZal 

224 

a 

*Th#4  1C 

225 

c 

1C  18  NEXT  NODE  TO  WHICH  HC28A0E  MUST 

22* 

c 

I0F8T  18  ITS  FINAL  DESTINATION 

2*7 

ICSalC 

222 

IF  ( ILCNTHtCQiO)  00  TO  400 

229 

c 

ILENTHaO  HEANS  DEPARTURE 

290 

IF  (LBHI I 1 tCQtJI  LSMIIlaO 

291 

IF  (I0C2T*EQ*n  00  TO  105 

292 

c 

NOW  WILL  SEE  IF  THIS  IN  A  NCTRV 

293 

00  2  Kal<4 

294 

IF  (LSI laK) tCQtaJ)  00  TO  3 

295 

2 

CONTINUE 

29* 

00  TO  4 

297 

3 

IF  (NS( ICItCQtO)  00  TO  150 

292 

c 

NSllClaO  IF  NODE  1C  IS  BLOCKEO 

299 

00  TO  200 

300 

4 

CALL  ROUTEl la IOFST# ICl 

14C 


301 

N! 

303 

30* 

308 

30* 

307 

308 

309 

310 

311 
313 
313 
31* 

318 
31* 
317 
311 

319 

sro 

381 

388 

383 

38* 

388 

38* 

387 

388 

389 

330 

331 
338 

333 

334 
338 
33* 

337 

338 

339 

340 

341 
348 

343 

344 

345 
34* 

347 

348 

349 

350 


KCMaO 

*  NOI I# J«8laND( I* J*8)»IC84]C 
J2MC-1 

IF  ( I jZtEOi  1 I •  OR# I JZ*EQ* («N8Cal 1)1  K»1 

IF  l(JZ#E0*-l>*6R#(JZ*EQ*(  N90>im  K»8 

IF  (<jz»ca«  NSMEItOR*  I  JZ#EO«  (•NSQaNStOE) ) )  K#3 

IF  I (uZaEOtaNSlOFI #0R# I JZtEO# I  NSQaNStOEIII  K«4 

IF  (K#LT«1#9R#K#0T#4>  06  TO  197 

NOI 1* J#3>  *0 

06  T6  199 

197  WRITE  ( 108# 198)  K 

198  FORMAT  ( 1H1*7HERR0R*X« 3X# 2HKa# IS I 
RETURN 

199  CONTINUE 

IF  (NOI ICItNEtOl  00  T8  7 

IF  (KCH#CQ#1 I  00  TO  8 

IF  (KCM#C0#8|  00  TO  7 

KCHal 

ICSalC 

IC4I01 

00  TO  * 

8  K-Ul 

IF  (K*QT#N8QI  KaK»N80 
IF  (N9IKI#EQ#0I  OR  TO  9 
17  ICSalC 
IC*K 
OR  TO  * 

9  Ratal 

IF  (K#LT#1)  KaKaNRO 
IF  (N8(K|#EO#OI  00  TO  10 
OR  TO  17 

10  KalaNSIOC 

IF  (K#0T#N80I  KrK>N80 
IF  |N8(K)#E0#0I  00  TO  11 

00  TO  17 

11  KalaNSIOC 

IF  (K#LT«1>  KaKaNOO 
IF  (NSiKitNEtO)  00  TO  17 
KCHaS 
OR  TO  17 
7  CONTINUF 

IF  (LSI l#K).NE*0>  OR  TO  100 
IF  (NS( 1C) ■NE*OI  OR  TO  800 
LS  I  I«K)  ar>J 
ND(I#J#3)a0 
OR  TO  150 
800  LS  (I#K»eJ 
IS  ITIMEaNOWalLENTHaMULT(K) 

14  ND(I#U#l)alTIME 
810  IF  (K#EQ#S)  00  TO  850 
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SSI 

MI 

M3 

18* 

MC 

38* 

M7 

MS 

Hi 

3*0 

3*1 

3*3 

3*3 

3** 

3*8 

3«« 

3*7 

3*8 

3*9 

370 

371 
378 
373 
37* 
378 
37* 

377 

378 

379 

380 

381 
388 
383 
38* 
388 
38* 

387 

388 

389 

390 

391 
398 
393 
39* 
398 
39* 

397 

398 

399 
•00 


e 


c 


8 

8 

8 


08  811  JZ*1«N 

18  !N0( IC« JZ«8) aEQvO)  88  T9  818 
811  eSt./INUE 

NRITE  (108*8131 

813  E8RMAT  (1H1«1*HCRR8R  IN  CVNT8I 
RETURN 

JZ  18  THE  NUH8ER  87  A  FREE  8RACC  IN  THE  QUEUE 
818  N8  (IC>aNB(IC<-l 
NO!  |C« JZ*1 )*ITIME*1 
NOC 1C* JZ*8I*K0< !« J«8I*IC 
no  <ie«jz«3i*o 

MSTORaNEXI IC*1) 

IQIalC 
I88*JZ 
CALL  8R0ERQ 


IE  <N8T8RtEQ*N£x<IC«in  88  TO  1,8 
REMOVE  IC  ER8H  LINKED  LIST  AH9N8  NODES 
IE  (INEXTtEQaO)  88  T8  819 
IE  ( ICaEQa I NEXT)  88  T8  819 
NSTBRalNEXT 
08  81*  Ml*l*«8 
IE  (N8T0R*E0*0)  09  T8  818 
IE  |M8T8R*E9alCl  88  T9  817 
H8T0Rl*MST0R 
81*  H8T0R *NXND  C  H8T9R I 
88  T9  818 
819  INExT*NXNO(l£l 
NXND(ICI*0 
88  TO  818 

817  NXNO( H8T8R llaNXNO(IC) 

NXNOdClaO 

818  lOlalC 

call  oroern 

818  CONTINUE 
880  NZ*N0<I*J«8I 

LI*  1  1 

LH«*  NZ* 1 

8TN«*  NZ 


ND(l«J«8|aNZ 
NZaNDl I* J«3) 
a  Li* i  i 

8  LH«*  NZ*1 

8  8Tn«*  I QUEUE 

NO ( I «  J«  3 ) a I OUEUE 
181*1 
I82aj 

CALL  8R0ERQ 
300  181*1 

CALL  OROERN 
return 
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\ 

\  : 
I  ' 


•01 

•88  CONTINUE 

•0! 

, 

IOESTaHM8(8) 

•01 

tl.£NTH*MFSS(l) 

•0* 

MCSG  ( 1 « IDEST ) *MFRG ( 1 « 1 0E8T ) *1 

•08 

c 

8 

»a>l8*06*S18  i 

•08 

1 T  IHE*N0W*  ILENTM*MULT  C  8 I 

•07 

08  88  JZ*1*N 

•08 

IF  (N0( 1* JZ«2)*EQ*0)  00  TO  70 

•08 

88 

CONTINUE 

•  10 

WRITE  (108*711  1 

•11 

71 

FORMAT  (1H1*7H£RR0R-A*8X*IS) 

•18 

RETURN 

•  11 

70 

N0( 1«JZ«1)*TTIME 

•  14 

8 

LN«8  ILENTM  \ 

•18 

8 

1  s»tm*b  .  •  1 

•1* 

■ 

LM«B  I0F8T 

•17 

8 

L!*l  8 

•  18 

6 

STS*  5  •*  1 

•19 

_  8 

STW«8  NZ 

•80 

NO! 1* J2*R)*NZ 

•81 

NO  (!*J7*1>*0 

•88 

\ 

NS(n>NS(l)*l 

•81 

1 

L8M( 1 )*J2  \ 

*8* 

1 

IS1  •  1 

•88 

IQ8*JZ 

•8* 

CALL  ORDERS 

•87 

00  TO  >100 

•88 

105  tC*S  1  ,  \ 

•89 

N0( I«j«lla0 

•10 

IF  (LSI l*KI *EO*OI  00  TO  200 

•11 

100  LSNER*LS(I»K>  ' 

•18 

IF  (LS(|«K)*LT*0)  LSREP*-I  i(l#KI 

•11 

KC*0 

•!• 

108  NZ*ND( t«LSREP*3) 

•IS 

KC*K£*1 

•18 

IF  (KC*OT* (N/8+1 ) * AnD*K*NE*5)  00  TO  110 

•17 

8 

Ll«l  1 

•18 

8 

Lh*4  N7*l 

•19 

S 

\ 

Sll »«•  IQUFUF 

••0 

i 

IF  ( ILUFlJEaEOaO)  GO  TO  107  , 

••1 

LRRERalQUEUE 

••8 

GO  TO  108  , 

••1 

107 

mt  |*LSRFP*3)*M)(I,I  8REP*a>«J 

••• 

NO  (I*w*l)a0  ' 

••S 

ISO 

ND(  I«<j*  0*0  , 

••6 

00  TO  100 

••7 

110 

K*X*1 

••8 

IF  (X»GT»*  i  K*1 

••9 

101*1 

•50 

rr.8*n 

\ 


\ 

! 

'  •  t 
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CALL  INVENT 
ICSvIC  1 
I C* 101 

NO) I# J#F)*N0( I# J»2I*ICS*IC 
•6  T8  7 

400  NUAIT*!0U£U£ 

OS  401  K*i«S 

IF  (LS( I#K)  #EQ# J)  OR  TS  401 

401  CONTINUE 
UNITE  (101*402) 

402  FORMAT  (1H1#*MERR8R2) 

NETUNN 

403  noci*j*i)*o 

NOI  I# J#7)*0 

NO!  I# J#3)*0  i 

LI  I  |#K)*0 

NSl I ) *NS( I ) ♦ 1 

IF  (NSm.NC.l)  SO  TO  440 

10141 

CALL  NCTNY  , 

440  IF  (NUAtT*NL#0)  SO  TO  4S0 
SO  TS  300 
450  JaNUAIT 

t*  (K.E0.5)  SO  TH  A AO 
IF  (NSMCI'STtO)  GO  TO  410 


I*  (K#EQ#5)  < 
IF  (NSC  tCt  *01 
LS(I«KI**J 
SO  TO  300 
4*0  N2*N0( I#J#2) 
LH*4 
ST4#4 
OR  TO \ 200 


N2 

ILENTH 


SUBROUTINE  RETRY  , 

COMMON  N«  I NEXT#  MULT(*I«  SIAM*,  SIOR#  RMU#  NSC**)#  LS(64#S)« 

1  NCXC *5<2)#  NXN0CA5)*  MESNI4I.  LSM( **)#  NNAN0(12l)«  LRANO, 

2  MRAND#  KRAnO#  NOW#  JRATE.  N0(*4«5l«3)f  MESGCA*.**) 

COMMON  IG1#  102#  103#  NSIOE#  N8Q#  ISCEO 

01  MENS  I  ON  Jl  (4I|#  K1  (*)#  JF  (41  #  KF(4> 

1*101 

00  10  L*l«4 

JF(LI*0 

KF(LI*0 

J1(L)*0 

K1ILI*0 

J*I*i 

IF  ( J*E0#0)  J*N80 

\  K*1  ■ 

I  M*0  \ 


I 


SOI 

sot 

soa 

504 

505 

so* 

B07 

SOI 

SOS 

510 

511 
sit 
sia 

514 

515 
SI* 
517 
SIS 
SIS 

sto 

Stl 
St  2 
sta 

524 

S8B 

St* 

S87 

StS 

StS 

sao 

sai 

sat 

saa 

sa* 

sas 

sa* 

sa7 

saa 

sas 

540 

841 

S4t 

S4a 

544 

545 
54* 
547 
54* 
S4S 
SSO 


IF  (LSI J*K) •QC'O)  OS  TS  1 
Mai 

JlIHlaJ 
Kl(H>aK 
I  B«1 

OS  TS  SO 

1  Jalal 

IF  (JtCOdNSQain  Jal 
Kat 

IF  (LSI J*K) »SC*0)  OS  TS  t 

MaMal 

Jl(M)aJ 

Xl(M)aK 

fS>t 

SS  TS  tO 

2  Jal.NSIDF 

IF  (J«LT*1)  JaJaNSO 
Ka3 

IF  (LSI J*<) •6E*0)  OS  TS  a 

MaMal 

JIIHIaJ 

KllHlaK 

iso 

08  TS  80 
a  jalaMSIK 

IF  (J.OT.NSO)  J.J-NSO 
Ka4 

IF  (LSI J*K ) aOCaOl  G8  TS  4 
MaMal 

JllMlaJ 

XKMIaK 

1S*4 

00  TO  80 

4  IF  (M.EO.O)  RETURN 
CALL  ORANO(RNUM) 

IF  (MF.OT.O)  08  TO  S 
M INT ( SNUMaFLS AT ( M 1 1 al 
JaJlIHai 
KaKltMai 

CALL  ORANOIRNUMI 
BIGRlaSlGR/FLBATIMt 
08  TO  * 

5  M3« INT ( RNUMaFLBATI MF t  )al 
•laJF  (M3! 

KaKF(M3t 

CALL  ORANOC'NUH) 
SIGRlaSIOR/FLSATIHF) 

*  |RETRVa«f ttT ( SlORlaALSGI RNJM) taiaNBM 

M4«*LS( J«K I 
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r 


BSI 

BS2 

ssa 

SB* 

BBS 

SB* 

BS7 

BBS 

BBS 

S*0 

B*1 

B*2 

s*a 

S*4 

B*S 

B** 

B*7 

B*B 

B*9 

B70 

S71 

B72 

B73 

S74 

B7B 

B7* 

B77 

B7S 

B79 

810 

BSI 

BS2 

BS3 

BB4 

BS6 

BS* 

BS7 

BBS 

BBS 

S90 

B9I 

B92 

B93 

B94 

S9B 

B9* 

S97 

S9B 

B99 

*00 


NRAVEaNDI J*N4* 1 ) 

NO (J*M4*l)a| RETRY 

If  ( NSAVE«EQ«C)  00  T9  21B 

NZ*NO( J*N4*3) 


S 

LH*  4 

NZ 

B 

STh*  4 

1DBUC 

S 

Ll*l 

1 

1 

L>*4 

NZ*  1 

s 

su*  4 

IOT 

If  ( NfcX  ( J*2 1 «E0»0)  06  TB  219 

If  |N4>C0>NEXIJ«2> 1  08  TB  219 

N8T8RaNEXIJ*2) 

oa  214  NlalfN 

If  (HSTBRtEO.O)  Oa  T8  21A 
If  (N£TBRtE0*N4)  OB  TB  217 
HSTBRlaHSTBR 
NZaND<J*N6T8R«3) 

S  LH*4  NZ 

S  STW*  4  HBT8R 

214  CONTINUE 
OR  TO  21B 

219  NEXI J*2) alDSUC 

If  < IDSUC»E0»0)  09  TB  220 
NEX(J*l)aNOlJ*lMUC*l) 

OR  TB  221 

220  NEXIJillaO 

221  ND<J*N4,3)aI0T 
OB  T8  21A 

217  NZ*ND( J»N8TBRl*  31 


S 

Li*l 

1 

s 

LH*4 

N7*l 

s 

LM*S 

n  sue 

a 

BTH*B 

4 

B 

STh*4 

NZ 

NO ( J* NSTBRl «  3 1 aNZ 
N0( J*N4*3 1* 10T 
21B  NEXaNLXI J*1 1 
lOlaJ 
102aM4 
CALL  BROERO 

If  |NEX.E0.NEX|J*1 ))  OB  TB  230 
If  ( J«EQ> I NEXT)  OB  TB  719 
If  IlNEXT.EO.C)  GO  TB  719 
NfiTBRalNCXT 
OB  714  MlaliAS 
If  (NSTBR«£Q*0)  OR  TB  71H 
If  (NSTBRtEO* J)  OB  TB  717 
NSTBRl aMBTBR 
714  NlTBRaNXND I NSTBRl 
OB  TB  711 
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01  719  INEXTaNXNO(J) 

02  NXNOfJIaO 

09  60  TO  71 • 

04  717  NXND(MST0R1»bnxN0(J» 

OB  NxMHJIaO 

00  710  tiilaj 

07  CALL  ORDERN 

00  090  CONTINUE 

09  RETURN 

10  <0  ITa.LSCJaKi 

11  IE  (ISEFOalQaC)  60  T0  (1*2*9*41*16 

10  NlaNOCJf IT#OI 

19  6  L  I* 1  2 

14  8  LB»4  N2a  1 

IB  6  STtaa4  IOERT 

1*  IE  ( lOESTaEQt 1 1  60  TO  21 

17  60  TO  < 1*2*9«4) a  !R 

10  01  HEanEal 

19  JECNE laj 

20  KECmE la« 

21  60  TO  ( 1*2*9*41 <  IS 

22  ENO 

>29  r 

>24  C 

>25  SUBROUTINE  OROERO 

26  COMMON  N,  I  NEXT  <  MULT(6>«  SIGMAa  SIOR.  RMU.  NS(64)a  L8I64.5I* 

27  1  NEX(6K*2>«  NXNO (40)*  MESSIAH  LSM(64>*  NRAND<120>«  LRANOa 

20  2  MRAND.  KRA*Oa  NflWa  JRATE*  N0(64*51*9>*  ME00(44*44) 

29  COMMON  |Q1«  I 02a  lG9a  NOIOEa  NOO 

90  laid 

91  JalG2 

92  IF  (NOIlaJall aEQaO I  RETURN 

99  NZaNOI I aja9> 

94  6  L  la  1  1 

96  S  LH«4  NZal 

96  0  STtoa4  I QUEUE 

97  MSToRaNEXI I #  2  > 

90  IE  (MSTORaEOaOl  60  TO  12 

99  IE  (ND( la MOTOR# 1 ) *GE#ND( I* J* 1 ) )  GO  To  14 

40  00  9  Mala51 

41  IE  (MSTORaEQiO  60  TO  11 

42  IE  I N0( I a MOTOR* 1)« BE a NO (I* Jail)  00  TO  10 

49  MBTORlaMSTOR 

44  NZ«NO( JaMST0Ra3» 

45  0  LHa4  NZ 

46  S  STn*4  M8T0R 

47  9  CONTINUE 

40  NRITE  ( 106*21 > 

49  21  EORMAT  { lHlal2HFRR0R>9R0ER0) 

50  NRITE  (106*1001  la JaNEXI Ia2l 
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Ml 


M> 

Ml 

*14 

MS 

M* 

M7 

MS 

*•9 

MO 

Ml 

MI 

MS 

M4 

MS 

MS 

M7 

MS 

MS 

*70 

*71 

*71 

*7S 

*7* 

*7S 

*7* 

*77 

*70 

*79 

MO 

Ml 

*S! 

MS 

M* 

MS 

M* 

*S7 

*ss 

*S9 

*90 

*91 

*91 

*9S 

*9* 

*99 

*9* 

*97 

MS 

*99 

700 


0*  101  MF«1«B1 

WRITE  (108*101)  MF«NOfl*NF«ll 
NZ*N0(I*HF«1) 

S  LM*8  NZ 

S  8TW*B  ILF 

S  Li* 1  I 

S  L0*8  NZ*1 

S  8TN*R  IOC 

S  LI*  t  S 

■  LS«S  N?»1 

S  8TN*B  ICC 

MRITC  I10S*10SI  ILF •  IOC*  ICC 
NZ»NOll«MF«S) 

S  LM.O  NZ 

S  STN.S  ISS 

S  Li* 1  1 

s  lm#b  NZ«1 

S  STw**  100 

101  MRITC  ( 108*104)  ISS* I 00 

100  FORMAT  (»X«3(IS«SX)I 

101  FRRMAT  ISX*  I1«SX«I9I 
10S  FORMAT  (10X*S( 18*1X1 1 
10*  FORMAT  |10X*llIS«IXt> 

return 

11  NCX( 1*11 aNOl I* J* 1 I 
NCXI I*|laJ 
NOC I* J* J)»I0UCUF 
mtuRm 

11  N0II*J*3|4l0UCUC 
80  TO  IS 
10  CONTINUC 

S  LM«4  I0UFUF 

S  LN«B  MOTOR 

•  sth*s  * 

S  STM**  NZ 

N0(I«J*S>«NZ 
IS  NZ"N0I1«MSTSR1*SI 
S  LI*1  1 

S  LH«*  NZ«1 

S  LM«B  J 

S  STM* 9  * 

S  STM**  NZ 

N0lI«MST0Rl*SlaNZ 

RETURN 

1*  NFX( I* 1 ) aNOl I*J*1) 

NFX(I«lt«J 

•  lm«*  i Queue 

•  LM«R  MOTOR 

S  STM* s  * 

•  STM«*  NZ 
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701 

N0(  I>U>3)aN7 

702 

m  Turn 

703 

fnd 

70* 

C 

70S 

C 

704 

SUBROUTINE  OROERN 

707 

COMMON  N*  I NEXT*  MULTI6)*  SIGMA*  SIGN*  RMU*  NS(6*)«  LS(6*«5)« 

701 

1  NEX(65«2)«  NXNDIAB)*  M£SS(*t«  LSM(6*)«  NRAN0C12*)* 

LRAND* 

709 

2  MRANO*  KRAnD*  NON*  URATE*  ND( 6*«51« 3 ) «  MERG(6*«6*) 

710 

COMMON  IG1*  IG2*  IG3*  NSIDE*  NSO 

711 

1*101 

712 

IF  (NEX(I*1)»E0.0)  RETURN 

713 

MOTOR* I NEXT 

71* 

IE  ( INEXT *E0*C )  GO  TO  12 

715 

IE  INEXI MOTOR* 1 )*GE*NEX(I«1I)  GO  TO  1* 

71* 

DO  9  M*1*6S 

717 

IE  (MST0R*EQ*0)  GO  TO  11 

71* 

IE  (NEX(MSTOR«ll.GE*NEX(!«llt  GO  TO  10 

719 

MSTORl-MSTOR 

720 

9 

MSTORaNXND(MSTOR) 

721 

WRITE  (10**21) 

722 

21 

FORMAT  (1H1«12HERROR*0ROERN) 

723 

RETURN 

72* 

12 

iNEXTalfll 

725 

NXN0II)*O 

72* 

RETURN 

727 

11 

NXN0IMST0R1I-IG1 

72S 

NXNOI I )*0 

729 

return 

730 

10 

NXNOIMSTORl ) *IG1 

731 

NxNOI I  )*C1ST0R 

732 

return 

733 

l* 

INEXT-IG1 

73* 

NXNOI  DaMSTOR 

735 

return 

736 

END 

737 

C 

73* 

c 

739 

SUBROUTINE  MESSAG 

7*0 

COMMON  N«  INEXT*  MULTI 6)*  SIGMA*  SIGR*  RMU*  NS(6*>*  LS(66*5)« 

7*1 

1  NEX ( A5*2 ) *  NXNOl AS) «  MESS!*)*  LSMI66)*  NRANDI12*)* 

LRAND* 

7*2 

2  MRANO*  KRANO*  NOW.  URATE*  ND(A**5l*3)*  MESG(**«**) 

7*3 

COMMON  TGI*  152*  IG3*  NSIDE*  NSO 

7** 

I  - 101 

765 

CALL  DRANO(RNUM) 

766 

NO ( I «  51  * 1 ) 'NO*. TNT ( SI GMAaALOG ( RNUM ) ) ♦ 1 

767 

|E  <(NSU>*EQ*0)*0R*(LSM(!)*NE*0))  RETURN 

76* 

1 

CALL  DR AND (RNUM) 

7*9 

JaINT (RNUMaELOAT(NSO) >*1 

75C 

IE  ( I *EQ*U)  GO  TO  1 
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751 

752 

753 

754 

755 
75* 
757 
751 
759 
7*0 
7*1 
7*« 
7*3 
7*4 
7*5 
7** 
7*7 
7*4 


t  CALL  ORANLIRNUHI 

ILlNTH»»INT<RHU#AL5a<9NUM>  1*1 
HEI8I2NJ 
HF.SSOlalLENTH 
RETURN 
FNO 
SYMBOL 

DEE  ISIT5N 

ISITBN  CAL2.1  0 

R0»0  0 

STCF  3 

CAL2j1  1 

SCS«3  4 

L*#4  13 

LN«13  1*4 

AND*  3  *13 

b  2*4 

END 


150 


In  the  main  body  of  this  work  we  consider  M/MA  queueing  systems, 
i.e. ,  stochastic  service  systems  which  experience  Markovian  arrivals  and 
in  which  customers  depart  after  receiving  an  amount  of  service  time  that 
is  exponentially  distributed  and  is  given  by  one  of  k  servers.  If 
there  are  n  customers  presently  in  the  system,  then  a  customer  will 
arrive  in  the  next  instant  of  time  At  with  probability  AnAt  +  0(At) 
and  a  customer  will  depart  in  the  next  instant  of  time  with  probability 
PnAt  +  O'At) . 

The  stationary  probability  of  finding  n  customers  in  the  system 
is  related  to  Pg  =  P  [erpty  system]  in  the  following  way: 


which  is  valid  for  all  n  >  0  if  we  define 


it-1- 


Then  is 


found  from  the  fact  that  if  this  is  to  be  a  valid  probability  distribu- 


E  pn  =  1 

n=0 


If  =  A  and  V>n  ®  M  for  all  n,  then 


p»  "  P°  v  =  P0(M> 


therefore 


p  =  (1  -  p)pn  for  p  <  1  where  p  =  ^ 
n  y 


is  called  the  "utilization  factor" 
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p  =  1  -  =  P[  system  is  busy] 

For  an  infinite  server  system  every  customer  has  his  own  server. 

We  take  X  =  A,  u  =  ny,  then 
n  n 


pn  =  p0 


0,Axn 


A  _ 

(i+l)y  nTV 


therefore 


-  e"X//y 
pn  -  ~nl  y 


p  =  1  -  pn  =  1  -  e“X/M 


A  busy  period  in  a  queueing  system  begins  when  a  customer  arrives 
to  an  enpty  system.  The  busy  period  continues  as  long  as  there  is  at 
lease  one  customer  in  the  system,  and  the  busy  period  ends  the  first  time 
that  a  customer  departs  leaving  behind  him  an  enpty  system.  For  the 
fVty'l  system  with  An  =  A,  y^  =  y  the  probability  density  of  the  length 
t  of  a  busy  period  is 

p(t)  =  i/t/p  e  +  I^(2ti/ay) 


where  again  p  =  —  and  I^(x)  is  the  modified  Bessel  function  of  the 
first  kind,  of  order  one  [19] .  The  average  length  of  the  busy  period  is 


simply 


1 

yd  -  p) 


For  an  excellent  treatment  of  queueing  theory  the  reader  is  direc¬ 
ted  to  the  book  by  Oox  and  Smith  [7] . 
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